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Abstract 



a: 

We investigate the issue of regularization/renormalization in the presence of a non- 
trivial background in the case of l+l-(supersymmetric) solitons. In particular we study 
and compare the commonly employed regularization methods (mode- energy/momentum- 
cutoff and derivative regularization). We show that the main point for a consistent reg- 
ularization/renormalization is to find a relation between the "cutoffs" in the vacuum and 
the nontrivial sector so that they can be related in a consistent manner. For each scheme 
we give a principle to derive this relation and to perform calculations in a consistent 
way. These principles are simple and not restricted to two dimensions or supersymmetric 
theories. 
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1 Introduction 



1.1 Historical 

In August 1834 the technician and marine-engineer John Scott Russell rode along the Edinburgh- 
Glasgow-channel and watched a boat which was pulled by two horses. Ten years later (1844) 
he described his observations in a report to the British Association for the Advancement of 
Science 0. He wrote that as the boat stopped a large single wave-amplitude with large ve- 
locity and constant shape was running down the channel. He followed the wave on his horse 
and lost it after two miles. This was his first encounter with this "singular and beautiful 
phenomena" which he called "wave of translation". This is the first known (at least to me) 
mention of what we today call "soliton" or a "solitary wave" and it was forgotten for more 
than a century. In 1895 Diederik Korteweg and his PhD-student Hendrik de Vries discovered 
a nonlinear equation, describing water waves, the so called KdV-equation. They showed that 
their equation posses a solitary wave as solution. But at this time it was seen as an accident 
that nonlinear equations could be solved explicitly and therefore their discovery was almost 
forgotten. 

In the twentieth century the invention of computers made it possible to investigate nonlin- 
ear problems with prospect of success. In 1965, seventy years after its discovery, Norman 
Zabusky and Martin Kruskal investigated the KdV-equation, which was found to describe 
different systems. In numerical calculations they discovered that a solitary wave can overtake 
a slower one and after a complicated, nonlinear interaction the two waves continue moving 
with unchanged velocity and shape. The residual effect of the interaction is a phase shift in 
the relative position of the two waves, an effect which is impossible for linear waves. Because 
of the individual character of these nonlinear waves Zabusky and Kruskal coined the notion 
soliton. 

After this work an intensive investigation of nonlinear soliton-bearing equations began and rich 
connections between different branches of physics and mathematics - scattering theory, lattice 
dynamic, Kac-Moody-algebras, Verma-modules, cohomology, topology, Pontrjagin numbers - 
were found. 



1.2 Current status of research 

During the last decade an enormous flurry of activity and also substantial progress has taken 
place in understanding non-perturbative effects in both supersymmetric field theories and su- 
perstring theories ^fl. Central to this is the occurrence of extended objects such as solitons and 
instantons ||, whose masses and actions, respectively, are inversely proportional to coupling 
constants so that they gain importance in the strongly coupled regime. As first observed in 
the two-dimensional sine-Gordon theory (0,[^,|26|), there is the possibility of an intriguing 



duality between the ordinary elementary quanta of quantum field theory and bound states of 
solitons. 

The consistent regularization and renormalization in the quantization of (supersymmetric) 
solitons is still an active area of research (see |33fl,[|$4]],||3o|] and references therein) with a 



number of not completely resolved issues. In some special cases certain properties of the 
quantum theory in the presence of a non-trivial background can be decided without doing 
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explicit calculations of the quantum corrections, such as the saturation of the BPS bound in 
two-dimensional minimal supersymmetric theories with kink-solitons (SUSY- sine-Gordon and 
4> A model) by the use of the residual supersymmetry [36|] or the quantum mass of the kink in the 
case of the minimal SUSY-sine-Gordon model obtained from Yang-Baxter equation assuming 
factorization of the S-matrix |37|] , [38|] . These methods provide highly non-trivial cross-checks, 



but depend on special properties of the considered theories and give no general insight into 
the impact of the presence of a nontrivial background on the renormalization procedure and 
thus different aspects of the Quantum Field Theory associated to nontrivial classical solutions. 
These aspects are overruled by "higher" knowledge such as supersymmetry. 

In this work we deliberately do not make use of supersymmetry at each step (which does 
not mean that we neglect or violate it unnecessarily). We mostly carry out the calculations 
for bosons and fermions separately to point out the different aspects of the influence of the 
nontrivial background on fermions and bosons. We investigate different regularization schemes 
(mode-, energy- cutoff and derivative regularization) within the renormalization procedure 
which are well known in standard perturbation theory and adapt it to the requirements of 
the presence of a non-trivial background. Besides the resolution of various subtleties we 
are able to show that the considered regularization schemes, which are very successful and 
popular in standard perturbation theory, with the necessary modifications are still well working 
techniques even in the presence of a non-trivial background. The investigations demonstrate 
that different schemes emphasize different aspects of the nontrivial background for quantum 
corrections but all of them eventually give the same unambiguous results. We thus are able 
to solve certain outstanding problems in the computation of the quantum mass of (SUSY) 
solitons. Over a long period of time different methods gave different answers ([[3£|, ||) an d 



there seemed to be no convergence in the results. More recent works |38|],[[3£j cleared up a lot 
of things in this discussion but also posed new questions, which are still in discussion f33|,[[|§ . 
Clearly, the resolution of the remaining open points is an important step for reliable further 
investigations. 

The modifications that have to be applied to the different regularization schemes are based on 
very simple principles which are not restricted to two dimensions or supersymmetric theories. 
The generality and simplicity of these principles thus pave the way for further investigations in 
more general cases than 2D-SUSY solitons. Nevertheless minimal 2D-SUSY solitons are still of 
particular interest. Firstly, the discussion on the consistent renormalization in the presence of 
exact static classical solutions is concentrated on these models, secondly there exists a "higher" 
knowledge due to supersymmetry or exact S-matrices which makes it possible to verify these 
principles "a posteriori" and above all because of their simplicity, so that one can focus on the 
problem of consistent regularization and renormalization in the presence of a static nontrivial 
background. 

The new principles are formulated and used for the calculation of the quantum correction to 
the soliton mass. When they are respected, all considered methods give the (one loop) mass 
corrections, now "accepted by all workers in this field" |34fl , (m = Vl/j,, where /i is the minimal 
renormalized mass parameter and Z = 1, 2 for SG, (p 4 ): 

nm 

(susy) - SG : AM B = 

7T 

AM F = 

2tt 

Tivn 

AM susy = AM F + AM B = -— 

Ztt 
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(susy) — 4 : AMb = hm 
AM F = hm 
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fl/TTl 

AM susy = AM F + AM B = -— 

In 

Although in the supersymmetric case the bosonic and fermionic corrections AM B , AM F have 
no physical meaning by themselves we have calculated them separately due to the reasons 
mentioned above. Only the sum AM susy has a physical meaning and that it is the same in 
both theories is related to its supersymmetric origin. 



1.3 Organization of this work and conventions 

In section || we review some properties of solitons. This section mostly follows reference [3|]. 
In section [3| we first discuss general principles of the quantization of (static) solitons and 
renormalization. This will be used to calculate the (one loop) quantum corrections to masses 
of the (p 4 - and sine-Gordon- kink-soliton solutions in section |j. The main point of this section 
is a consistent regularization in different "topological" sectors. In section |5| we consider solitons 
coupled to fermions. Especially the supersymmetric extensions of the bosonic theories of the 
foregoing sections are considered and the additional (one loop) fermionic quantum corrections 
to the kink masses are calculated. 

Throughout this work we use, except stated otherwise, units so that c = 1 and H ^ 1, because 
h will be our main perturbative parameter. The metric signature is (+, — ). 

2 Classical Solitons 
2.1 Introduction 

[|First we consider the massless Klein-Gordon equation in 1+1 dimensions 

n<j>{x,t) = (^%-%)<j>{x,t) = 
This equation and its solutions have well known properties 

* It is linear and dispersionless. 

* Each "well behaved" function of the form f(x± ct) is a solution. 

* It is a second order equation, and the plane waves cos(fcx + out) and sin(A; ± tot) 
with us = kc form a basis in the space of general solutions. 

lr This is the only section where we use units in which c ^ 1. 
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Thus each "well behaved" function can be expanded according to this basis: 

f(x — ct) = J dk[a(k) cos(kx — cut) + b{k) sin(foc — cut)). (1) 

For proper functions a(k) and b(k) this is a wave packet traveling in positive x -direction with 
the velocity c and since all modes have the same velocity the shape of the wave packet is 
stable, i.e. constant in time, and thus dispersionless. 

Because of the linearity of the massless Klein-Gordon equation a linear combination of solutions 
is again a solution. Thus one can construct a solution built of several wave packets which can 
travel with different (opposite) velocities. Consider, for instance, two wave packets: 

fs(x, t) = fx(x - ct) + f 2 (x + ct) 

This solution has following properties: 

t — > oo : Two widely separated wave packets. 
t = finite : Collision of the wave packets. 

t — > — oo : Again two widely separated wave packets with the same shapes and velocities as 
before the collision. 

Solutions with several wave packets have analogous properties. For the massless Klein-Gordon 
equation we conclude, that it is linear and dispersionless and from this follows: 

(i) Shape and velocity conservation of a wave packet. 

(ii) Asymptotic shape and velocity conservation after collision of several wave packets. 



2.1.1 Klein-Gordon equation in D=l+1 

The Klein-Gordon equation, 

(□ + m 2 c 2 )(f)(x,t) = 0, 
is also linear and a solution basis is again given by plane waves 

cos(kx ± ujt) and sin(fc±co>t). (2) 
But now for (Q) being a solution ui and k must fulfill the following equation: 

2 (22, 24 

u = k c + m c . 

From this follows that different modes (different k's in ([|)) of the wave packet move with 
different (phase) velocities 



.,. Lu(k) / m 2 c 2 

and a wave packet with a certain shape at time t — 0, 

f(x, 0) — J dk[a(k) cos(kx) + b(k) sin(A;2r)], 

will spread as time moves on. The Klein-Gordon-equation is dispersive. Thus solutions of 
the Klein-Gordon- equation have neither property (i) nor property (ii) of the solutions of the 
massless Klein-Gordon-equation. 
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2.1.2 Nonlinear equations 



If we neglect the dispersive mass term of the Klein-Gordon-equation and add a nonlinear term 
instead to get something new, for example 

D<f>(x,t) +(p\x,t) = 0, (3) 

we also get wave packet solutions which spread with time. For @ this can be observed by 
numerical calculations. 

In equations with dispersive and nonlinear terms, these two effects can balance each other, 
so that special solutions occur that have the property (i) or even the properties (i) and (ii). 
Solutions of nonlinear equations with property (i) are called solitary waves. Solutions of 
nonlinear equations with properties (i) and (ii) are called solitons. It is common to call both, 
solitons and solitary waves, solitons (or "lumps" ||). We will give a more precise definition 
in the next chapter. These non-dissipative solutions which do not spread out with time form 
lumps of energy holding themselves together by their own self-interaction. 



2.1.3 Solitons and solitary waves 

The definition of solitons and solitary waves vary from author to author, but the several 
definitions are very similar, they differ only for special cases in the characterization of solutions 
of nonlinear wave equations. We give a definition which is appropriate for our interests and 
equivalent to other definitions in the cases treated here. 

We characterize localized solution by the energy density (Hamiltonian density) of the field 
configuration (composite function of the fields) 

e(x,t) = F((f>i(x,t)). 

Connected to the energy density is the total energy (Hamiltonian) of the system by 

space 

For physical systems the energy is bounded below and we can shift the minimum to zero, i.e. 
E m in = . With this normalization we define localized solution as follows: 

Definition We call solutions of nonlinear wave equations localized solutions if the associated 
energy densities e(x, t) have following properties: 

e(x,t) is localized in space for finite times t, i.e. 

(i) ]im.e(x,t) — > fast enough to be integrable 

x— >oo 

(ii) e(x,t) is finite in finite regions of space 

For systems with E[<f>i] = iff <pi(x,t) = this definition of localized solutions is equivalent 
to localized fields, i.e. if one requires lim t) = = lim <9^0(x, t). 
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Definition A solitary wave is a non-singular localized solution of nonlinear field equations 
whose energy density has a space-time dependency of the form 

e(x, t) = e(x — ut) 

where u is an arbitrary velocity(vector). 

This means that the energy density moves with constant velocity and constant shape, i.e. 
undistorted. From this follows that any static localized solution is a solitary wave with u — 0. 
For relativistic or Galilean invariant systems one obtains moving solutions by boosting the 
static ones. Therefore static nontrivial solutions with localized and finite energy will be of 
central interest for us. 

Definition A solution of nonlinear field equations with N solitary waves with energy densities 
e (x — ut) is called a soliton if the energy density has following properties: 

N 

(i) e(x, t) — > ^2eo(x — di — u{t) as t — > — oo 

i=i 

N 

(ii) e(x, t) Yl £ o( x — — u{t + Si) as t — > +00 

i=i 

Where Oj and Ui are the initial positions and velocities and Si are constants or constant vectors 
for higher dimensions. 

So solitons are solitary waves whose velocities and shapes of energy densities are asymptotically 
(t — > +00 ) restored to their initial, i.e. pre-collision, ones. The constants (vectors) Si represent 
the displacement of the pre-collision trajectories and should be the sole residual effect of the 
collision. It is clear that all solitons are solitary waves but not vice versa. 

2.2 General properties of scalar solitons in D=l+1 

We give a qualitative discussion of possible solutions of nonlinear field equations for scalar 
fields. We consider the simplest cases, i.e. only one field in D — 1 + 1 dimensions. Of 
special interest are static solutions with localized energy, which are transformed by a boost 
into moving localized energy-lumps. This fits to the concept of a particle, but with finite 
extension. The quantum theory of these objects (see part 2) will validate the particle-picture 
of these extended objects. 

The dynamics of such a simple system is described by the Lagrangian (density), 

C=\{dci>y-U{ct>) = \i?- l -ct>i 2 -U{ct>) (4) 
for which we set up following assumptions 

(i) C is 2D-Lorentz-invariant 

(ii) U (0) is a positive semidefinite function of and does not depend on the derivatives 
d(f> of the field. 
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(iii) absolute minima of U are zero, i.e. U m i n = 

The equations of motion (e.o.m.) are obtained by a variation principle: 

dx 2 C = => n<f){x,t) + — {(f)) = 0. (5) 

In this case the variation of the fields and its derivatives vanish per definition at possible 
boundaries of the considered space-time interval (for an unbounded space-time l" this is 
always automatically true). This variation principle is consistent with the second order e.o.m. 
This is not true for first order systems like fermions, where one needs a modified variation 
principle (see below). 

The Hamiltonian (energy) of systems like (M is given by 

E[<j>] = f dxe(x,t) = J dx(^ -Q = J <M~0 2 + \<t>' 2 + U(4>)\ (6) 



where time independent Lagrangian's give time independent Hamiltonians,i.e. Q£ = =>- 

dE = Q 

dt u - 

From (H) one can see that the energy E[<j>] is the sum of positive definite terms. Thus to get 
the minimum of the energy, i.e. the ground state (the vacuum in the quantized theory) of 
the theory, each term must be minimal. This is achieved by constant fields (f)(x, t) = for 
which the potential U((f>) has an absolute minimum, i.e. U(g^) = U m in = because of our 
assumptions. 

2.2.1 Boundary conditions 

The requirement of localized and finite energy solutions implies strong restrictions on the 
spatial boundary conditions for these fields. For a localized energy density e(x,t) i.e. for 
finite energy E[4>] the energy density e(x,t) must vanish (fast enough) with |x| — > oo for the 
considered field-configuration. From (|6|) one can see that this is only possible if asymptotically 
(\x\ —> oo) d^(p = and U((p) = 0, since all quantities in e(x,t) and E[4>], respectively, are 
positive. 

Let 4> — 9^\i — 1-..M > 1 be the absolute minima of potential, i.e. U(4> = grW) = 0. 
A necessary condition for e(x,t) to vanish asymptotically is that cf> approaches one of the 
(constant) minima . Thus the localized energy condition e(x — > ±oo) — > implies for the 
field the following solitary wave-boundary conditions: 

lim 0(x,t)=/ ± ) (7) 
lim dMx, t) = (8) 

x— >±oo 

One has to distinguish between the case of one minimum (M = 1) and the case of several 
minima (M > 1). If M = 1 then for both limits (x — > +oo and x — > — oo) the field converges 
to the same value g( % ^ = g. For M > 1 one can have different limits g l+ ^ g l or the equal 
ones, g l+ = g l . 
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\ V(q) 



Figure 1: The potential U(tp) (solid) of the field theory and its (one degree of freedom)-mechanical 
analogue V(q) (dashed) in the case of static fields in D = 1 + 1 dimensions, (a) For a field theory 
with a unique ground state and (b) for degenerated absolute minima of the potential U(<j>). 



2.2.2 Static solutions and the mechanical analogue 

For static solutions the e.o.m. (||) simplify, in our D — 1 + 1 case, to an ordinary differential 
equation: 



dU 

D(f)(x,t)\ s tatic = ~<f> fa) = — ^jl 



X 



(9) 



This equation is analogous to the Newton equation of motion of a unit mass in a potential 
V = —U (see fig.|l]), if one considers as the coordinate and x as the time, i.e. <f>(x) = q(t). 
Therefore we discuss the familiar mechanical analogue instead of the original field system and 
solve the e.o.m. (||) simply by quadrature. This is only possible in D = 1 + 1 dimensions . At 
the end of this section we make some comments on the situation in more general cases. The 
mechanical analogue has the following properties: 



Energy 
Boundary conditions 

Action 



W 



2 

lim W 



2 K dx- 



+±oo 



S[q] = 



+ V(q) 
= 



U((j)) = constant 



E\ 



\ \static 



(10) 

(11) 

(12) 



The boundary condition (|TT]) follows from (|7|) and (§), (0) f°ll° ws from 
Virial theorem. Integrating @ with J^ L dx(f)'(x) gives: 



1 ,2 



u 



(13) 



This is the virial theorem for the mechanical analogue and in connection with non-trivial static 
field configurations it is called the Bogomol'nyi equation. 

Thus static solutions of the field system are trajectories with finite action, (0), and zero 
energy, fllPf) and ([□]), of the mechanical analogue. We consider these trajectories for two 
classes of potentials: 

a) From fig.|l] one can see that for potentials U with a unique minimum there exists no non- 
singular, nontrivial trajectory with boundary condition (|TT|). A particle starting at the "time" 
x = — oo at 0i = will never return. The only solution is the trivial one <p(x) = 0i, i.e. there 
exists no static solitary wave (the field is constant in space). 

b) In the case of several degenerated minima of U according to (|Tl|) the particle must start at 
one of the minima of U, <pi, and move to one of the neighboring minima 0j±i(the field varies in 
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space). It cannot return or go further since all derivatives of 0, i.e. "velocities", "accelerations", 
etc, vanish at the fa's due to the equation of motion (||) and (|i~3|) and boundary the conditions 

u(fa) = o => fa' = o 

dU 

— ((pi) = = 0/' a.s.o. 

From the mechanical analogue one concludes for the existence of static solitons for theories of 
the form C= |(^0) 2 - U{4>) (one field): 

1. U(4>) has an unique minimum =^ 3 one trivial static solution fa and $ static solitary 
waves. 

2. U(fa) has n minima =>- 3 2{n — 1) nontrivial static solutions with (p(x — > — oo) = 0fc and 
4>(x — > oo) = or 0fc_i, and n trivial solutions fa. 

From the mechanical analogue one can also see that the "particle" moves monotonically from 
one minimum oil! (top of the hill) to a neighboring one. Therefore the static solitary wave is a 
monotonicly increasing or decreasing function. The above considerations are not restricted to 
a special shape of the potential U(fa). The main point is the existence of several (at least more 
than one) degenerated ground states (absolute minima of U(fa)) which can be accompanied 
by spontaneous breakdown of a symmetry (see below). 

Solving by quadrature. As mentioned above, in the simple D = 1 + 1 case static solutions 
are obtained by quadrature. Integration of the virial theorem gives 

f m deb , x 

x - x = ± / y (14) 

Because of the boundary conditions the integrand is regular except for x - > — °o and x — > oo, 
where Xq is the integration constant. 

Next we consider two special models and investigate further features of solitons (solitary waves) 
on the basis of them. These two models will also be considered in the quantization procedure 
(section |3|). 



2.3 (J) 4 - theory, the kink 

We consider the 4 - theory in D = 1 + 1 dimensions (not dimensionally reduced) with a 
mass (quadratic) term which causes "spontaneous symmetry breaking". The 4 -model with 
the opposite sign of the mass term generates only one unique minimum so that there exist 
no static solitary solutions, as mentioned above. The Lagrangian, which also fulfills our 
assumptions ( ^2|) is given by 

c=\{d^y-u{<i>) , (15) 

U{fa) = f(0 2 - ^) 2 - (16) 
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The potential (0) has the shape of the potential (b) in fig.[l|. The e.o.m. of this system are 

□0(x, t) - [i 2 <\) + A0 3 = 0, (17) 

and for static solutions they read 

-f - fi 2 4> + A0 3 = o. (18) 

The minima of the potential have the value zero for 0±, where the ground state configurations 
= <f>± of the system are 

From our solitary waves- (localized energy) boundary conditions (|7|) follows that the field must 
asymptotically approach these values, i.e. 0(x — > ±oo) — > <fi±. Thus we have two possible 
non-trivial static solutions. One evolving with x from 0_ to 0+ and a second one in the 
opposite direction. 



Integration, static localized solutions With the potential (|i"6|) the general integral flTil) 
reads (er = ±1) 



x — x = —a 



cf>(x) 



Hxo) v^/2(0 2 - /i 2 A) 



This is a standard integral. For fields satisfying | ^0 |< 1 and by setting the integration 
constant 0(x o ) = 0, one obtains by elementary integration and solving for 0: 

<PkA x ) = ^-^= tanh[— (a; - x )}. (20) 

For later considerations we have introduced the sign variable a. One can easily prove by 
inserting in (111) that this is a solution. For fields satisfying | —0 |> 1 one gets functions 



involving coth. These cannot satisfy the boundary conditions. Thus we have two nontrivial 
static solutions which are called kink (cr = +) and antikink (a = — ), and shown in fig.|[ 
These two solutions can be traced back to the two signs of the Bogomol'nyi equation (|13|): 

0V = -<ry/2U(fa). 

As one can see, these solutions are singular for A — > 0. Thus they cannot be obtained by 
perturbation theory starting from the linear equations (A = 0). Thus the kink 0^ + and the 
antikink 0x_ are non-perturbative results. 



2.3.1 Spontaneous symmetry breaking 

We give some comments on the symmetry and its spontaneous breakdown of this model. 
Spontaneous symmetry breaking will be of importance in the quantization procedure but it 
is a feature which is already present at the classical level and does not come from quantiza- 
tion. In (quantum) perturbation theory one generally expands the Lagrangian around a fixed 
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(classical) field-configuration (p c i (see below). Thus the Lagrangian becomes a function of the 
perturbations 77 around the fixed configuration C /, i.e.0 

C(<P) = C(<p d + r 1 ) = C(r 1 ). 

Now assume that the theory has a linear symmetry, so that the Lagrangian is invariant (up 
to total divergences) under the linear transformations T0 of the field, i.e 

C(T<j>) = £(<f>) = £(<p cl + v ) = C( V ). (21) 

On the other hand, if the "ground state" <fi c i is not invariant under this transformation (if the 
boundary conditions do not respect the symmetry), i.e T(j) c i = <j) c i 7^ (f) c i, one has 

£(T0) = C{T<p cl + T V ) = £(0 C/ + T V ) ^ C(T V ) =}► C(T V ) ^ £{ v ). 

This effect is called spontaneous symmetry breaking. Actually it is just hidden symmetry, since 
by writing the Lagrangian as a function of the perturbations r\ the symmetry of the system 
is no longer visible although still present as one can see in (|2l|). The effect of spontaneous 
symmetry breaking gives rise to a rich structure in quantum field theory and particle physics 
(Higgs effect, Goldstone theorem). In the quantum theory the "ground states" C / are usually 
one of the minima of the potential, i.e. the configuration with the lowest energy, and the 
associated quantum mechanical state is the vacuum |0). Therefore in quantum theory on calls 
a symmetry spontaneously broken if the vacuum state is not annihilated by the symmetry 
transformation, i.e. 

T|0)^0. 

A less trivial "ground state" is a non-trivial classical solution like our kinks. Their quantum 
theory will be the main part of this work. 

Let us return to the 4 model. The Lagrangian ( |l~5|) and the associated action is symmetric 
under the parity transformations x — > Px = —x and separately for the Z 2 (gauge) transfor- 
mation — > Z(j) = —(p. By the Z 2 (gauge) transformation the two minima 0± fll9D are not 
invariant but transformed into each other, i.e. Z(fi± = <fi T . This is the "classical" situation 
of spontaneous symmetry breaking. But from (|20|) one can see that for the kink and antikink 
4>K„ both symmetries transform the two solutions into each other: 

parity : <PkA Px ) = <f>K.(-x) = -<j>K a { x ) = <Pk_A x ) 
^2 : Z(j) Ka (x) = -(f>K„( x ) = 0K-A x ) 
2 AW equations resp. inequalities are to be understood modulo total divergences. 
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Figure 3: The energy density for the kink and antikink. It is localized around the center xq of the 
(anti)kink an characterized by its half- width di = j-. 



This is typical for systems with spontaneous symmetry breaking, since the localized energy 
solutions connect the different vacua. That the (anti)kink is not invariant also under parity 
transformation (antisymmetric), will have interesting consequences for the fermionic boundary 
conditions in the quantum theory of supersymmetric solitons. 



Energy density, classical kink mass By the use of the virial theorem ((L3|) and the 
(anti) kink-solution one obtains for the energy density 

which is regular for real x and satisfies the conditions of the definition for solitary waves (see 
fig HD- So the kink and antikink are solitary waves. In analogy to the rest mass of a particle 
the total energy E[<f)] for the (anti)kink is called the classical kink mass, since they are static. 
It is given by 



2.3.2 Moving kinks 

By a 2D-Lorentz boost to a system moving with velocity u one obtains a moving (anti) kink. 
Since is a scalar field we only have to transform the coordinate: 

X — > j{x — ut) 



In these new coordinates the kink reads as follows 



0»(s,t) = -^tanh[4 > J^ Ut -) , u e (-1,1) (23) 



v^A ' V2 vT 



u 



It is easy to check that this is a solution of the equation of motion (|17D. The energy density 
and kink mass are obtained by inserting (p3|) into @: 
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For the total energy this gives with a variable substitution x — > ^~/f^=T^ i n the integral of 

] = j^/M— f dx[- 1 sinh2xi 1 2 ^ 3 Md 



2 



This is a very nice result, since it is the relativistic energy-mass-relation for a particle. So one 
can expect that the quantized theory will give one particle states associated with the classical 
extended object. 



2.4 Topological indices and topological conservation laws 

It is often possible to make a topological classification of the solutions of a system. Specifically 
if one can define a topological index which is conserved in time it will play an important role 
of a quantum number in the quantum theory as other conserved quantities. But its origin 
is different from that of familiar conserved quantities. Again we consider a special class of 
theories in 2D, namely 

where U(<f>) has a discrete set of absolute minima U{<p) m i n = 0. We are interested in nonsingular 
solutions with finite energy which also include solitary waves and solitons. The requirement 
of finite energy at any time t implies the following boundary conditions: 

lim 00, t ) = <K±oo, t ) = <f>± (24) 

x— >±oo 

where at <p± the potential U has an absolute minimum, i.e. U{4>±) — 0. Since this must 
be true for all t (all terms in E[<p] are positive, see (||)) and the discrete minima of U are 
separated, the field 4>(oo,t) = <f>± must be stationary, i.e. conserved: 

d t (j>{±oo,t) = 0. (25) 

Thus we can divide the space of all non-singular finite-energy solutions into sectors, charac- 
terized by two time independent indices, namely </>(— oo) and 0(oo). These sectors are not 
topological connected. Fields of one sector cannot be deformed continuously into fields of 
another sector without violating the finite-energy-condition. This emphasizes the difference 
between the conserved indices (|25|) and familiar conservation laws which are a consequences 
of continuous symmetries of the theory. The fact that different sectors are not connected is 
a consequence of the topological property of the space of non-singular finite energy solutions. 
For this reason (|25|) is called a topological conservation law. 

One can show that the existence of a topological conservation law is sufficient for the existence 
of non-dissipative solutions. This is important in more complicated theories, for which the 
direct integration is not so easy as for a single scalar field in D = 1 + 1 dimensions (this will 
be shortly discussed in section |2.7| ) . By means of the 4 model we want give an idea how this 
works. Instead of nonsingular solutions of finite energy we consider non-singular initial-value 
data 4>(x, t ) and d t <p(x, to) at some fixed time to (for the existence of nonsingular solutions for 
this initial value problem we refer to the reference in |§). For these initial- value data, just as 
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for the time-independent solutions (the kinks) the finite energy condition implies the relations 
(IHD and (|25|) . If U has more than one absolute minimum equation (p5|) is non-trivial. Now 
one can show || that any solution with the conserved indices 

</>(oo,t) = —</>(—oo,t), 

is non-dissipative, i.e. that the energy density does not spread indefinitely with time ( limmax e(x, 

t— »00 X 

0). By continuity in x, for any t, there must be some Xq for which <fi(xo, t) = 0. At this point 
the energy density @ is 

e(x o ,t)>U{0) = £ 
Thus for all times the maximum of the energy density is unequal zero, 

max £(x, t) > —, 
x ~ 4A 

and therefore the energy density does not dissipate but stays localized. In an analogous way 
the existence of nontrivial topological conservation laws can be used to prove the existence of 
non-dissipative solutions. 



Topological indices of the (anti)kink 

The potential U(4>) = 1/4A(0 2 — /U 2 /A) 2 has two minima at (j) m in = ±/x/\/A . This gives rise 
to four topological sectors of non-singular finite-energy solutions with the following indices set 
(writing only the signs) 

{(0(oc), 0(-oo))} = {(-, +), (+, -), (-, -), (+, +)}. 

The kink, the antikink and the two trivial solutions = ±tt| are elements of the four sec- 
tors respectively. Another example is a kink from x — > — oo and a antikink from x — > oo 
approaching each other. This field configuration lies in the (trivial) (— — 77=) sector. Even 
though one cannot easily calculate the collision, we know that the resulting field will always 
stay in that sector. In fact the (anti)kink is only a solitary wave and not a soliton. These 
topological constraints also stabilize the (anti)kink and because of this these nontrivial solu- 
tions will become fundamental particles in the quantum theory, since they cannot decay. For 
a decay the (anti)kink would be deformed into a trivial topological sector, which would need 
an infinite amount of energy. For the existence of such topological solutions it is necessary 
that the quantum corrections do not lift the degeneration of neighbouring ground states and 
thus spoil the requirements for theire existence. This effect usually occurs in the case of an 
accidental degeneration, i.e. if the minima of the potential are not related to a spontaneously 
broken symmetry. Since the classical theory is just a limit of the quantized theory no such 
solutions exists in that case [|J. 

This fact, that the (anti) kink cannot be deformed continuously into the trivial sector without 
violating the finite energy condition is the origin of misunderstandings in the use of boundary 
conditions in the trivial and nontrivial sector during the quantization procedure. Also because 
of this one temporarily uses the kink-antikink configuration for the quantization procedure 
which is an intractable trick in more complicated cases 

( HI , lll.il)- This will be clarified 

up later. 
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2.4.1 Topological charge and conserved current 

Although the conserved topological indices come from the finite-energy-condition and not 
from a continuous symmetry, one can define a conserved current and a corresponding charge 
connected to the topological indices 

Q-=c [((f)(x = 00)- ((f)(x = -00)] , k" := c e^d v (f> (26) 

where is the antisymmetric epsilon symbol and c is an arbitrary constant. This is trivial 
in D=l+1. With these definitions one has 

a^Jfe" = and Q = J dxk . (27) 

Note that plane waves e lP/tXl> do not change the topological charge. Assume a field configuration 
4>to P with a definite topological charge Qt op , like the kink or the vacuum. Then the topological 
current k p in (|2~6"|) and the topological charge Q (|2~7|) get an additional contribution from the 
plane wave for the field fa op + as follows 

5k" = e^d u (e ip ^) = e pu i Pu (e ip ^) (28) 
5Q = i Pl e tpot J dx e~ ipiX = Pl S( Pl )2m e ipot = 0. (29) 

Thus small (quantum) fluctuations will not change the topological charge and the topological 
sector of the classical finite energy field configuration with a certain topological charge. 

To classify the topological sector one needs <p(x = 00) and fax = —00), so that the knowledge 
of Q is not enough, but for quantities which depend only on the difference of the conserved 
indices Q is sufficient. 

For our 4 -theory we set c = VX/m so that Q G {—1, 0, 1}. Solitary waves are called topological 
if Q ^ , otherwise non-topological. This terminology says that the nontrivial solutions of the 
4 -theory are topological. 

Symmetry breaking and topological indices. Suppose the Lagrangian C is invariant 
under some transformation T of the fields. Then one can distinguish two cases: 

1. U has a unique minimum at 0o =^ 0o itself must be invariant under T, i.e. T(f) = <fi , 
since a symmetry transformation does not change the energy of a solution. 

2. U has several degenerated minima at fa , i — 1, ...M > 1 =^ the full set {0j} must be 
invariant under T, i.e. Tfa G {fa}, but not each fa itself. If not each fa itself is invariant 
under T one calls this a spontaneously broken symmetry (see above). 

In order to get non-trivial topological sectors the existence of more than one degenerated 
minimum of U is necessary and sufficient (see ref. in [||] for the existence of nonsingular 
solutions for nonsingular initial- value data of finite energy). Thus a spontaneous symmetry 
breaking gives rise to nontrivial topological sectors. The converse is not always true, since the 
fa's could nevertheless be invariant under T. 
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2.5 The sine - Gordon system in D=l+1 



The 4 - theory discussed above yields only solitary waves, but not solitons. The sine - Gordon 
system also yields solitons, as we will see. The sine-Gordon equation has a long story. In the 
last quarter of the nineteenth century it was extensively studied by geometers since it describes 



a two dimensional Riemann-surface of constant negative Gaussian curvature K = — 1 [10|| . It 
entered particle physics through works of Skyrme (1958,1960) who studied simple nonlinear 
field theories. Its name is a pun and seems to belonging to either Finkelstein and Rubinstein 
(Klein-Gordon — > sine-Gordon) or Kruskal who investigated numerical solutions of nonlinear 
field equations and also discovered solitonic solutions of the Korteweg-de Vries equation (see 



1T| and references therein. The history of this name is not completely clear but it has prevailed 
over other names). A mechanical system which is also described by the sine-Gordon equation 
is realized by a continuous chain of elastic connected pendular on a horizontal line in a constant 
gravity field (or an infinite ribbon with a load at one edge [|ll|). The field <f>(x,t) in this case 
describes the angular amplitude of the pendulum. The sine-Gordon equation also describes an 
infinite Jefferson contact |12j. As one can see a lot of systems are described by the sine-Gordon 
equation. 

The Lagrangian of the sine-Gordon system is given by 

£ = i(<9,0) 2 + ^ [cosCvW - 1] . 
A series expansion shows that approximations of this Lagrangian are well known systems 



For 7 — > this is just the free Klein - Gordon field and including the 0(7) - term one has the 
(attractive) 4 - theory. The equations of motions are 

□0+-^sin(VT0) = O. 



With a change in variables x a — > jx a and the rescaling — > ^=0 the Lagrangian and equation 
of motion writes as 

C= ^[-(<90) 2 + cos0- 1] (30) 
72 

□0 + sin = 0. (31) 

2 

As one can see, in principle the system can be described by only one parameter (3 :— which 
does not enter the classical e.o.m. It is a generic property that the classical field equations are 
independent of the coupling 7. Also for the 4 theory (|15D this can be achieved by rescaling 



the field as = A0. This can also be seen by the fact that in classical physics 7 (A) is a 
dimensionful parameter and thus can be used to set the scale. Of course, 7 (A) is relevant in 
quantum physics, since quantum physics contains a new constant, H, and the important object 
in quantum theory is ^£ and the relevant dimensionless parameter is ^ or hX, respectively 
(see section. |). 

Energy and energy density are given by (note that dx — > dx) 

E[(f>] = J dxe(x,t) = J dx^cj) 2 + ^0' 2 + (1 - cos0)] (32) 
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The Lagrangian C and so the equations of motion have the following discrete symmetries 



-> -0 and -> + 2vriV, N E Z (33) 
Since the overall factor /i 2 /7 does not enter the classical e.o.m. (|31~1) we consider the potential 

[/(0) = (1 -COS0), 
whose minima are given by (note that U > 0) 

^ = s in0 = M = 2vrM, MgZ 

Thus we have a countably infinite set of absolute minima for which the energy E[4>m] vanishes, 
since U(</>m) = 0. The minima 0m are transformed into each other by the discrete symmetries 
(|33|) and therefore this symmetry, except for a Z 2 transformation, is also spontaneously broken. 



For finite energy - solutions, i.e. lim e(x,t) = 0, we get the following boundary condition 

x— >±oo 

0(x -> ±oo) = 2nM± (34) 

Thus, according to (|3|) we can characterize our topological sectors by the conserved pair of 
integer indices (M+,M_). If only fields modulo 2n are physically relevant (this depends on 
the "interpretation" of the particle states ||), then only the topological charge 

q = M+ - M_ = 1. [ dxd x <j) 

111 J 

matters. 

From our analysis of the mechanical analogue we know that a nontrivial finite-energy solution 
must move with x from one absolute minimum of U to a neighboring one, i.e. they must carry 
the charge Q — ±1. From (|i~4|) we get the static solution as follows: 



X — Xq = ± 



4>(x) 



B0 ) a/2(1 - COS0) 



with |(1 — cos?/) = sin 2 | and setting the integration constant 0(x o ) = tt one obtains by 
solving fotfjff) 

5± ( x ) = 4 arctan[e ±(a; - Xo) ] + 4A;tt , fcGZ (35) 

The solution 0s + with the (+) sign is called the soliton, the solution 0s_ with (— ) sign is 
called the antisoliton of the system. Their graphs are plotted in fig.| and are very similar to 
the kink and antikink of the 4 model (for both models we will often call them simply kinks 
4>k)- As one can see, for both, 0s and 0s_, there exists an infinite set of solutions which 
connect different neighboring minima 0m of the potential. Their topological charges are 

Qs + = t^-[0s + (z = oo) - 5+ (x = -oo)] = (2k + l)-2k = l 
Qs- = ^-[<f>S-(x = oo) - s _(a: = -oo)] = 2k - (2k + 1) = -1. 

Z7T 
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Figure 4: The "fundamental" (k = 1) soliton 4>s+ an d antisoliton 0s_of the sine-Gordon system. 



For the mechanical realization this solution describes a chain of pendular which is turned 
around once at the position Xq. The sign of the topological charge describes the orientation of 
the winding and is opposite for soliton and antisoliton. one obtains the classical mass of each 



(anti) soliton using the virial theorem (|i~3|) and ([32 



I dxe(x) ]static = ^ I dx[h' 2 + U\ = ^ / dx—jl = ^ (36) 

J 1 J 2 7 J cosh z (x - x ) 7 



2.5.1 Time dependent solitons 

Again one obtains moving solutions by a 2D Lorentz-boost, i.e. a coordinate transformation 
(the 7 used here has nothing to do with coupling, it is the usual parameter of relativistic 
kinematics): 

x — ut 

x — > j{x — ut) 



From this we get the moving (anti) soliton 

<f>s±,u(x,t) = 4arctan[e ±7(x - a:o - ut) ] + Akn , k G Z. 

We already called these solutions (anti) solitons because the sine - Gordon system provides 
several soliton-solutions according to our definition of solitons. We only can mention some 
examples here: 

Soliton-antisoliton-scattering. This solution of (pl| ) is 



i / \ A r MCOsh(7x), r smh(7u£) , . n 

<PsA{x,t) = 4arctan — = 4arctan — — - ± In 

sinh(7«t) u cosh(7a;) 

where 7 is defined above and the sign in the second form depends on the quadrant of the 
argument of arctan. That this is an exact solution can be proved by insertion in the e.o.m. 
The asymptotic behavior shows that this is a soliton solution (with the abbr. 5 — " ' 1 



4>SA(x,t) — > 4 arctan [2m cosh(7x)e 7 "*] 
t— *— 00 



4 arctan [ e 7(*+«(t+*))] = 5+ ( 7 [ x + u (t + §)]) 
4 arctan [ e -70=-«(t+*))] = 5 _ ( 7 [ x - u (t + 5)}) 



3 In the "old" literature one uses different branches of the inverse tangens. We follow || and use the unique 
arctan to avoid misunderstandings. 
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> X 



Figure 5: Soliton-antisoliton scattering at the times t « 0, t 0_ and t » 



Thus for t — ► — oo we have a soliton moving with velocity — u from x — ► oo, i.e. to the center 
x — 0, and a antisoliton moving with velocity m from x — > — oo, i.e. an soliton and antisoliton 
approaching each other. For the asymptotic future we have 



<Psa{x, t) — ♦ 4 arctan^-f^-v] + 2vr 



2« cosh(7x) J 



— ► 4arctan[e- 7 ( x -"(*- ,5 »] + 2tt = 5 „(7[x - u(t - 5)}) 

x^oo 

— * 4 arctan[e^ +u ( f - 5 »] + 2tt = 5+ ( 7 [x + u(t - 6)]) 

x—>—oo 

which are again a soliton and antisoliton with the same shape and velocities, but now departing 
from each other. The only change from the initial condition is the time delay 5 which remains 
the sole residual effect of the collision. Since in our units u < 1, the delay 5 is negative. This 
indicates that the soliton and antisoliton attract each other. This can be seen very illustrative 
by the mechanical realization, where the two different windings attract each other. 

At t — 0, 4>sa vanishes (the two opposite windings come together and unwind each other), 
i.e. the approaching (anti)solitons tend to annihilate each other until t — 0, but the field 
re-emerges for growing t and asymptotically restore the soliton and antisoliton (fig.||). Since 
every solution of the sine-Gordon system is a solution modulo 2n the graph can be shifted up 
and down by 27r-steps. 

Soliton-soliton scattering. This solution is given by 

r «sinh(7:c). 



(pss(x,t) = 4arctan[ 



cosh (71^) 



An analogous procedure as above shows that asymptotically two solitons approach each other 
for t — > —00 and departing from each other with same shape and speed, but with opposite 
velocity and a time delay. Thus they bounce back (backward scattering). At any instant of 
time the field ranges from — 2n to +2tt as x goes from —00 to 00. So the solution lies in 
the Q = 2 -sector (= total winding number of the pendular-chain). If we do not distinguish 
between fields modulo 27r, then there is no difference asymptotically between backward and 
forward scattering (see fig|]). The discrete symmetry under — > —cj) gives us an analogous 
solution for two antisolitons, namely <\>aa — —§ss- 

The doublet or breather solution. This solution is obtained by setting u = iv in the 
solution (psA (omitting the modulo-constant 2ti) 

sin( vt ^ 



4> v (x,t) = 4arctan[ ; ] (37) 

^cosh( 7 ^ f ) 
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Figure 6: The soliton-soliton scattering at t < 0, t = 0, i.e. their maximum approach and t > 0. 
After this they disappear from each other again in the opposite direction of their approach. 




Figure 7: The famous sine-Gordon breather for different times ii > t2 > £3 > £4. It is periodic in 
time with the period r. 

this is still a real exact solution of sine-Gordon system. The parameter v is now not a velocity 
but it is connected to the period r of this periodic solution 

T = Z7C . 

V 

The breather solution can be considered as a soliton and antisoliton oscillating about each 
other (see fig.[7|). In contrast to the scattering solution (f>sA the soliton and antisoliton input 
does not separate arbitrarily far apart as t — > ±00 but rather separate only up to a finite 
distance and never become fully free and undistorted. Thus this can be seen as a bound 
solution. The solution is given in its rest frame, i.e. it is centered around x = for all time. 
Moving breather-solutions are again obtained by a 2D-Lorentz boost. But also in its rest 
frame the breather has nontrivial time dependence in contrast to the former solutions which 
are static in their rest frame. The breather does not fit into our definition of solitary waves 
but the field and the energy density is localized. 



2.6 Stability and zero modes 
2.6.1 Static non-trivial solutions 

We have already mentioned that the static finite energy solutions will become new fundamental 
particles(states) in the quantum theory. Therefore we investigate the stability of these objects 
to see if they "survive" the quantization procedure. We have already seen that the existence of 
a topological conservation law leads to non-dissipative solutions, i.e. solutions for which the 
energy density does not dissipate to zero. Now we investigate the solutions explicitly according 
to their behavior under small perturbations. The equations of motion for our models are 

□0 + L7 / (0) = O. (38) 
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Next we consider a small perturbation of a static solution c /(x) of this e.o.m. 

(p(x,t) = (p d (x) + r](x,t), (39) 

and investigate the development with time of the perturbation T](x,t), which is determined by 
e.o.m. Inserting (|39|) into (|38|) one gets the following e.o.m for the perturbation rj: 

U V + U"[<p d (x)}r ] = Q + 0(r ] 2 ). (40) 

The linearized e.o.m. for the perturbation is separable for static solutions (invariant under 
time translations) and the general solution for rj is a superposition of normal modes: 

V (x,t) = ReJ2ane iu;n %(x) (41) 

The sum has to be treated as an integral for continuum modes. The coefficients a n are arbitrary 
complex numbers and the frequencies u) n and modes £ n have to solve the eigenvalue problem 
(inserting the ansatz fl4l|) into the linearized equation (|40|)): 

{-d 2 x + U"(<p d )) £ n = (42) 

This is a Schrodinger equation with the potential U"((j> c i) and in this context called stability 
equation. If one can explicitly solve this equation, as it is possible for the SG - and 4 - model, 
one obtains (in the linear approximation) a rich set of solutions of the e.o.m. In (f|l]) one can 
make the following classification for the perturbation w.r.t. the eigen- values u 2 of (ff2"|): 



u 2 > : The mode £ n stays oscillatory in time. 

u> 2 < : The mode £ n grows exponentially with time and thus also the perturbation, the 
linearized equation is only valid for short times. 

on 2 . = : The mode £ n is a so called zero mode, the perturbation is constant in time. 



First of all we can see that the static solution (fi c i(%) can only be stable if all eigenvalues of 
the stability equation (|4^) are positive, i.e. u n > for all n. The occurrence of zero modes is 
connected with symmetries of the system. In our special case it is the translation symmetry 
of the 4 and the SG- kinksQ. The explicit solution of the stability equation of the 4 and 
the SG- kink are given in the appendix (section [74] ) . In general one can show that because 
of the spatial translation invariance all eigenvalues are positive and therefore the solution 4> c i 
is stable ||. B ecause of the spatial translation invariance, also the boundary conditions (0J|) 
are translational invariant; also 4> c i{x + Xq) is a solution if (j>d{ x ) is one - F° r infinitesimal 
translations one has 

4> d (x + Xq) = 4>d(x) + Xod x <pd(x) + O(Xq), 



inserting this into the e.o.m. (|38|) one gets as zero mode of the stability equation (|42|) , i.e. 
the eigen-function to the eigen- value u 2 = 0, the function d x <p c i(x). As mentioned above the 
static finite energy solutions c z, connecting neighboring minima of the potential U (</>), are 
monotonic functions. Therefore the derivative d x (p c i has no nodes. It is well known that for 
one-dimensional Schrodinger equations with arbitrary potentials the eigen-function with no 



4 Here we use as mentioned above for both models the notion kink for the static finite energy solutions. 
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nodes is the eigen-function with the lowest energy. Since this eigenvalue is equal to zero, all 
eigenvalues u 2 are positive, and thus the perturbation rj stays oscillatory in time and the static 
finite energy solutions are stable. 

The occurrence of the zero mode will cause troubles in the quantization procedure. Another 
way of writing the stability-operator in (|42| ) and (|40|) is to express it through the action 
S[4>) = J dtdxC. For the classical solution (p c i(x) the action is stationary, thus the e.o.m. 
writes as: 

SS, 



0. 



(43) 



Expanding the action around c j one obtains^ 

S[<j )cl +7 1 } = S(<f )d ) + 



1 f 


~5 2 S 


+ - / dtdx rj 


_g ( f?\*<*_ 



rj + 0{rf). 



(44) 



The linear term is absent because of the e.o.m. ([43]) and the spatial part of the operator 
f^f l^is exactly the operator in ([42]). The second term in flli] ) will be the central object in 
the quantization procedure. 

To see that the occurrence of zero modes is in general connected with symmetries we consider 
a very general theory with an arbitrary number of fields {4>i}'. 



(45) 



For special constant matrices A^J ,3^ and potentials U(4>i) this theory also includes fermions 
and gauge fields . The equations of motion are given as 

9U{(j) h ) x 



V 



+ 



0. 



(46) 



where XV, is the matrix- valued differential operator, including first and/or second order deriva- 



tives in general, of the coupled system (|46l) . For fermions (Grassmann fields) one has to take 
care of signs of course. Assume that the field "vector" {<f>f (x)} = 4> d is a static solution of the 
e.o.m. ( |46"|) for certain boundary conditions, i.e. 



with 



(±oo) = a 



where C± is some constant vector in the field-space. Now assume the existence of a continuous 
(internal or space-time) symmetry 1Z of the theory ( ^5| ) which does not involve the time and 
leaves the boundary conditions invariant. Since under a symmetry transformation the e.o.m. 
are invariant, the field configuration 0, = TZ(pf is also a static solution. Since 71 is a continuous 
symmetry it is also possible to consider infinitesimal transformations. By an infinitesimal 
transformation one obtains a static solution of the form 



n4 



id 



()■ 



TZ' 



Since (pi is a classical solution, the deviation Sn't'k fulfills the full stability equation (|40|), and 
since 0j and (j)f is static also the deviation 5^0^ must be static and therefore the frequency of 



'For the moment we neglect subtleties connected with surface terms. For the exact treatment see section 
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this mode in (|4l| ) must be zero. Thus Snfak xs a zero m °de of the Schrodinger equation (|42|) . 
The invariance of the boundary condition C± ensures that fa is also a finite energy solution. 
This restriction is not really needed since because of the topological conservation law it is 
impossible for a continuous (symmetry) transformation to change the topological sector of a 
field configuration. If there are more, independent symmetries present which fulfill the above 
requirements, then each of them has its own zero mode. We call two symmetries independent 
if the Poisson bracket of their Noether charges vanishes. It is clear that in the presence of 
such a symmetry one does not only have one classical solution but a continuous set of such 
classical solutions. 

If the classical solution <pf is an absolute minimum of the potential U(fa), i.e. the trivial 
(vacuum) solution, the associated "zero-mode" 5n<j>t is the Goldstone mode, which is connected 
with spontaneous symmetry breaking and arises from the "valley" of continuous degenerated 
vacua. We want to distinguish this case from the occurrence of zero eigen-values connected 
with nontrivial solutions. Therefore we reserve the notion zero modes for the nontrivial case. 

For our static kinks the continuous symmetry is the space-translation invariance and the 
continuous family of solutions is parametrized by the position (center) of the kink Xq. As can 
be seen in the appendix this leads to a zero mode and the fact that the classical solution is 
not isolated will spoil some requirements in the quantization procedure (see below). 



2.6.2 Periodic time dependent solutions 



For time-dependent solitons things are not as simple as for the static ones and therefore one 
needs somewhat more advanced mathematical techniques. We consider non-trivial solutions 
(j)T(x,t) which are periodic in time, even in their rest-frame, with the period T: 

Ufar + U'{far) = with (j) T {x, t + T) = (j) T {x, t). (47) 

For a small perturbation r](x,t) the linearized e.o.m. again gives the stability equation 

[D + U"(4> T )] V (x,t) = 0, (48) 



but now the stability operator in (|48|) is no longer separable (invariant under time translations) , 
since the potential U"{4>t) is now time-dependent through 4>t(x, t). Because of the periodicity 
of 4>t the stability equation (|48"|) is periodic with the same period T and thus invariant under 
time translations t — > t+T. Therefore (|48|) is the field theoretical analogue of the Hill equation, 
known from point-mechanics in connection with the stability of periodic orbits |i~3|| . Because 
of the residual time translation invariance the solutions of the stability equations have special 



properties, described by Floquet's theorem [Tj|. The solutions of the stability equation are of 
the form (no summation over indices) 



e^Ux, t) with £„(x, t + T) = £ n {x, t) 



(49) 



where T is the period of the stability operator. Because of the reality of the stability equation 
the solutions come in complex conjugated pairs. A perturbative mode ([49"D after a time T has 
the form 



■q n (x,t + T) = e 



±i(Hnt+tJ,„T) 



£n(x,t) 



Vn(x,t) 
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The phases u n := jj, n T are called stability angles fl26j, characteristic exponents |15| or phase 
advance [13|], and in general consist of a discrete and a continuous set, like the eigenvalues uj\ 
in the static case. They are the generalized analogue of the frequencies u and characterize the 
periodic "orbit" 4>t(x, t). If one knows all stability angles z/„, this means to know all solutions 
of the stability equation, and one can decide for the stability of the periodic solution as follows: 

v n = real: The mode £ n stays oscillatory in time. 

v n = complex: The mode £ n grows exponentially with time and thus also the perturbation, 

the linearized equation is only valid for short times. 
v n = : The mode is a so called zero mode, the perturbation 

is constant in time. 

This classification is quite analogous to the above one for static solutions. As one can see, 
for the classical solution (pT(x,t) to be stable all stability angles v n must be real. Modes with 
vanishing stability angles v n = are called zero modes and their occurrence again leads to 
problems in the quantization procedure. As in the static case, these zero modes are connected 
with the symmetries of the system: 

Assume the existence of a continuous symmetry 7Z of the theory which leaves the e.o.m. and 
the boundary conditions (|47|) invariant. Thus if (pT(x,t) is a solution of (|47]) then (fix = Ti<t>T 
is also a solution, i.e. 

□0 T + U'{4> T ) = with 4> T (x, t + T) = 4> T {x, t). 

Since 1Z is a continuous symmetry there exists an infinitesimal transformation of 0t which is 
again a solution of fl4*7] ) and has the form 

4> T (x, t) = TZ(j) T = <f) T (x, T) + 5k4>t- 



Since 0t is a classical solution the deviation 5k4>t fulfills the stability equation (|4|) and since 
4>t and 4>t are periodic with the period T the deviation is also periodic with the period 

T. By Floquets theorem b~n4>T must be of the form (|49|) and since it is periodic with period 
T the phase \i must be zero. Thus 5k4>t has a zero stability angle v = \iT and is therefore 
a zero mode. For all independent symmetries fulfilling the above requirements there exists a 
separate zero mode. In the presence of such a symmetry there exists not an isolated but a 
continuous set of periodic solutions of periodicity T. 

An illustrative example is the Kepler problem (the above statements are of course also true 
for discrete systems). Because of the rotational symmetry one can rotate the Kepler ellipse 
in the plane and each ellipse is a periodic solution with the same period and same angular 
momentum. So one has a continuous set of ellipses and not an isolated periodic orbit of 
given period. With the time translationally invariant classical e.o.m. and BC ([47]), the time 
derivative d t (pT{x.t) is always a zero mode. 

As a summary one can say that if one can explicitly solve the stability equations one has a 
rich set of solutions of the (linearized) equations of motions. For the static case the explicit 
solutions for the 4 - and SG- model are given in the appendix. For the periodic time-dependent 



SG- breather solution ( p7|) this was done by Dashen, Hasslacher and Neveu in p6[| . Of special 
interest are solutions with discrete eigenvalues u 2 (or stability angles z/), other than the zero 
mode. The associated solution <fi c i + r\ yields a solution to the (linearized) e.o.m. which is of 
finite energy and periodic in time, also for static solitons. In quantum theory one can think of 
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this situation as a meson bound to a soliton. From the continuum eigen-functions one can only 
form a solution of finite energy by forming a wave packet (this is possible for the linearized 
e.o.m.). In quantum theory this wave packet can be seen as a meson scattering off a soliton. 

2.7 Existence of non-singular finite-energy solutions 

An obvious extension of the models investigated above is to consider scalar field theories in 
more then one spatial dimension and perhaps with more than one scalar field. Unfortunately 
this does not lead to new static non-trivial solutions which is expressed in the no-go theorem 

I: 

Derrick's theorem. Let = {fa} be a set of scalar fields in D = 1 + d dimensions whose 
dynamics is described by 

C=±d lt $-P$-U($), 

and let U be positive and zero for the ground state(s) (minima) of the theory. Then the only 
non-singular time-independent solutions of finite energy are the ground states. 

Proof. Define 

Vx 
V 2 

Both functionals V\ and V 2 are non-negative and are simultaneously zero only for the ground 
states. Let <p(x) be a static solution. Consider the one-parameter family of field configurations 
defined by 

0a(x) := 0(Ax), 

where A is a positive number. For this family the energy is given as 

E x = \^V 1 ($ x ) + \- d V 2 ($ x ). 

Since 0(x) is per assumption a solution oft the e.o.m., it follows by Hamilton's principle that 
the energy E\ must be stationary at A = 1. Thus, 

(d-2)V 1 ($(x)) + dV 2 ($(x)) = 0. 

For d > 2 this implies that both V\ and V 2 vanish, which is only possible for the trivial ground 
states </>(x) = const. For two spatial dimensions only V 2 must vanish. Since U is per definition 
a positive function this implies again that 0(x) is the trivial ground state for which U and 
thus V 2 vanishes, q.e.d. 

Well, Derrick's theorem only denies the existence of static finite energy solutions in more 
than one spatial dimension. But these kind of solutions are of special interest since they 
appear as new particle states in the quantum theory, as mentioned above. For scalar fields the 
D = 1 + 1 dimensions are very special^. The existence of nonsingular finite energy solutions 

6 1+1 dimensions are of course in general very special. 
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is also connected with the existence of topological conserved quantities (see above). In one 
spatial dimension there is "no way out" for the spatial asymptotic field values and thus leads 
to conserved topological indices. 

A way to circumvent Derrick's no-go theorem is to consider gauge fields. For completeness 
we just give some comments^ but for the rest of this work this topic is beyond our scope. As 
discussed in sect. the existence of non-trivial finite energy solutions can be proved by the 
use of topological conservation laws and is connected with spontaneous symmetry braking. 
The main results for gauge theories in three dimensions are: 



1. If the theory has no spontaneously broken gauge symmetry, the space of non-singular 
finite-energy solutions has only one component, and there are no non-trivial topological 
conservation laws. 

2. The same situation as in (1) is present if the symmetry breakdown is total, i.e. if no 
massless gauge mesons survive. 

3. If only one massless gauge boson survives (photon), the space of non-singular finite- 
energy solutions has an infinite number of components and there are non-trivial topo- 
logical conservation laws, except when the gauge group contains a U(l) factor whose 
generator enters into the expression of the electrical charge (e.g. Weinberg-Salam model). 

4. Similar results as in (3) hold if many massless gauge bosons survive symmetry break- 
down. 



Note that topological conservation laws enable us to establish the existence of non-dissipative 
solutions, not necessarily time-independent ones. Also topological conservation laws are suffi- 
cient but not necessary conditions for the existence of such solutions. It is quite possible that 
there exist non-dissipative solutions even when there are no non-trivial topological conserva- 
tion laws. Nevertheless topological conservation laws give us important informations without 
exactly solving the e.o.m., which in higher dimensions is generally not as simple as for a static 
D — 1 + 1 scalar field, where it can be obtained by quadrature. 



3 Quantization of Solitons 

In standard quantum field theory the perturbative approach starts with solutions of the free 
field equation, i.e. solutions of linear equations. Quantum effects around these free solutions 
are calculated order by order. In the case of solitons we start even classically with solutions 
of the non linear equations and then quantize around these non-perturbative solutions. The 
appropriate formalism to implement any classical fields into the quantization procedure is the 
path integral. In standard perturbation theory, by this point of view, one quantizes around 
the trivial classical solution, i.e. the solution of lowest energy called the vacuum. 

7 For more see for example || and [§J, 
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3.1 The path integral 



To calculate the quantum corrections to the energy spectrum of classical nontrivial solutions 
we will use for static solitons the trace of the time-evolution operator or the "propagator". For 
periodic time-dependent solitons one uses the trace of the Green function (WKB method). 
Therefore we shortly review some fundamental relations. 

3.1.1 Green functions, propagators and the spectral function 

The time evolution of a quantum system is determined by the time dependent Schrodinger 
equation (Schrodinger picture) 

H m)) = ikdt M*)> • (50) 

Equivalent to the knowledge of the states \if)(t)) is the knowledge of the unitary time- evolution 
operator U(t,t ) of the system, which fulfills the initial data problem 

(H - ihd t )U(t,t ) = 
U(to,t ) = 1, 

and the composition law 

U{t" ,t') =U{t" ,t)U{t,t'). (51) 

Knowing U(t, t ) means having a solution of the time-dependent Schrodinger equation ( |50| ) in 
the sense that for a given initial state \ip(t )), the state of the system at the time t is given by 

|V(t)> =W(t,t )hK*o)>. (52) 

Closely related to the time-evolution operator is the propagation kernel (short: kernel) /C = 
9(t — to)U, which fulfills the inhomogeneous equation 

(H - ihd t )JC{t, t ) = -ih8(t - t )l 
lim/C(Mo) = 1 

t—*tQ 

Thus the kernel /C is a Green-operator of the Schrodinger equation (|50|). Since /C and U differ 
only by a step function, the following relations are very similar for U. For time independent 
Hamiltoniansf] 7i for the kernel one immediately obtains the explicit solution 

/C(M ) = e{t-to)e-r^ t - to \ (53) 

For time-dependent Hamiltonians K, is a time-ordered product of infinitesimal versions of (|53|). 
For time-independent Hamiltonians, K, or U only depends on the difference T := t — t . There 
are many systems for which it is more comfortable to work with the Fourier transformed kernel 

(e > 0) 

Q{E) :=UdT ei^ T K = U°° dT e^'^ T = - 1 - , (54) 
H J hj (H-E-ie) 

8 In the Schrodinger picture TL is always time-independent for fundamental theories. This not true for the 
Dirac picture. 
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which fulfills the inhomogeneous Schrodinger equation 



(H — E)Q(E) = 1. 

In mathematics, the operator (A — -z) -1 , z e C \ spec(^4), is called the resolvent of a given 
operator A. Therefore Q is the resolvent kernel of 7i and its analytical structure gives the 
spectrum of 7i. We have added a small imaginary part to the energy ("pole-prescription") to 
ensure convergence of the integral. The singularities of Q are at e = 0, which means that the 
spectrum spec(7i) of the Hamiltonian is real. In view of (|54|) Q can equivalently be obtained 
by a Laplace transformation of the time-evolution operator U. 

Coordinate-representation. For notational simplicity and since the generalization to more 
degrees of freedom (DOF) is obvious, we first consider a one-dimensional quantum mechanical 
system, described by the Hamiltonian 

n=^-p 2 + V(q). (55) 
2m 

For practical calculations one works in the so called coordinate- or q- representation of the 
abstract Hilbert space. In this representation the spectrum of the operator q is defined as 

Q \q) — Q \o) 9^K, and (q\ 1 \q) = 5(q — q) (56) 

where we have assumed that the one-dimensional motion of the particle takes place on the 
whole real line, without additional topological constraints^. Because of the normalization 
in (|56|), both K and U fulfill first-order equations with distributional initial conditions. This 
representation of the kernel gives the (retarded) Feynman kernel also called propagator (this is 
not the "Feynman propagator" which occurs in field theory) as follows (q' = q(f), q" = q{t")) m - 

K(q",T\q') = 9(T) (q"\U(t",t') \q') = 9(T) (q"\ e^ nT \q') (57) 

= e(T) H (q",t"\q',t') H . (58) 

In the last equation we have expressed the kernel in terms of Heisenberg coordinate-states, 
which are eigen-states of the time-dependent Heisenberg operators 

QH{t) \q,t) H = q H (t) \q,t) H . 



Thus K(q' ; , T\q') is the amplitude that a particle starting at the position q' at the time t' is at 
the position q" at the time t" . With (|52|) , the time evolution of the Schrodinger wave function, 
defined by 

#?,<):= <g|V(*)> (59) 
is given by the kernel as follows (T = t" — t' > 0) 

W,t")= / dq'K(q",T\q')iP(q>,t'). 

The composition law (|51]) implies for the kernel K (t" > t\ > t') 

K{q",t"\q',t , )= f dq 1 K(q",t"\q 1 ,t 1 )K(q 1 ,t 1 \q',t'), 



9 In the regularization procedure of the field theory we will consider such topological constraints. 
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the law for the composition of amplitudes for events which occur successively in time |24|| . 

One can derive the path integral representation of the kernel (amplitude) K by multiple 
insertion of unity 



j dq(ti) \q,U) (q,U\ =: J dqi |<&) (^j 



between the Heisenberg states in ({58| ) on a time-lattice {U}. With the abbreviation U+i— U = e 
and Ne = T and after renaming positions, i.e. the initial (q = q') and the final (q^ = q") 
one, fl58|) can be written as (for details see e.g. ||22|, [^| ) 



» . N-l 

K(q",T | q>) = lim B N (T) / dq l ...dq N _ 1 ex P 4W^(^±^) 2 - V( ft )]} 

17 i=0 



/''W^"™. /"Ww (60) 

Jq',t' Jq>,t> 



The factor B(T), which determines the "measure" of the path integral, will be adjusted to 
make the integral finite and suitably normalized as iV — > oo. Thus the path integral is 
defined as a limiting process of a discrete lattice calculation. This "derivation" suffers from 



certain problems (for details see ^3|). (i) The u H(p : g)-symbol" [g^] in general depends on the 
ordering-prescription for the operators q and p (this is no problem for Hamiltonians of the form 
(|55|)) and can differ from the classical Hamiltonian and therefore the classical action, by terms 
proportional to H. This problem exists already in the operator formalism where the derivation 
starts, (ii) There also exists a time ordering ambiguity for the momentum integrations, which 



are already carried out in (60). For a different prescription of the integrations one would have to 
change the ordering prescription for the operators, i.e. choose a different if -symbol, to get the 
same amplitude. Thus the ordering ambiguities when passing from the classical Hamiltonian 
H to the operator Tt do not disappear, (iii) The existence of the (complex) measure for the 
functional integration. The problem can be skipped by evaluating the Euclidean path integral 
and analytic continuation, if it exists, to the (real-time) Minkowski space. But the existence 
of the measure is still in question. 

Nevertheless we will see that the path integral is a comfortable tool for our calculations. 
The theories which we consider here do not suffer from the operator ordering problem in the 
derivation above and thus are free from the evaluation prescription problems on the lattice. 

The q- representation of the Fourier transformed kernel, which depends on the period T instead 
of the energy E, gives the (outgoing) Green function (also called propagator)^ 

O^W,E ) = W \ H _ l E _ u) W) (61) 
Thus the singularities (for e — > 0) are poles of the Green function. The kernel K is related to 



the Green function G by the inverse Fourier transformation of (|6l|): 



K(q",T\q')= f° ^ e^ ET G{q"\q\ E). 



10 The names for K. and Q are author dependent. We use the notion kernel for /C and Green function for the 
Fourier transformed version Q. 
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Spectral representation and energy levels. In general the spectrum of the Hamiltonian 
will consist of a discrete part (bound states) and a continuous part (scattering states), i.e 

H\n) = E n \n) n = 0,1,2... N (62) 

H \p) = E(p) \ P ) pen, (63) 

where the "domains" N and B will depend on the considered system, usually is B = R . The 
states in (|62| , |63D are orthonormal and fulfill the completeness relation 

- f 

V |n) (n| + dp \p) (p\ = 1 (64) 
n Jm ' 

with the wave functions 

u n (x) := (q\n) and <p(p,x) := (q\p) . 



By inserting the completeness relation ( p4|) into the expressions for the kernel (|58|) and the 
Green function ( |61| ) one obtains a sum of the discrete and the continuous part: 

N r 

Kiq'^TW) = J2MQ"X(Q') e-^ T 6(T) + / dpip(p, q")tp*(p, q') e~^ E ^ T 0{T) (65) 

Jm 



# £(p) — E — ie J B * E(p) -E-ie 

As one can see, for the discrete spectrum we identify the poles with the bound state energies 
and the residues with the bound state wave functions. For the continuous spectrum of TC the 
Green function G has a branch cut. In the following we will use a more symbolic notation and 
write a single infinite sum for both, the discrete and the continuous spectrum. 

The generalization of the above formulas to more DOF is straightforward. The coordinate q, 
q(t) simply becomes q m , q m (t). For the Schrodinger wave functions (|59| ) this means 

ip(qi, ■ ■ ■ , q*r, t) = (qi, ■ ■ ■ , gjvM*)) • (67) 
In field theory, which can formally be obtained by the limit lim {qi(t), . . . , ?Ar(i)} — > <j>{x, t), 

N^oo 



the wave function fl67j) becomes a functional of the field ("state functional") , an exceedingly 
complicated object, especially for nontrivial field configurations such as solitons. Therefore 
we will not calculate the kernel K (U) or the Green function G, but their trace, so one never 
has to construct any state functionals. These methods are based on the work of Dashen et.al 

m. 



Spectral function and Feynman-Kac-formula. We are interested in the energy-spectrum 
of a system, especially the lowest energy in the presence of a nontrivial classical solution, to 
calculate the correction to the classical (kink) masses. In the case of static classical solutions 
this is done by investigating the trace of the time-evolution operator U or the kernel K (note 
that for T > 0, which we assume in the following, U — K). If in ( |65j) we set q" = q' := q , i.e. 
we consider closed paths, and integrate over all possible initial conditions go, one obtains for 
the traceQ , using (|6o|) 



Tr e-* HT = J dq K(q Q ,T | q ) = J^e'^ J dq | u n (q 



(68) 



ii 



We use the symbolic notation for the discrete and continuous spectrum, as mentioned above. 
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This is the spectral function of the theory, defined by H. Analytic continuation, which is related 
to the Euclidean path integral, to complex variable r — it and taking the limit r — > oo picks 
out the ground state energy of the sum in (|6S| ), thus 

lime £oT / dq K{q ,-ihT |) = k, 

T^OO J 

where k is the multiplicity of the ground state, i.e. the degree of degeneration. For a non- 
degenerated ground state (k = 1) one obtains for the limit of the logarithm 



1 f , -1 

E = lim — In / dxK(q Q , —ihr \ q ) = lim — In Tr e" 



■Ht 



This is the Feynman-Kac-formula and it allows to calculate the ground state energy without 
detailed knowledge of K. Nevertheless we will directly calculate the spectral function (|68|), in 
a perturbative calculation, and read off the energy spectrum. 



The trace in ( pq ) can be written as a path integral ( p0|) for closed paths with an additional 



integration over the initial=final position q = q^: 

/r rqo,t" 
dq K(q ,T | q ) = / dq / V q{t)ei S ^ T \ (69) 

To evaluate the trace we use approximation techniques for the path integral and obtain in this 
way approximate energy levels of the system. 



3.1.2 Stationary phase approximation (SPA) and perturbation theory 

In general it is not possible to exactly evaluate the kernel K. Thus one has to use some 
approximation techniques (perturbation theory) to calculate K. As can be seen from the path 
integral representation of K, one has to deal with functional generalization integrals of the 
form: 

F(« = /^> (70) 

This is also true in field theory. The idea is to get the dominant contribution of these integrals 
for (3 — > oo. In quantum theory the perturbation parameter is (3 = -z or by a rescaling of the 
fields the dimensionless parameter j^, respectively (see section |275|). The limit (3 — >• oo either 
corresponds to limit h — > 0, which means that the action S 3> h [semi- classical expansion) 
or A — > (weak coupling), which is the situation in standard perturbation theory, and there 
is no need for the action to be S ^> h. In this limit the integrand oscillates very fast and by 
the lemma of Riemann-Lebesgue the integral vanishes. The leading contribution comes from 
the stationary region of the phase f(x) (this corresponds to 5S = 0), i.e. from those values of 
x near xo, with f'(xo) = 0. In a first approximation we expand f(x) around xq and neglect 
terms of 0((x — xo) 3 ). This gives for ([70|) 

^ ) ^ W( "'Vm5 +0( ^ <71) 
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where we have assumed that f"(to) 7^ 0. The f"(to) = - case needs an extra examination^. 
That regions of x for which f'(x) 7^ only give contributions of O(i) can be seen as follows: 
Let f'(x) 7^ for a < x < b, then we can change to the variable z — f(x) in (f70|) . Thus 



dx e w(x) = — 
i(3 



7 e 



i/3z\f(b) 
I /(a) 



dz e tf*JL(_p-i 

/(a) dz 



Hence F a b goes to zero like ^ as (3 — > oo, where regions having /' = are of order ^ 
and therefore dominate in this limit. To illustrate the difference and the connection between 
standard perturbation theory, respectively, i.e. perturbation theory around the vacuum, and 
perturbation theory around a nontrivial stationary point (soliton) we examine exponents of 
the form 

f(x) = x 2 + v(x; A) with v(x; A) = —v(Xx) and v = 0{> x 3 ). 

A 

v(x] A) is the "interaction-term" (for example v(x; A) = —Ax 4 ). We assume that f(x) has two 
stationary points. The trivial one ("vacuum"), x = xy = with f(xv) = being the absolute 
minimum and a nontrivial one, x s ^ and x s = 0(4=) ("soliton"), for which f(x s ) is large 

relatively to the scale @ CV')- Expanding the exponent around the trivial stationary point 
x = ("the vacuum") one obtains (y = x — xy) 



F~ dy e if3y2 



n=0 



n ! 



(p(x)y 



(72) 



where P(x) is a polynomial which one gets by expanding the "interaction" v(x; A) around xy 
and factoring out the coupling A. This is a perturbative expansion in the coupling A which 
is reasonable in the weak-coupling regime (| <C 1) for which also the integrand is oscillating 
very fast, although the "action" f(x) for the vacuum is zero as mentioned above, and thus the 
stationary phase approximation is applicable. The symbol ~ indicates that this perturbative 
expansion is in general only an asymptotic series and not a convergent one ||22|| . Expanding 
now the "action" around the nontrivial stationary point x s one obtains (y = x — x s ) 



oo , n 

dy e WW £1. [pp {y) 



n=0 



where P(y) is a polynomial of 0(> y 3 ), obtained by expanding the "interaction" v around x s . 
Again this perturbative expansion is in general only an asymptotic series. This expansion is 
reasonable if j3P is small, i.e. the contributions of the deviations y of the nontrivial classical 
solution x s to the action are small relative to the scale (3 ("i")- That the integrand oscillates 
very fast is due to the nontrivial solution x s for which the action f(x s ) is large relative to the 
scale (3 ("|"), i.e that (3f(x s ) > 1 ^> (3 -> oo. 

We can collect the ingredients of these perturbation expansion as follows 



1. Standard perturbation theory (= weak coupling). The action is expanded around the 
trivial classical solution x v = 0(A°), the action f(x v ) and the energy of this solution is 

12 This case is related to the zero mode problem. 
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zero =>- x v is the vacuum. For a weak coupling A (dimensionless 1 /3 — > oo) the 
integrand is oscillating very fast =>■ both, stationary phase approximation (expansion 
around f(x v ) ) as higher order perturbative expansion (= expansion in A) is reasonable. 
The perturbative expansion is an expansion around the free "field modes" (the Gaussian 
integration in (|72|) respects only the quadratic action (A = 0) which gives the free 
Feynman propagator in field theory). 

2. Non-trivial perturbation theory (= "large action"). The action is expanded around a 
nontrivial classical solution x s = 0(-^=) =>- non perturbative, the action f(x s ) is large 

relative to the scale (3 = * =>• f(x s )/3 1 ^> (3 — > oo, therefore the stationary phase 
approximation is reasonable. The contributions of the "paths" nearby the non-trivial 
classical solution to the action are small, i.e with f(x) = f(x s ) + f"(x s )y 2 + Af(y) is 
PAf(y) = ^fi 1 < 1, therefore the ' 'semi-classical perturbative expansion" (= expansion 
in h) is reasonable, even in a strong coupling regime as long as this does not violate the 
above requirements. The asymptotic states are not free fields as in standard perturbation 
theory. 

3. SPA, semi- classical approximation. For the quadratic term f^ 2 \y) := f"(x s )y 2 being 
the dominant correction, it should be of order 4 = h, i.e. Pf^ = 0(1). In this case 
(ffT|) gives the leading correction and is called stationary phase approximation (SPA) or 
semi-classical approximation^. 



We will be mostly interested in the second case, where the nontrivial stationary points of the 
action will be the solitonic solutions of section ||. 

In the case of multiple integrals where / depends on N variables the expansion of / around 
an extremum at q = a writes as (yi = qi — eij) 

f(q) = /(«) + + 0(y 3 ) , Ay = ^7(3) (73) 

and the integral gives 

dq l ...dq N e i ^ = e™® {2m) N l 2 {\) N l 2 ^L= + 0{\) N . 

P VdetA P 

Here again we have assumed that no eigen- value of the matrix A is zero. Closely related to the 
method of stationary phase is Laplace's method for integrands of the form exp(— /3g(t)), where 
g(t) is bounded from below. Laplace's method takes the place of the SPA in the Euclidean 
path integral formulation. 

Stationary phase approximation for the path integral. The "phase" in the path integral 
is the action S[q\. Thus an approximation around the stationary phase means an approxima- 
tion around classical paths q c i(t) for which the action is stationary, 



^Ig^ = with q d (t') = q and q d (t") = q" , 



13 In the literature also the notion WKB method is used, but we want to reserve this notion for special cases 
of the SPA. 
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and the values q', q" are the initial and the final position for the kernel (|60|) . For simplicity we 
consider D = 1 quantum mechanics with the particle-action 

S[q] = £dt[±?-V(q)]. 

Expanding this action around the classical path q c i(t) {shifting method) gives 
S[q] = S[q d + rj\ 



S(q c i) + qv\o + 



dt v (-d*-v"\ qd {t))v + vv\l 



1 N 1 1 

=: S(q d ) + SS\ gcl + -S 2 S\ qcl + Efc!^ = S ^ + 6S \^ + + AS 

k=3 

The surface terms ^Vvlo an d Qv\o vanish if the classical path connects the initial and final 
position q',q" in the kernel (|60D, since in this case is 77(0) = = r](T). This is not always 
true0. The classical action S(q c i) is the action evaluated for the classical path q c i and thus an 
ordinary function of T. Therefore we will often write S c i(T) for this term. The first variation 
SS is of course zero (up to possible boundary terms) for the classical path. The operator 
O = — <9 f 2 — V"(t) is the analogue of the matrix Aij of (|73|) . The term ASi gives higher order 
corrections. Thus we approximate the path integral of the kernel ( ]60|) as follows: 

"<?",*" rO,t" 

lq',t> 

= e^'^B'(T)^^. (75) 
WdeW 

The pre-factors and the measure constant B{T) are absorbed in the new constant B'(T). 
This approximation is also called semi- classical approximation since the sum over all paths is 
approximated by the sum over the classical path q c i(t) and paths in its neighborhood. The 
quantum effects (corrections) are included in the factor 

B\T) ' 



r q",t" f 0,t" 

K(q",T \q')= / Vq{t)e^ T] w e^ i(T) / Vrj(t) e^^o v ( 74 ) 



VdetO 



When calculating the determinant detO one has to respect boundary conditions, in this case 
the homogeneous one, rj(0) = T)(T) = 0. In more general cases, i.e. if the classical path does 
not exactly connect the positions q',q", one has to choose the boundary conditions in a way, 
so that the set of fluctuations {rj} form a linear space, in which the operator O acts. This is 
the advantage of the "shifting method", that the "path integration domain" (PID) becomes a 
linear space. If no classical solution is available, the PID can be turned into a linear space 
by the following variable transformation in the path integral ||21|| 

m~#) - fT-gy-'* ]. 



14 When using functional derivatives one has to be very careful if they really exist, this is only true if such 
boundary terms do not occur. Otherwise one loses automatically these boundary terms. This is also different 
from a variation principle, which is defined by fixing the variation at the endpoints and corresponds to the 
boundary conditions for the e.o.m. In standard perturbation theory one does not have to care about surface 
terms, since one considers an unbounded time interval (— 00,00). 
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Again for this approximation we have assumed that no eigen-value of the operator O vanishes. 
Vanishing eigenvalues will lead us to the zero-mode problem and spoil the conditions for the 
validity of the SPA. 

If the action has several stationary points then each gives an additional separate contribution, 
provided the paths which make the action stationary are not too "close" to each other, since 
otherwise the condition for the SPA, that paths near the classical one, the fluctuations 77, 
give only small contributions to the action (-^ <C 1, see point 2. above), is not fulfilled. A 
characteristic length for the validity of the SPA is heuristically obtained as follows: 

If the quadratic term dominates, then we have 

^M = 0(1) as/wO 
n 

This means that 77 is on the order of yh, i.e. 

v = o(Vh) as ft ^ 



■11 •••) 



This is a relevant length for the SPA. Assuming that there exist several classical paths q a , qp, 
their distance must be larger than the characteristic length. Well, the distance is measured 
by the action. Let be q = q a + 77 = qp, then an expansion of the action gives 

S[q a + v] = S(q a ) + S 2 S a + AS a = Sfap). 

In order not to spoil the conditions of the SPA, for the difference in the actions of two classical 
paths must satisfy: 

S(qp) ~ S(q a ) = 6 2 S a + AS a » h (76) 

/5 2 S 
dt(q a - qp)j^\q a (q Q ~ qp) > ft- (77) 

Otherwise the paths are near a focal point which also leads to the zero-mode problem. Espe- 
cially in the presence of continuous symmetries these problems occur (see section |2.6|) , since 
an infinitesimal transformation of a classical path q a can give a continuous set of neighboring 
paths with the same action. 

Higher order corrections. The higher order corrections one gets from perturbation theory 
using the "rest" of the action. The interaction which is treated perturbative reads 



N , N 



AS =: S[q] - Sspa = J2y 6kS = E/ dt^% cl (t)(v) k 



(78) 



fc=3 k=3 

With this definition the path integral can be written as 

fq",T 



K(q",T\q')= Vq(t) e* SspA 

Jq',0 

= I. Vm * (-) 



I AS 



Thus for further perturbation theory AS < K must be valid. This leads to generalized Feynman 
graphs, but of course more complicated since the "couplings" V^ k '\qd(t) are time dependent. 
In field theory this was done in (2[ 
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3.1.3 One exactly solvable problem, the harmonic oscillator 

We now calculate the spectral function (|69|), i.e. the trace of the kernel, for the harmonic 
oscillator. Although this is a well known and trivial system we treat it in some more detail. 
Especially we are interested in the explicit expression of the measure, since we will end up 
with the harmonic oscillator every time. The Lagrangian and the e.o.m. (obtained by a 
variation principle, i.e. an extremum of the action, with vanishing variation of the endpoints 
Sq(0) = Sq(T) = 0) are given by: 

L = l(q 2 -u 2 q 2 ) 
q + u) 2 q = 

The classical solution with the dosed path BC q(0) = q(T) = q Q for the trace, and the 
associated action-function Sd(qo, T) for these trajectories are given by (assuming ujT ^ nil, n G 
N)0 

Qd (t) = qo (cos ut + 2 ^Jf sin ut J (79) 

S d (qo, T) = \ j Q T dt{ql - u 2 q 2 d ) = -2uq 2 ^§> (80) 

Here and and in the following we use a mode-expansion method to evaluate the path integral, 
rather than a lattice calculation. For this we expand the exponent in (|69|), i.e the action 

S[q] = ~£dt[f(t)-uY(t)] (81) 

around the classical solution (|79|), i.e. q{t) = q c i(t) +r](t). The expansion terminates after the 
second order since (|8l|) is quadratic in q{t). So the following is exact and so is the semi-classical 
calculationF]: 

S[q] = S[q d + V ] = S d (q ,T) + \ ^ dt^(t){-d 2 - u 2 } V (t) (82) 
with 77(0) = = 7](T) . . . closed path BC (83) 

Because of the closed path BC the fluctuations fulfill Dirichlet boundary conditions. The 
advantage of the shifting of the path integration over the paths q(t) to a path integration over 
fluctuations r](t) around the classical path is that the set of fluctuations i](t) form a linear 
space, as long the boundary conditions are linear relations, in contrast to the class of paths 
from g(0) = qo to q(T) = qx (the sum of two paths goes from 2q to 2qx). This is necessary 
for the functional integration of the exponent (|83|) , since this is done by diagonalization of 
the differential operator O = —d 2 — oo 2 in the path integration domain PID. Thus the path 
integration domain must be a linear space. 

The diagonalization is done by solving the homogeneous (BC) eigenvalue problem for the 
differential operator O. O is a Schrodinger-like operator and thus has an ordered spectrum 
(Ai < A2 • • • ) and a complete orthonormal set (in the sense of the space L 2 (IR) ) of eigen- 
functions. The homogeneous (BC) eigenvalue problem reads 

(-d 2 - U 2 )iJ n = e„^n , ^n(O) = VnCO = (84) 

15 For our purposes this singular situation can always be avoided, since we can choose the period T arbitrary. 
16 Surface terms again vanishes because of the closed path BC 
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The solution is easily obtained: 

T17T 

t n = k\-J 1 k n = — 7 n = 1, 2, . . . (85) 

if) n (t) = 9{T — t)J f sin k n t / dtipnipn* = 5 n>n > (86) 

The set of fluctuations PID = {q 6 C[0,T] | r/(0) = r/(T) = 0} is "larger" than the space 
{L 2 ([0,T], homogenous Be)} in which the set {ip n \ n = 1,2,...} forms a basis, since one 
has to consider all fluctuations r], even those which are not square-integrable on the interval 
[0, T]. Anyhow, we expand each deviation T](t) according to the basis {ip n {t)}. The correctness 
of the final result seems to say that the set of fluctuations not included in this expansion is 
of "measure" zero. The reason for this is that the paths far away from the stationary point 
interfere destructively (see also the comments on the characteristic length of the SPA, which 
is exact here). To have well defined expressions we make a finite expansion, which is called a 
mode regularization: 

N 

V (t) = J2 a nMt) V(0) = V(T) = 0. (87) 

71=1 

As one can see, because of the implementation of the (linear) BC on the individual modes ip n 
the full fluctuation field r\ automatically fulfills the required BC. The coefficients are completely 
free since the information on the BC is encoded in the eigen- values e n . By this expansion 
according to the fixed basis {ip n } a variation in the function r)(t) means a variation in its 
coefficients a n . Thus a path integration over r)(t) means an integration over the coefficients 
a n . For the action ( jB2|) we get 



%] = S[q d + y] = S d (q , T) + ^e n ^. 

n=l 

For the trace of the time-evolution-operator (=Kernel for T > 0) one obtains 

dq / Vq et%] = / dq e^ 1 ^ / V V (a n ) e*£ 

Jqn.O J J0,0 



The measure Vr] is given as 

JV 

Vr):=B N (T) Y[da n 



n=l 



where B(T) is an appropriate normalization constant which will be defined below. It is 
important to note that the mode-expansion-evaluation is independent of the lattice calculation 
although there exists a one-to-one correspondence at least for bosons, between them. But the 
expansion (|87[) is not an ordinary variable transformation from the lattice points q(U) = qi 
to the mode coefficients a n (24]] • Therefore the measure B(T) must be defined by proper 
normalization conditions independent of the measure given on the lattice; it cannot be obtained 
by the Jacobian of a transformation qi — > a n . The trace is now 

N N 

Tre~^ T = I dq ei scl ^B N (T)l[ f da n e^<= I dq^ 1 ^ B N {T){mh)^ f[ -j=. 

n=l n=l » n 
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The finite product of eigen-values e„ of the operator — d 2 — u 2 is the regularized functional 
determinant of this operator in the space spanned by {ip n }- For the eigenvalues of the harmonic 
oscillator (|86|) the product can written down in closed form 



N 



B N (T)(iirh)*Y[ 

n=l 

.x 



1 



B N (T) 



7T 



N\ 



n=l 



B N (T)(mh)- 
VdetO BC 

u 2 T 2 



N 



B N {T){mh)^\\ 



n=l 



n 2 TT 2 



lim 

N—fOC 



B N {T) 



(n 2 7i 2 






-u 2 


V T 2 




rpN' 




yTT J 


ivT 



smujT 



For the limit N — > oo we have used that both factors exist by themselves. For the second 
product, involving the dynamics through ui, this is a standard formula 0. For the first product 
we have assumed that the measure B N (T) is chosen in such a way that the product also exists, 
and in the final normalization we will see that this is true. Finally one obtains for the trace 



Tre"« 



HT 



lim 



B N {T) 



j},\TT N 

ivT 



7T 



d qo e-^ si ^ T ^ 



lim 



uT 
sintuT 

1 



2isin(a;T/2) 

Writing the sine as exponentials and using the formula for the geometric series one obtains 

-iuj{u+\)T 



Tre~^ HT = lim 



(88) 



u=0 



Finally we have to fix the measure by a normalization condition. For obvious reasons we 
choose 



lim 



, m\ t t n 



1 



B N {T) 



7T\f N\ 




Here we can see that the suggestive notation B(T) for the measure constant is justified, 
since it does not depend on the dynamics, i.e. on u, and is thus purely kinetic and for all 
harmonic oscillators (different a/s) the same. With this normalization we can read off the 
energy spectrum of the harmonic oscillator from its spectral function (|88|) as follows 



E v = hiij{y + 



0,1, 



This is the well known spectrum of the harmonic oscillator. The normalization-condition 
for the measure (|89|) is unique up to factors of the form e~ lcT , where c is a constant, since 
this would shift the energy spectrum by the constant c which corresponds to the freedom 
of choosing the ground state energy. Up to this freedom the normalization is unique and 



corresponds to the wave-function normalization, as can be seen from (|68|). As one can also see 
the measure constant B{T) does not exist by itself, but this is a well known situation which 
also occurs in Wiener integrals and the reason for this is that the exponential of the velocity 



term in the action is part of the functional measure [32 . 

We have treated these fundamental, perhaps trivial affairs like ID quantum mechanics and 
the harmonic oscillator in such detail because the more complicated (field) systems, which we 
will consider in the following, will always be traced back to these "simple" foundations. 
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3.1.4 Field theory 



In field theory one changes from functions q(t), depending on one parameter, to functions 
4>(x,t), depending on the parameters {x,t}. Or one can say that instead of the correlation of 
one number q(t) for each t one has an infinite set of numbers 4>(x, t) for each t. This "view" of 
quantum field theory is sometimes very helpful (see e.g. the lattice resp. the mode-expansion 
formulation of the path integral) . But of course relativistic field theory is more than the formal 
limit to an infinite number of degrees of freedom (DOF). 

Field representation. The formulas of the above sections are straightforward to generalize 
to field theory. The analogue of the coordinate representation ( |56|) , defined by eigen-states of 
the position operators in the Schrodinger picture, is the field representation, defined by the 
field operator in the Schrodinger picture (i.e. at a fixed time): 

(f)(x) \(f)(x)} = (j)(x) \(f)(x)} , 

or alternatively one cane use (time-dependent) Heisenberg operators 

Mx,t) \<j)(x),t) H = <j)(x,t) \(f)(x),t) H 

These two pictures are connected as usually by 

|0(f), t) H = ei m |0(f)). 

The kernel (propagator) is now the amplitude that the system evolves from a field configuration 
4> a (x) at t = t' to a field configuration 0b(x) at a (later) time t = t" (T = t" — t r ), and reads: 

ir(0 6 (f),T|0 a (f)) =(<f )a (x)\e-^ HT \<p b (x))= V<f)(x,t)e^^ (90) 

= H f\<j>{Z),l!) H (91) 

where the action is given by 

S[(p,T} = J dt J dx£(4>(x,t)) 

space 

The field representation is useful only for general (formal) considerations. For a lattice cal- 
culation the path integral is now defined on a space-time lattice. To show that there exists 
a unique Lorentz-invariant limit on the space-time lattice is of course a nontrivial problem, 
since there are a lot of possible kinds of lattice-structures in D > 1 dimensions. We will again 
use the mode-expansion method and assume that the functional integral exists uniquely. 

For the following paragraph we use x = x. 

Spectral function. For the trace formula ( |68|) one needs one more integration over the 
initial=final field configuration. This is again a functional integral in the case of fields. For 
the trace we again evaluate the kernel for closed paths. In field theory this means 0(x, t') = 
<p(x,t") =: (p a (x), and integrate in addition over the initial=final field configuration (f) a (x). As 
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before we insert a complete set of energy eigen-states (symbolic notation) in (|90|). So we get: 

I V<j> a (x)Y^ (<f>a(x)\n) (n\<f> a (x)) e~^ E " T (92) 

= J2( n \ (/ \n)e~i E " T (93) 

= Y t e-* E » T = Tre-l n . (94) 

n 

Here we have used the completeness relation 

V(f> a (x)\<f> a (x))(<f> a {x)\ = l, (95) 



which can be obtained by the limit lim of the unit in the Hilbert space of N degree of freedom 

dqi...dq n \qi...q N >< qi...q N \ for all t 
by identifying the q^s with the values fa on the space lattice and the continuum state \4>(x)) = 



lim \(fii...(f) N ) . One can also read this in a different way. fl93|) also shows the normalization 

of the state functional, i.e. the wave function in the field representation, the exceedingly 
complicated object which one wants to get rid of: 



V<f>(x)^f*[<f>(x)]^[<f>(x)] = 1. 
Thus one obtains for the spectral function: 



T 

x,t)e* s ^. (96) 



K(T) = Tr e-^ n = [ P0 a (x)K(0 a , T\ a , 0) = f V<j> a (x) [ " V<j>{ 

J J J<j>afl 

For the calculation of the functional integral we again use the stationary phase approximation. 
By the integration over the field space one also has to respect spatial boundary conditions for 
the fields. 

Spatial boundary conditions. We shortly examine the influence of the spatial boundary 
conditions on (|96|) for theories in D=l+1 of the form 

£ = ^0) 2 -W) (97) 

where for all minima U m in — is valid, 
(i) unbroken symmetry: 

In the case of a unique minimum of U the minimum should lie at = 0, which can always 
be reached by shifting the field. From (||) one can see, that for = also the energy is zero. 
Thus we expect the quantum vacuum state \vac)^ at <f) = 0, i.e. \vac)^ = \<p(x, t) = 0). The 
boundary conditions for finite energy solutions are (0) 

(j)(x -> ±oo, t) = (98) 
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and so the Fock space should only exists over functions satisfying fl98| ), i.e. we have (Fock) 
states which are located around functions <p(x — > ±00) = 0. Correspondingly in the 

functional integral one has to integrate only over fluctuations around a stationary point sat- 
isfying fl98T ) for all t. Thus in the case of an unbroken symmetry one can only perform a 
perturbation theory around the vacuum, i.e. standard perturbation theory in our framework. 

(ii) spontaneously broken symmetry: 

In this case one has several minima U((f>i) = , i = 1...M which gives rise to nontrivial 
topological sectors (see section |2T]) . The boundary conditions for finite energy solutions 
classify the fields topological and are denoted by (0) 

4>(x — > ±00, t) = (f>i±. 

Since the different topological sectors are not connected the trace and thus the functional 
integral fl96|) has to be evaluated for each sector separately. The topological charge, which is 
not changed by "quantum fluctuations" (see (|29l)), acts as a super-selection quantum number. 
One has to integrate over fields appropriate to the sector, i.e. over fluctuations around a 
stationary point of definite topological charge. The completeness relation ( |95| ) holds in the 
subspace according to the topological sector. This will be justified below, where we construct 
the Hilbert space of the nontrivial sector. 

3.1.5 Quantum energy levels for static solitons 

We now consider the SPA for theories of the form fl97l) which permits one static soliton solution 
4>d — 4>ci( x ) in a topological sector S, i.e. we neglect the zero-mode problem^. We also restrict 
our considerations to a finite space region, which has to do with the regularization procedure 
(see below). Thus the action is given by 

S[<j>] = / dtdxC(cj)), 
Jb 

where B = LxT is the finite space-time region. The spectral function is given by the 

Trace-cum-path integral 

K{T) = Tr s e~i HT = I V[<j> a {x)\ I ^ V[<f>(x,t)\<* m - 

J J<t>a,0 

top. sector 

The first integration sums up all contributions of closed paths with start- and end- point 
4> a { x ) m the considered topological sector. The second one adds all these contributions for 
all starting (end) points in the topological sector. The SPA will pick out the contribution 
of fields in the neighborhood of C ;. This set is characterized by the closed path condition 
(f)(x, 0) = (f)(x,T) = 4> a (x) and that only second order deviations O(0(x, t) — (p c i) 2 count. 

Stationary phase approximation. We approximate the action (the phase), which is sta- 
tionary for c /(x), around this classical solution to evaluate the functional integral and the 
trace. Therefore we consider fields 

4>(x,t) = 4>d{x) + r)(x,t). 

17 If there are more, well separated classical solutions in a topological sector, the spectral function is the sum 
of each contribution. This case has to be distinguished from the existence of zero modes. 
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Since 4> d (x) is static, the closed-path condition implies for the fluctuations r\: 

<f)(x, 0) = (f)(x, T) = <f) a (x) r}(x, 0) = r](x, T) = r} a (x). 

The spatial BC will be determined below. They are essential ingredients of the regularization 
process. Expanding the action around the stationary point (j) c i, according the considered 
topological sector, one obtains 

S[(f>, T] = S[<f> d +v] = S{(j) ch T)- \ f dtdx r](n + U"(<f> d )) rj (99) 

^ J LxT 



+ 



{-drj + d(f) d )r) 



\8B 



+ o(\\v\\% (100) 



The boundary terms are not vanishing, in contrast to the above sections, since the classical 
solution (pd are not exactly identical with the initial=final field configuration (p a (x). Since <\> c \ 
is static, the classical part of the action gives (using (||)) 

S{<j> d , T ) = -J Q dtj dx[^<f> d ' 2 + U{<j> d )\ = -J dtE[<j> d ] = -E[(j) d ]T = —M d T. (101) 

With the translation 0(x, t) — > r](x, t) = (f)(x, t) — 4> d of the integration variable we get for the 
trace 

K{T) SPA = e* s{M V[ Va {x)] / V[r](x,t)}e-iy^Tdtdx V (a+u")v+-_ ( 102 ) 

The dots stands for the boundary term. The operator in the exponent, 

0(x,t) = dl -dl + U% cl {x) =: d 2 t + SO{x) 

is separable since <fi c i depends only on x. Therefore we expand, analogous to the harmonic oscil- 
lator, the paths according to eigen- functions of the spatial part, but now with time dependent 
coefficients. So we have 

{-d 2 x + U"{<j> d )) £ n (x) = ulUx) With J dxCM = 6m,n (103) 

r}(x,t) = y^c n (t)£ n (x) and r} a (x) = ^c a ^ n {x). (104) 

n n 

The operator in (|103|) is a Schrodinger operator, and thus the eigen-functions {£ n } form a 
complete set. The spatial BC i.e., £(— L/2), £(L/2), will be specified below. Also we leave 
the explicit form of the sums in (|104|) open, since this will also be part of the regularization 
procedure. So we get for the spatial part of the exponent in ( P-02j ) 

J dxr)(x,t)Or)(x,t) = £)c?(*) (J dx$(x)(% + SO)U%)\ c k (t) (105) 

= ^(*)(^ + w ?)^)- ( 106 ) 
i 

0For the boundary term we assume that the spatial boundary conditions do not introduce 
any contributions, which will be justified in concrete calculations e.g. by the use of topological 

18 In the second line we gave set c\ = c;. From the eigen-functions in the appendix (7.1) one can see that 
the reality condition for the field r\ for the continuum modes is c[ = c_;. But by a u nitary transformation 
Q — > UikCk, which leaves the path integral invariant, one gets real oscillators q. In (106) and in the following, 
it is assumed that this transformation is already carried out, after the spatial integrations. 
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boundary conditions (see below). Nevertheless the time-like boundaries induce contributions, 
because the fluctuations are only closed paths and not periodic ones. With this assumption 
and the time-independence of (f> c i the boundary contribution in (|100|) is 



1 1 f^' z 

dtdx df,[d^((f) c i + —77)77] = - / dxr) a (x)[r)(t , x) - 77 (t , x)] 

Z Z J-L/2 



Together with (|106| ) this gives for the exponent in the path integral (|102|) 

-\t I dt cl (t)(d?+co 2 ) Cl (t)+^2 1 dtd t ( Cl ci) = It I * (^)-^w) 

j Jt l Jt l Jt 



This is the sum of harmonic oscillators, each with the action as given by (|81|) and with the 
closed path BC q(0) = q(T) = c a> i which is the analogue of g(0) = q(T) = q and therefore 
no further boundary contributions occur as shown for the harmonic oscillator above. The 
measure is therefore given by a product of harmonic oscillator measures 



Vr ] (x,t) = l[Dc l (t) = ll 



N 



B N (T)Y[dci tn 



n=l 



where for each oscillator the measure constant is the same, since it is independent of the 
oscillator frequency uji as discussed above. The occurrence of a zero mode, i.e. oji = must be 
treated separately, since in this case the action of this mode is no longer that of an harmonic 
oscillator but that of a free propagating particle in one dimension. 



So we get for the trace (|102|) 



K SPA {T) = Tr e-i HT \sPA 



n / dc ^ 



Vci(t) eTiiT 



(107) 



where the sum in the exponent is now written as product of exponentials. 

Comments: (I) In the SPA the trace is a product of the classical part with an "infinite" set 
of harmonic oscillators. (II) the system was made discrete by the introduction of spatial BC 
which has to be specified and will be an essential part of the regularization procedure. They 
should be chosen in such a way that no spatial boundary contributions ( p.00|) to the action 
occur. This means that they should be "topological" (see below). (Ill) The derivation above 
is also valid for constant, i.e. trivial, classical solutions (fid = (fiy = const, and therefore the 
trace for the vacuum sector gives analogous results. The only difference is the eigen-value 



problem (|103| ), which in the case of a trivial solution is of course much simpler and gives 
different eigen- values uj 2 . (IV) We have excluded the possibility of a zero mode uj 1 = 0, which 
leads to subtle problems and must be investigated separately. 

Energy levels. Assuming that no eigenvalue vanishes (no zero mode) we get with the result 
of the harmonic oscillator (|88|) for the spectral function ( [1071) in the topological sector S 



l ,HT\ 



\SPA 



IIE' 



-iuif{ui+\)T 
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where / is the mode index and vi is the excitation index of the /'th mode. Thus a general state 
has the energy-spectrum^] 

EslM] = EM + ftj>f fa + \) + 0(h 2 ). (108) 

i 

This energy-spectrum formula is valid for trivial solutions (f) c i, like the vacuum, and nontrivial 
solutions, i.e. solitons. In both cases the lowest energy-level is given by the state where no 
mode / is excited, i.e. v\ — for all I. These are the ground states in the considered sectors S 
and given as 

E s = E[<j&+^£d* + 0(h 2 ). (109) 
i 

We have now all ingredients to calculate the energy correction for solitons except the renormal- 
ization contributions to the quantum-action. As we will see, the semi-classical approximation 
(=SPA) in the nontrivial sector (the soliton) is already a one-loop result, i.e. order h, and thus 
one has to renormalize the theory to get control of UV-divergences. This will be considered 
in one of the following sections. 

In a last comment we want to outline the quantum nature of the SP A-correction, beside the 
occurrence of h. Since the quadratic part of the action (|103| ) is exactly the stability equation 



fl42|) , the functions (pd + c n (t)£ n (x) are wei l behaved nearby classical solutions (see |1T|) if the 
eigen-value uj\ is positive. The quantum nature of these fluctuations, beside that their action 
occurs as a phase, is that the oscillations are not classical oscillations ~ e luJt but treated as 
quantum-oscillators through the path integral in ( |107| ). This can also be seen from the full 
stability equation (|40D , where the classical nearby-solutions are those with "eigen-value" zero 
of this equation and for the diagonalization of the quadratic action in the path integral we are 
treating the eigen-modes of this equation with non-zero eigen- values. 



3.1.6 The zero mode 

We have shown above (section [2^ ) that the occurrence of zero modes is connected with 
symmetries of the system and the classical solution. For our special models, the SG and 4> 4 
theory, such a zero mode occurs in the kink sector. It is the lowest eigen-value of the stability 
equation (|7.1|) . These zero modes are connected with the translational symmetry of the kink 
solutions ( |35l , [20| ) . For any position of the kinks Xq the equations of motions are fulfilled. In 
both cases the zero modes are proportional to the spatial derivative of the kinks, i.e. 

7]q(x) ~ d x (f) K (x), 

and thus are a result of a small translation of the kink. In the path integral quantization this 
results in zero frequency u = in (|107|) and thus the associated degree of freedom c is not a 



harmonic oscillator but rather like a free particle of unit mass. The reason for this is that the 
kink solution is only a local minimum of the action and the potential energy(density) U ((/)), 



19 That also with the nontrivial sector, i.e. for 4> c i a soliton solution, are quantum states associated will be 
considered below. 

20 The notion "energy-spectrum" should make clear that this is not the energy of the state |{^}), since there 
are some more ingredients for the energy like renormalization and zero-point energy (see below). 
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respectively. Thus a fluctuation in the symmetry direction does not change the energy and 
feels no restoring "force". Fluctuations transverse to this symmetry direction (in field space) 
feel the restoring force of the increasing potential, which is in first order that of a harmonic 
oscillator. 

The trace integration of the zero mode gives a divergent result. The only closed path with 
initial = finial position c ao for a free propagating particle is the constant solution c c \ = c ao . 
The classical action J dt-%- is zero for this solution. Thus the trace integration of the zero 



mode c in (|107|) gives 

B N (T) J dc ao -> oo. 

The factor B N {T) is the measure constant of a free unit mass particle. The breakdown of the 
SPA is no surprise, since for its validity we had required that two classical solutions are not 
too close to each other (|T7|). But because of the translational invariance of the kink solutions 
we have a continuous family of solutions, parametrized by the kink position xq. The free 
propagation of the zero mode degree of freedom corresponds to a collective motion of the 
kink and its internal quantum fluctuations. It is customary to treat the zero mode by the 
use of appropriate coordinates to describe the symmetry, called collective coordinates. This is 
analogous to for example atomic physics where the collective center of mass motion is separated 
from the internal motions described by relative coordinates. The idea is to find coordinates 
which describe the motion in the symmetry direction, i.e in the "valley" (surface) in field space 
which forms the relative minimum of the action, these are the collective coordinates. The 
method of SPA is only applicable to the residual coordinates, which describes the internal 
motion. The integration for the collective coordinates has to be carried out exactly. We 
demonstrate this for a simple integral: 

r d n x e g2 - x ^ )2 . (110) 
The exponent is stationary for the "classical" solution 

x d (l - 2XxV) = =>• \x\ ' 



2A 

This is an n — 1 parameter family of solutions, each vector with the length This cor- 

responds to the 0{n) symmetry of the exponent in ( |110| ). If we single out one stationary 
(saddle) point and evaluate its contribution in a Gaussian approximation we would get n — 1 
zero eigenvalues and thus a catastrophic, divergent result. The solution to this problem is to 
use angular variables and integrate them exactly. Only for the radial variable one expands 
around the stationary point. Thus one obtains 

dtt^r^dr e r2 - Xr4 = V S n-x [ drr n ~ x e r2 ~ Xr \ 



The radial integral can now be evaluated using Gaussian approximation. The angular integral 
is evaluated exactly and gives the volume of a n — 1 dimensional sphere. Here the angular 
variables are the analogue of collective coordinates or cyclic coordinates, as they are called in 
classical mechanics. 
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The proper collective coordinate for the kinks is of course the position of the kink, which 
we will call X(t). This means that we change from the coordinates {ci(t);l = 0,...} to 
coordinates {X(t), Ci(t); I = 1 . . . } to get rid of the problematic zero mode Co(t). This is done 
by expanding the field according to 

<t>{x, t) = <j> K (x - X{t)) + J2ci{t)ti{x - X(t)), (111) 

1=1 

where <\>k and £ are the same functions, i.e. the kink functions and the eigen-functions of 
the stability equation except for the zero mode, as before. That the new coordinate X(t) 
consistently replaces Co and is thus independent of the other coordinates can be seen as follows: 
A small variation of X(t) adds to 4>(x,t) a term proportional to d x (f>K(x — X(t)), i.e. 

5 x< P{x,t) ^ d x (t>K{x - X{t))5X. 

But the derivative of the kink is proportional to the zero mode and thus 

dxd x (j) K {x - X(t))&(x - X(t)) = V I. 



Therefore a variation in the collective coordinate X(t) is orthogonal to the other coordinate 
directions and the set {X(t), C[(t);l = 1 . . . } consists only of independent variables. Because 
of the translational invariance the new expansion ( |1 1 1| ) only changes the kinetic part of the 
action and the Lagrange function, respectively, i.e. the spatial integral of the Lagrangian: 



L = J dxC = J dx [{dtcj)) 2 - (d x <p) 2 - U 
= LKin(X, q) + Lp ot {ci). 



Thus the Lagrange function is independent of the coordinate X(t), it depends only on the 
velocity X(t). This is true for the full Lagrange function and action. Especially the quadratic 
part of Lp ot gives the same quadratic action as in (|107|) except for the zero mode Co(t). Since 



L depends only on the derivative X(t), the collective coordinate is a cyclic coordinate [|T^] and 
thus the proper coordinate to describe symmetries of the system. From the Euler-Lagrange 
e.o.m. it follows that the canonical momenta of cyclic coordinates are conserved, i.e. on a 
classical trajectory (on shell): 

dL , d ( dL\ d 

— = Vi cyclic - Pw , = 0. 

Thus the canonical conjugated momenta of cyclic coordinates are the Noether charges of the 
associated symmetries. And as expected the canonical momentum of X(t) is equal to the 
conserved total field momentum P fl3|, i.e. the Noether charge of the spatial translation 
symmetry: 

Tlx = = P with ^-P = 0. 

dX dt 

Since we are considering relativistic theories the classical energy of the kinks has the form 



E{P) = JMl + P*. (112) 
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Therefore the integration of the collective coordinate in the path integral gives only kinetic 
contributions which one can neglect in the considered order in (|107|) ||]. This means that the 
kink is effectively at rest. This is a reasonable approximation especially for the calculation 
of the quantum mass of the kink in this order . In higher orders the different modes in the 
spectral function are no longer independent. For ( |107| ) this means that the fluctuations for 
/ > interact with the zero mode fluctuation which is then no longer a free (zero) mode. 
But for the one-loop (= SPA) calculation of the kink masses we simply omit the integration 
over this mode. Thus, because of the translational zero mode, one has one mode less in the 
kink spectrum. This has a completely different origin than the discrete excited mode £1 of the 
4 -kink (see appendix (|7.1|) ) but an similar consequences in the bosonic case, as we will see. 

The interplay between cyclic (collective) coordinates and the conserved associated momenta is 
to be expected to be more subtle in the case of constrained systems like fermions. And indeed 
we will see that this naive counting of the zero mode in the energy spectrum leads to wrong 
results in the case of fermions. 



3.2 Standard perturbation theory and renormalization 

We shortly survey the main points in the renormalization procedure. For simplicity we do 
this mostly for the 4 - model. One crucial point for a consistent renormalization is that one 
renormalizes the theory only once (at a given perturbation order) and uses this then fixed 
renormalized theory to calculate the desired quantities also in different topological sectors. It 
is customary to renormalize the theory by setting up renormalization conditions in standard 
perturbation theory, i.e. relations between scattering process contributions. This associates 
the parameters of the theory with certain measurable physical processes. So one can (must!) 
determine them by experiments. In renormalizable theories a finite number of such conditions 
also eliminates the divergences in the standard perturbation theory. 

The standard perturbation theory is tailored for calculation of amplitudes in scattering pro- 
cesses, which are related to (vacuum) n-point correlation (Green) functions or their Fourier 
transformation, respectively, which are given as (x = x^) 

G n (xi, ...,x n ) = T(j) H (xi) . . . (j) H {x n ) = (0\S\0) ' ^ ' 

The index H/D stands for the Heisenberg/ Dirac -picture and T is the time-ordering symbol. 
S is the S'-matrix (operator) which is related to the time-evolution operator of the Schrodinger 
equation (|50|), and reads in the Dirac-picture as follows 

S := lim U(-oo, oo) = lim T e -*f-°° dtn h^ . 

The limit e — > 0+ indicates the adiabatic "switch on" of the interaction, i.e. H £ d '■= Ti.D^~ 6 ^, 
and corresponds to the boundary conditions (preparation) for scattering processes (asymptotic 
free particles (fields)). For details see e.g. pSfl . This is the reason why in standard perturbation 
theory soliton contributions are not seen. Of particular interest is the two-point function 

G{x,y) = {n\T<f> H {x)<f> H (y) 
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Figure 8: The spectrum of TL plotted against the spectrum of "P. At the origin sits the vacuum state 
On the ordinate lie the "masses" of the discrete one-particle- and bound- states, as well as a 
continuum of multi-particle- states (above the highest hyperboloid). For non-zero momentum they 

form hyperboloids according to the relativistic energy relation E(p) = \J p 1 + m\ , which asymptoti- 
cally approach to light cones. There may be also more bound states bellow the threshold of two free 
particles-creation. 



3.2.1 Analytical structure of G(x,y) and field strength renormalization 

First we consider the spectrum of the Hamilton operator TL and the momentum operator Vi 
(we consider a D = 1 + 3 space-time, so that Vi is a three- vector and V th = V = (TL,Vi) T ). 
Since they commute, i.e. [TL,Vi] = 0, they have common eigen-states. The vacuum state \Q) 
is the eigen-state to the eigenvalue zero, i.e. TL \ — = Vi |fi). Let |A ) be eigen-states of 
zero momentum, i.e. 

ft|Ao)=m A |A > (114) 
Vi |A > = 0, (115) 

then the boosted states Ub oost (p) |A ) = \\p) are also eigen-states, but with momentum p and, 
because of relativistic invariance, energy E(p) = ^Jp 1 + m\ . Thus the eigen values m\ are 
the energies in the rest-frame. In general the spectrum consists of the vacuum, the one- 
particle state (nix = Tn, particle mass), possible bound states (m\ =mg) and a continuum of 
multi-particle states (see fig.||). Thus the completeness relation, expressed in relativistically 
normalized states, in the Hilbert space reads as |3( 



A 



For the following we assume that x° > y°. Inserting the unit ( pL 1 6| ) we thus get for the 
two-point function 

G{x, y ) = (n\Mx)My)\n) (117) 

dp 3 1 



E 

A 



(0| 4> H (x) \\p) (Xp\ 4>h(v) 1^) , (118) 



where we have dropped the uninteresting constant (0| 4>h(x) |^) (^| 4>H(y) |^), which is usually 
zero |30|| . The matrix elements in (|118| ) can be written as follows 

(fi| <j> H {x) |V> = m ^M^^ \\p) 

= (n\Mo)\^p)e- t H P °=E(p ) 
= (n\Mo)\^)e- i n P »=E(p). 
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Figure 9: The spectral density for a typical interacting theory. The one-particle states contribute 
a delta function at m 2 . Multi-particle states have a continuous spectrum, starting at (2m) 2 . There 
may also discrete contributions of bound states. 

In the first line we have written the field operator at space-time position x as a translation 
of the operator at the space-time origin. In the second line we have used the translational 
invariance, (Q\ e lV ' x = (fl\ e°, of the vacuum and the Lorentz invariance of 0#(O) and (fl\, i.e. 

UbooitffiMWtoLip) = <M0)B so that (Q\ <t> H {0)U bO ost{p) |A ) = U boost {$<t> H {0) |A > = 
(Q| 0^(0) |A >. 

Kallen-Lehmann spectral representation. By introducing an p° -integration the momen- 
tum p becomes off-shell and the two point function can be written as (x° > y°) 

(0| m*)4>h(v) W = £/ ^Y v 2_ ml + l( \ ^ <M°) |A ) I 2 - (H9) 

A 

Here the Feynman propagator Dp(x — y, m 2 ) appears but with m\ instead of only the particle 
mass m. An analogous expression holds for x° < y° so that the full two point function is given 

by 

G(x, y) = G(x, y) = (0| T4> H {x)4> H (y) \Q) = J —p{M 2 )D F {x - y, M 2 ). 

This is the Kallen-Lehman spectral representation, where p(M 2 ) is a positive spectral density 
function, 

p(M 2 ) = ^2vr5(M 2 - m\)\ (Q\ H (O) |A ) | 2 , 

A 

For a typical interacting theory it is given by fig.^. Stable one-particle states |A™) contribute 
an isolated delta function to the spectral density: 

p(M 2 ) = 2n5(M 2 - m 2 )Z + terms with M 2 > m 2 B , 

where Z — \ (Q| 0^(0) |A™) | 2 is the field-strength renormalization factor. The quantity m is 
the exact mass of a single particle, since it is the exact energy eigen- value of the full interacting 
theory, as can seen by (|114|) . This quantity in general differs from the mass-parameter used 
in the Lagrangian (see below) and we refer to it the physical mass of the 0-boson. 

21 For fields with spin one has to respect the nontrivial internal transformations of the field, which are of the 
form UtfiaU^ 1 = S a /3<f>f}. For spinor fields this gives the matrix structure + m in the propagator. 
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A Fourier transformation of the two-point function (|119|) gives 



00 dM 2 , . _ 2 % 



r r°° dM 2 

J dx 4 e ip - x (fi| TfaWMO) l^> = J ~^P( M 



p 2 - M 2 + ie 

~ p 2 -m 2 + ie + J^ m 2 B 2vr P ^ > 2 -M 2 + ie 

The analytical structure of the Fourier transformed two-point function is as follows: The first 
term gives a simple pole at p 2 = m 2 with the residue Z, while the second term contributes 
a branch cut beginning at p 2 = (2m) 2 and additional poles for possible bound states below 
the cut. Thus contributions of from one-particle and multi-particle intermediate states can 
be distinguished by the strength of their analytic singularities. This analysis relies only on 
general principles of relativity and quantum mechanics , it does not depend on the nature of 
interaction or on perturbation theory, except that we have used scalar fields due to notational 
simplicity. The only input of standard perturbation theory is that we have considered the 
vacuum-correlation function and thus soliton contributions do not occur, since the Hilbert 
space built around the soliton is not connected with the "vacuum" Hilbert space. There 
exists no operator which connects states between these two sectors (see below). This analysis 
generalizes to higher order n-point functions. The analytical structure shows that the (Fourier 
transformed) n-point functions are the multi-particle (field theoretical) analogue of the kernel 
resp. the Green function of section |3.1.1| . 

In the case of free fields, or zeroth order perturbation theory, the Fourier transformed two 
point function writes as 



/ 



dx^r* (0| T0 D (x)0 D (O) |0> = - 2 — . 

P 2 - m free + " 



For x° > this can be interpreted as the amplitude that a particle created at the space- 
time position y = propagates to x. It is similar to the full two-point function except two 
differences: There are no multi-particle contributions since free fields can only create single 
particle states. The field strength renormalization constant Z = \ (Q\<pjj(0) \\™) | 2 , i.e. the 
probability for 0#(O) to create an exact one particle state, is in the free case equal to one, i.e. 
(p\ </>d(0) |0) = 1. Note that by a renormalization of the field strength this probability can be 
normalized to one also for the interacting theory: 



in - := ^=<P H 



I <n| 0no) |as»> I 2 = \\ (o| mo) K) I 2 = § = 1. 

With this renormalization also the residue of the single-particle pole in the two-point function 
(propagator) is normalized to one: 

/ dx*jr* (n\ T<j> r ™(x)(f>T(0) \n) = \ [ dx 4 e ip - x (Q\ T<Ms)<M0) = ~ + • • • 

J Z J p z — mr + ie 

3.2.2 The systematics of renormalization 

Primarily the renormalization has nothing to do with the occurrence of divergences (in per- 
turbation theory), but to express the theory (action) in terms of measurable quantities, i..e. 
in terms of parameters which are related to certain (reference) experiments. In he following 
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we assume that the theory is regularized in some way, so that all considered quantities are 
well defined. This we will indicate with an index A. We indicate the unrenormalized (not ad- 
justed to the reference-experiments) parameters with an index to distinguish them from the 
renormalized ones (parameters which are related to the reference experiments). The selected 
reference experiments and their relations to the renormalized parameters define the values of 
these parameters in the measurement and are called renormalization conditions. For simplicity 
we consider the 4 theory without spontaneous symmetry breaking, The fields can be viewed 
as ordinary functions, since we will use the path integral for further calculations: 



(120) 



This Lagrangian and the associated action has all symmetries which are not destroyed by the 
regularization. We have seen that we can renormalize the field-strength. Therefore we set 



>o ^ 



(121) 



where the constant Z\ need not be equal to Z above, although this would be a natural choice. 
Its actual value will be defined by a renormalization condition. One can also use the freedom 
in fixing Z\ to make calculations simpler instead of getting a simple relation to the reference 
experiment. With this first renormalization the Lagrangian writes as 

4 = ^)-^-^. 



We have just inserted (|121|) in the Lagrangian, thus it is still the same theory. Next we split 



the mass and coupling parameter into a renormalized (measured) part and a part which will 
be fixed by the renormalization conditions in concrete calculations, i.e. we divide them into a 
part which will be measured and a part which will be calculated. We also write Z\ as such a 
sum, i.e we define 



m 2 Z 1 
\Z\ 

Inserting this into the Lagrangian we get 



m 2 + 5m 2 
X + 5X 
l + 8Z x . 



(122) 
(123) 
(124) 



(125) 



Of course this Lagrangian has the same symmetries as ( [120]) , i.e. those symmetries which de- 
spite regularization are still present. The Lagrangian consists of the "classical" part, expressed 
in terms of the renormalized parameters, and counter terms (the 5-terms) which are treated 
as interaction terms in the perturbation theory. Thus classical properties, as for example the 
classical kink masses (|22|) in the model considered above, are not affected. Also the funda- 
mental ingredient of the standard perturbation theory, the Feynman propagator (see below), 
is not affected, since it is derived from the quadratic part of the classical action. In this view 
the counter terms are the quantum contributions to the full quantum action (|125| ) relative to 
the classical action which is expressed in renormalized parameters. The additional interaction 
terms, the counter terms, lead to additional Feynman graphs in the perturbation theory: 
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= i(p 2 5Z 1 - 5m 2 ) 
= -iSX 

Until now the splittings of the parameters ( |122| ) - (|124f) are purely formal. To give them a 
physical meaning we have to set up renormalization conditions which relate the renormalized 
parameters to certain scattering processes so that they can be determined by experiments. As 
an example we choose: 



-O- 



+ (terms regular at p 2 = m 2 ) 5m 2 , 5Z 




\s=4m 2 , t=u=0 



-iX =>- 5\ 



These relations for the "scattering amplitudes" displayed by the graphs above determine the 
explicit form of the counter terms. Especially the first renormalization condition is very 
"natural", since it renormalizes the residue (the probability of a single particle creation of the 
vacuum) to one and sets the renormalized mass parameter m to the pole of the single particle 
contribution of the two point function. As we have seen above, this is the exact physical 
particle mass, i.e. the eigen- value to the one-particle eigen-state of the full Hamiltonian (|114|). 



But there is no need to choose these "physical" renormalization conditions. One can also choose 
less simple relations between the two-point function and the renormalized parameters, which 
are more comfortable for other calculations for example. The second renormalization condition 
has no "natural" or obviously best definition. The renormalization condition is set up for the 
scattering amplitude at zero momentum, which means in Mandelstam variables s = Am 2 and 
t = u = 0. The values of the variables p,s,t,u at which the renormalization conditions are 
defined are called renormalization point. Different renormalization points leads to different 
renormalization schemes. They are of course all equivalent, since a different renormalization 
point only changes the division between the renormalized parameters and the counter term 



constants in (122) - (124). One can always put a finite piece of e.g. m into 5m . Thus a 



different renormalization point (different renormalization condition) only results in 

m 2 + 5m 2 — > m 2 + 5m 2 = (m 2 + A) + (5m 2 - A). 
This of course gives the same results but expressed in terms of m 2 . 

Renormalizable theories. We were cheating as we said that the renormalization procedure 
has nothing to do with the occurrence of divergences, since this procedure was developed to 



control divergences in the perturbation theory. The classical part of the Lagrangian ( |125|) 
produces divergent contributions (if the regularization is removed) in the perturbation theory. 
Thus also the counter terms must be divergent (without regularization) to compensate the 
divergent contributions by the renormalization conditions. A theory is called renormalizable 
if all divergences in the perturbation theory can be compensated by a finite number of renor- 
malization conditions, thus the results are finite even if the regularization is removed at the 
end of the calculation. In a perturbative expansion the counter terms must be determined 
order by order, this means in each order one has to solve the renormalization conditions, given 
above. 
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Symmetries. If all divergences can be compensated by counter terms obtained by renor- 
malization of the parameters, as above, the quantum corrections do not violate any of the 
symmetries which are left in the regularized theory C\. The question is weather the sym- 
metries which are broken by the regularization are reestablished when the regularization is 
removed after calculations. If a symmetry can not be reestablished one calls this an anomaly. 
In a perturbative expansion in K (loop-expansion) each symmetry of the theory, which does 
not involve h, is fulfilled in each order, i.e. if one respects all contribution up to the considered 
order. If not all divergences are canceled by counter terms of the form of the Lagrangian 
one needs additional counter terms, which must be added "by hand", e.g. a term C(p 3 . The 
additional parameters must be determined by experiment and the symmetry of the theory is 
probably broken. 



3.2.3 Renormalization of (h 4 - and SG- model 



Generating functional. The n-point vacuum correlation functions (|113|) can be written as 
a path integral as follows |3( 



x [V(/)<t>(x 1 )...<t>(x n ) eif-T*** , , 

G n (x 1 ,...,x n )= km . (126) 

The denominator is the analogue of (0|«S|0) in ( |113|) , which connects the Dirac- vacuum |0) 
with the Heisenberg- vacuum The path integrals have to be evaluated in the vacuum 

sector, this means that 

//•<j>n,T r<t>v,T 
V<p = V(p= lim / V<j>. (127) 

J<ki-T r-»oo(l-fc) ./^ _ T 

The n-point function ( |126|) can be written as functional derivatives of a generating functional 

G n (xi, . . .x n ) = — t-t I - J — — r — — - \ j=0 , (128) 

Z[0] \ij 5j(x n )...5j(x 1 ) 



where the generating functional is given by 

and must be calculated perturbatively for nontrivial theories. The ie-term is the analog of the 
small imaginary part of the time in ( |126|) to ensure vacuum boundary conditions. It acts as a 
damping factor and gives the pole-description for the Feynman propagator. In the following 
we will suppress this term in our notation. 

4 -model. The full quantum action of the 4 model (|15"|) is the Lagrangian expressed in 
unrenormalized parameters 

U = l(d<f>o)-^(<f>l-^) 2 . (130) 

Since in two dimensions only the mass receives a divergent contribution (see below, "educated 
guessing") it is enough to choose a minimal renormalization scheme: 

SZ = =^ 0o = 0, 5X = => A = A, $ = /j 2 + 5/i 2 . 
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Since we have to evaluate the path integral (|126|) in the vacuum sector, and to fulfill the 
asymptotic boundary conditions ( |127| ), we expand the Lagrangian (action) around one of the 
classical vacua (pT9|) . We choose 0y = ^= so that = + 77. Thus by this perturbation 
theory one will never "see" soliton contributions. Also inserting [Aq = li 2 + 6fi 2 we get0 

C A = \[(d V f - 2 M V] - vV/ 3 - \^ - \^ 2 (V + ^=77) + 0(ft 2 ). (131) 

We have put the constant into the higher order terms indicated by 0{h 2 ), since we will 

determine 5fx 2 only in one loop (H) order. The physical boson mass at tree level is m = 
as one can read off the quadratic part of the Lagrangian, and has the correct sign. The vacuum 
boundary condition (pole prescription) in (|129| ) is respected implicitly by a small imaginary 
part of the squared mass, i.e. m 2 = m 2 — ie. With this Lagrangian we get for the generating 
functional (|129[) , but with the source coupled to the "physical" field 77 

Z[j] = J Vr] ekS^iUdvf-m^+jv+Cr) (132) 

e f / dj*($ l(a v ) 2 -mW\+fn) _ ^ 33) 



_k=o 

This is the perturbative expansion, where Cj is the "quantum" interaction Lagrangian 

C, = -VX^v 3 - \^ - \^ 2 (r, 2 + ^ + 0(h 2 ). 

The perturbative expansion can also be written as functional derivatives of the free generating 
functional 

Z\j) = eM^'rtftZoW with Z [j] = fvn e t/^(i[W— V]+i,) ( (134) 

J A 

where the index A indicates that the free generating functional has to be evaluated in a 
regularized way. The regularization takes place in the set of considered fluctuations 77, i.e. the 
"path integration domain" PID. We consider two possibilities which are very similar in the 
trivial sector. 

Energy- momentum cutoff. We restrict the PID to a regularized one, which is characterized 
by the cutoff A or the projector, respectively, 

/dk 2 
—eiA-ihDe^, 

which is the (regularized) unit in the regularized domain PlD^ eg = {r/ \ 5 ■ 77 = r/}. Also the 
sources j(x) must be functions in this domain, i.e. f dy 2 5(x — y)j(y) = j(x). The regularized 
set of fluctuations PlD^ eg are functions which have a compact support in the spatial Fourier 
transformed variable, i.e. 



2 

ik-x 



/dk 2 
^(A-1^1)^) e 



22 Since we are considering an unbounded space-time, no surface terms contribute to the action in this 
expansion. 
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The free generating functional (|134j ) can be evaluated in several ways. One is to expand the 
fluctuations r\ around the configuration rj , i.e. 77 — ^ 770 + ?7, where r] fulfills: 

(□ + m 2 )r/ (x) = j(x) (135) 
=> r] (x) = - J dy 2 A r F eg (x - y)j(y) with (□ + m 2 )A r ^ 9 {x) = 8(x). (136) 



The differential equation for the free field in (|135|) is solved by the method of Green func 



tions. The Green function A r p 9 (x) in ( |136| ) is called Feynman propagator and is characterized 
by the vacuum boundary conditions which is encoded in the pole prescription. Its Fourier 
representation is given by 



/JU2 -ik-x 



We have written the pole prescription explicitly. Inserting (|136|) into the free generating 
functional (|134|) one obtains 

Z [j] = N e~rAl d - 2d y 2 H*)*F'(*-v)Hv) i ( 137 ) 



where the constant N is a number which will be canceled by the denominator in ( |128|) . From 



(|134|) and with the form of the interaction Lagrangian (|131|) a general n-point function is a 
composition of Feynman graphs which are given by: 



z 



y ihA(x — y) propagator 
J dz 2 3- vertex 

— I A / dz 2 4- vertex 



T 



z —j:5n 2 fdz 2 seagull-counter-term 

j i -^5fi 2 J dz 2 tadpole-counter-term 



Because of the spontaneous symmetry breaking also a three-vertex interaction occurs. The 
integrations for the vertices lead to momentum conservation in momentum space. The discon- 
nected graphs are canceled by the denominator Z[0] in ( |128| ). This gives the usual Feynman 
rules (see for example (3l|). The order in h of an amputated graph is as follows (external 
lines are replaced by wave functions in the S-matrix) p^] : Because of the expansion of the 
exponential function in (|133| ), each vertex V contributes a factor h~ x . From each internal 



line J, ~ h 2 fz\fjAj, comes a factor h. Therefore \W\j] = 0(^ y ), where Z[j] = e^ w ^ . 
The functional W[j] is the generating functional for connected Green functions. Therefore the 
contribution of an connected amputated graph to W[j) is proportional to 

h T - v+1 = h\ 

where t is the loop number of this graph. 

Renormalization conditions. With the Feynman rules the two-point function has the 
following graphical representation: 
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G(x,y) 



I 



0(h 2 ) (138) 



In two dimensions only the seagull (third graph) and the tadpole-pole (fifth) graphs are diver- 
gent if the regularization is removed. We fix S/i 2 by the requirement that the tadpole vanishes 
exactly. This means 



= 



=> ih3\^XfxA(0) + -^Sfx 2 = (139) 
V A 

Thus the mass-counter term is given by 

,2 .«=oxa/^ f dk2 0(A-\k x \) 3Xh f A dk , . _ v 

Sn 2 = -ih3XA(0) = -ih3X / - — -— - = / - (140) 

^ 1 ^ 7 (2vr) 2 k 2 - m 2 + 2tt i Vfc 2 + m 2 1 J 

where in the last step we have carried out the fc -integration and set k 1 = k. This is the one 
loop counter term and from (|140|) one can see that it is of order 0(h) and depends on the 
regularization cutoff A. By this renormalization condition also the seagull-self-energy graph is 



canceled. The sum of the seagull- and the mass-counter- graph (graph three and four in |138|) 
reads as 

(ih) 2 3XA(0)+th5^ 2 = (ih) 2 3XA(0) - (ih) 2 3XA(0) = 0. 

This renormalization scheme is the most simple one from the technical point of view, but one 
must not forget that the second graph in ( |138| ), a three vertices self energy diagram, gives a 
finite contribution to the pole of the one loop propagator. Thus the renormalized parameter 
m, defined by the renormalization condition (|139|) is not the pole of the propagator. For the 
physical one-loop mass, i.e. the pole of the propagator, one must take this finite contribution 
into account. Summing up the series in (|138|) as usual |3l| the loop of the three vertex self 



energy diagram contributes to the pole. Thus the physical one loop mass is given by |39| 

9 2 \ i- f* dk 2 m 2 2 v^,. 

m P = m + 9Xih / jtt^jt^ 2 . ■ w7j ^ 2—— t p^ m 2 = m - —h\. 

J {2tt) z (k 2 — m 2 + ie)((k — p) z — m z + le) 2 

We also could have renormalized the coupling, e.g. through low-energy limit of scattering 
amplitudes, i.e. an analogous renormalization condition as in ( 3.2.2| ), to get simpler relations 



between the renormalized parameters and physical measured quantities. But our renormal- 
ization condition ( |139| ) is an equivalent scheme and for concrete calculations is the most 
comfortable one. 



The renormalization condition (|139|) gives the explicit expression for the counter term param- 
eter Sfx 2 and defines the renormalized parameter // resp. m = a/2/i. Its concrete value must be 
determined by experiment of course. But all quantities expressed through this parameter get 
their explicit meaning by the renormalization condition. With the explicit expression for 8jjL 2 
and the definition of fi through the renormalization condition also the one-loop renormalized 
Lagrangian (|130|) has an explicit meaning (<p = 0, Ao = A, /ijj = /i 2 + 5fi 2 : 

£ A = C(fi) + SC A , 



56 



where £(//) is the classical Lagrangian expressed through the renormalized parameter, and 
5C\is the (one-loop) counter term Lagrangian, the quantum corrections, which in our scheme 
is given by 

5£ A =^(V-C)+0(fc 2 )- (141) 



2 V A, 

The classical action can be thought of as renormalized at zero loop level, with the renormal- 
ization condition 5fi 2 = 0. Therefore all classical quantities are expressed through the finite, 
regularization independent, parameter /z 2 . This parameter must of course also be determined 
by experiment, for example by classical scattering amplitudes of solitons. In general this is a 
hypothetic issue, since the classical meaning of a quantum field theory is not always evident 
(see for example fermions). 

Mode number cutoff. The regularized evaluation of the fundamental ingredient of the 
perturbation theory, the free generating functional Z [j] and thus the Feynman propagator 
Ap 9 , can also be carried out in a discrete manner. For this one compactifies the spatial 
dimension to a circle of (large) perimeter L or considers a compact interval of length L and 
introduces appropriate boundary conditions (periodic or antiperiodic are proper ones, see 
below). The path integration domain PIDn eg is then defined by a finite Fourier expansion 
according to the discrete Fourier modes. In the vacuum sector there is no big difference between 
an energy momentum cutoff (EMC) and a mode number cutoff (MNC). Only the finite spatial 
momentum integrations are replaced by finite sums, i.e. one has the correspondence 

r A dk A 

/ 7rf( k ) M Ysfikn) with Lk n = (2n + A)n. 

J -A 27T ^ 

The value of A = 0, 1 depends on the boundary conditions. The two expressions coincide up 
to order O(t) and the two cutoffs EMC and MNC are related as 

These relations between EMC-integrals and MNC-sums are not so simple in the soliton sector 
as we shall see. But the counter term parameter 5fi 2 (|14Q ), defined in the vacuum sector is 
simply 

r 2 o^l 1 1 3XH f A dk A 1- 

5u 2 = 3\h — > = / +Q(-,-). 

2L ^ A) V k l + m2 2rc 'o Vk 2 + m 2 V L 'A ; 

Therefore we will always write the counter term in the integral representation, independent of 
the used regularization scheme. 

Renormalization of the sine-Gordon model. Since we were so explicit above we can treat 
the S'G-model relatively quickly. The full quantum Lagrangian, i.e. expressed in unrenormal- 
ized parameters, is given by 

^A = ^(^) 2 + -[COS( V ^0)-1], 

2 7 
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where we have used that only the mass gets divergent contributions (if the regularization is 
removed), as we shall see. Thus we have the following minimal renormalization scheme: 

5Z = 0, 5 7 = 0, ,t4 = /i 2 + 5fx 2 . 

Expanding the Lagrangian around the vacuum = and inserting /Xq = ji 2 + 5ji 2 one obtains 

U = ~ [(d v ) 2 - /iV] + ^ - ^ ■■■- 6 -^v 2 + .... (142) 

Only the seagull loops, which are for example of the form 




are divergent [3£|, and thus our minimal renormalization scheme is enough to get rid of the 
divergences. The tree-level boson mass is fi as one can see from the quadratic part of (|142|) . 
At one-loop the renormalization condition, that the seagull loops vanishes reads as 

t^-f + / ■ / = 



o fc ^7 f dk 1 

^v/p + zx 2 ' 



8fx 2 = h 



Since at one-loop order the one-loop seagull graph is the only contribution, the mass parameter 
/i, defined by the above renormalization condition, coincides at one loop with the pole of the 
propagator, i.e. the physical mass of the boson. The one loop counter term Lagrangian is 
given by 



,2 



5C A = -ii- [cos(v/70) - 1] + 0(h 2 ). (143) 

7 

The renormalization condition that all seagull loops vanish gives the complete counter term 
as follows [0] 

5£a = (e^ 2 - l) [cos(VT0) " 1] • 
3.2.4 Quantum action for solitons 

With the counter term contributions 5C\ and the explicit expressions the counter term pa- 
rameters <5/i we can write down the renormalized (quantum) action in the soliton sector. The 



classical kink masses (|22| , |56D and all other quantities are expressed through the renormalized 
parameters /i, which are defined by the renormalization conditions from above. It is important 
to use the action renormalized in the vacuum sector, also to calculate the desired quantities 
in the soliton sector, since this ensures that if one compares quantities of the two sectors one 
talks about the same things, i.e. one uses the same renormalized parameters, defined by the 
same renormalization conditions. The one-loop mass counter terms for the SG- and 4 - model, 
as calculated above are given by 

2 f^ m2 f A ^ (144) 
. 2 ,3A f A dk nA . 

^ = W Tra' (145) 
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where m is the tree-level boson mass, and related to the renormalized mass parameter /i by 
(I = 1/2 for SG/0 4 )0 

m = \fl\x. 

The one-loop counter term Lagrangians are given by (|143|) , (|141|) 

5C SG (5fi) = -H- [cos(^) - 1] + 0(h 2 ) (146) 

6C,(6») = ^^-^J+0(h% (147) 

these counter term Lagrangians lead to additional contributions to the action also in the kink 
sector. To calculate these contributions we expand ( |146|) ,( |i~47j) around the kink solutions 
(20,35) of the two models (cr = 1/ — 1 for kink/antikink), 

m(x — Xq) 



TTl 

W = cr , tanh 

V2X 



(148) 



SG: (t) K = arctanfe^-^)], (149) 

and keep terms up to order 0(h). For the SG kink we have re-introduced the original co- 
ordinate and field relatively to (|35|). For nontrivial classical solutions the classical action is 
non-zero, in contrast to the vacuum solutions, and therefore already gives contributions of 
order 0(h°). Thus the quadratic fluctuations (the semi-classical amplitude) are already of or- 
der 0(h). So the nontrivial contribution of the classical action S[4>k] shifts the order by one, 
relatively to the vacuum sector, in the expansion of quantum fluctuations. Thus the one-loop 
counter term contribution in the kink sector is given by 

5£(4> K + ri; 5fi) = SC{<f> K ) 6fj) + 0(h 2 ). 

The interaction of the fluctuations with the counter term ^5/1— 0(h) give higher order 
contributions. Thus the one-loop quantum corrections to the action in the kink sector are 
given by 

5S SG (5fi) = — f dtdx [cos(V7</0 - 1] + 0(h 2 ) = -T—5ti 2 + 0(h 2 ) 
SS^Sfi) = ^ J^dtdx - |Q + 0(h>) = -T^V + 0(tf). 



For SG we have used that SCsg = —^tU(4>) and the mass formula for static solutions. In 



the SPA calculation of the energy spectrum this gives additional contributions to the mass 
of static solutions. In addition to the classical mass of static solutions (|101|) , such as kinks, 
the counter terms give, besides the quantum fluctuations, contributions to the quantum mass. 
For SG and 4 kinks this gives 

S cl (<t> K ) + SS^k) = -T (M d + 5M(5fi)) with : (150) 

SG : 5M SG = Sfi 2 (151) 

771 

4 : 6M^ = -6fi 2 . (152) 
A 



23 



For the SG-model it is even the one-loop physical (pole) boson mass. 
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Now we have all ingredients to calculate the one-loop quantum corrections to the classical 
masses of static solitons like the SG and 4 kinks. 



4 Quantum masses of static solitons 

We have now all ingredients to calculate the quantum corrections to the classical kink masses 
(|22| , |36|) in a semi-classical approximation, i.e. at one loop order. The one loop corrections 
to the masses, especially for the supersymmetric extension of the here considered models (see 
below), are of particular interest, since they are connected with the possible occurrence of 
an anomaly in the supersymmetry algebra. The main question is whether the Bogomolnyi 
bound stays saturated by N = 1 supersymmetric solitons including quantum corrections resp. 



whether there exists an anomaly in the central charge of the SUSY algebra (||39|, [pq| ,p6|). 
This can be decided already at one-loop level. 

As one can see from (|109|) the ground state energies include for both, the vacuum and the 
kink sector, divergent sums over the mode energies uf. This corresponds to the ambiguity 
(freedom) in the choice of the energy zero point in ordinary (no gravity) quantum field theory 
(only energy differences are measurable). The absolute energy-zero point must be fixed for the 
vacuum ground-state. The natural choice is to normalize vacuum ground-state to zero energy. 
Thus we have to subtract the zero-point energy of the vacuum (no state is excited) from the 
energy calculated for any state. As long as one considers only the (undistorted) vacuum sector 
this procedure is trivial and respected by the normal ordered Hamiltonian, i.e. 

:W:=W-<0|W|0) (0|:W: |0) = 0. 

Thus the vacuum ground-state |0) has zero energy. The same reference-point for the energy 
must be used for all other states, also for the solitonQ, i.e. 

(sol\ : H : \sol) = (sol\ H \sol) - (0| H |0) . 

To evaluate this difference in the presence of a nontrivial background such as a kink is a 



highly nontrivial issue and object of controversial discussions for years (see |39|] and references 
therein). Thus the one-loop kink ground-state energy is given by 

E K — E v , 



where the ground-state energies Ex,Ey are given by (|109|) . In the discrete version, as given 



in (|109|) , one has to evaluate the difference of two infinite, divergent sums. Thus they (the 
theory) must be regularized in a consistent manner. This will be the subject of the following 
sections. The calculation of the difference of these ground-state energies is very similar to the 
Casimir effect but much more involved, since the kink background is of course much more 
complicated than conducting plates, which set up certain boundary conditions. 



4.1 Renormalized spectral function 

To evaluate the spectral function ( |102| ) in one-loop order one has to do a semi-classical expan- 
sion of the one-loop renormalized action^ S = J TxL dtdx (C(m) + 5C(6fjL)) around the classical 

24 That with the soliton also a quantum state is associated will be discussed below. 

25 Here we introduced the interval L for the spatial integration. The range of L depends on the regularization 
scheme and will be finite or equal K 
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solution, i.e. (f>(x,t) = <fi c i + n(x,t) up to order 0(h). The classical quantities are expressed in 
terms of renormalized parameters. For static solutions C / one obtains 

S[0 K + r](x, t)} = - (M(</> K ) + 5M(6fi)) T (153) 

~ J dtdx r](x, t) (□ + U"(<p d )) n(x, t) + 0(h 2 ). (154) 

TxL 



The first two terms in fll54j) are the classical energy (i.e. the classical mass for static solu- 
tions) and the counter-term contribution and given by the renormalized Lagrangian with the 
counter terms (|141] , |143| ). In the vacuum sector (0 C / = (py = const) both terms vanish. In the 
kink sector one obtains the classical kink mass M d and the counter-term contribution to the 



quantum correction of the mass (|151| , |152|) 



45m 2 ^ f A dz 1 . v 

SM SG = = 2h ^^== 155 

m5m 2 ^ f A dz 1 . 

In the integrals (|156|) we have transformed to the variable z = ^ so that the now dimensionless 
cutoff A is large (or small) relative to the scale m. The linear term in ( |154| ) is absent since 
4>d is a classical solution. In dimensionless spatial coordinates z = yx the spatial part of the 
operator in (|154j ) for the vacuum and the kink sector, respectively, is given as (/ = 1/2 for 



SG/^) 

Vacuum : Vl (z) = (-d 2 + I 2 ), (157) 

Kink : OkXz) = (-d 2 z +l 2 - ^ . (158) 

V cosh z J 

Therefore in the kink sector one has exactly the stability operator (see appendix). To evaluate 
the spectral function one has to diagonalize these two operators and one is left with the 
following expression 

Tr s [e-f, HT ] = e-fA M ^ s ^ T V Va V V e tIl + (h 2 ) (159) 

J JVaM 



{PID reg } 



The exponent of the first factor in (|159|) is only non-zero in the kink sector. The only difference 



to (|102|) is the counter term contribution. With {PID Reg } (for "path integration domain") we 
indicate that the set of considered paths and therefore the spectral function depends on the 
regularization which will be used. This will be discussed in detail in the next sections. In the 
quadratic part of the action the field- degrees of freedom are the fluctuations n, the classical 
solution (j) c i is a fixed background which is nontrivial in the kink sector. 



4.2 Mode regularization (MNC) for bosonic kinks 

Mode regularization proceeds by making the system discrete by a finite volume L and thus 
countable (so one is not involved with functional analytical subtleties as in the continuum) 

26 In principle there are also surface terms because of the finite time interval T (see above). But at the end 
they are transformed back to total derivative to get the harmonic oscillator action. This is a trivial step and 
we suppress it here. 
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and introduces a mode number cutoff (MNC) to regularize otherwise divergent expressions. 
This is very analogous to a lattice regularization and, at least for bosons, there exists a one- 
to-one correspondence even though both schemes are not related by an ordinary coordinate 



transformation in the path integral [24]]. There a two main points in a MNC-scheme that 
become crucial in the presence of a nontrivial background: (I) the boundary conditions should 
not induce effects that do not vanish in the limit L — > oo. (II) The number of discrete states 
and the correct evaluation of the sums over mode energies. The relation between the "cutoffs" 
in different topological sectors is given by the requirement that in both sectors an equal number 
of modes is taken into account. Thus the two Hilbert spaces PID]^ and PIDf™ k have the 
same dimension and therefore this two spaces are isomorphic. Thus the correct counting of 
the discrete states is essential for mode regularization. 



(I) Boundary Conditions 

In those cases where boundary conditions are essential in the regularization process, we adopt 
a new principle, which is closely related to but less restrictive than the "topological boundary 
conditions" of |^]. To ensure that the BC do not introduce a force which contributes to 
the energy we compactify the spatial direction to a circle of perimeter L (L — yL in our 
coordinates). The fields must therefore fulfill a matching condition which leads to certain BC, 
depending on the topology of the line bundle one chooses. According to ( |159| ) the fields that 
have to fulfill this matching condition are the fluctuations 77. First of all the resulting BC must 
be a linear relation so that the considered set of paths in ( |159| ) form a linear space. We require 
now that the BC must be chosen in a way that the transport of the quadratic Lagrangian 
£W(rj) around the compactified dimension leaves C^ 2 \rj) invariant. This means 

z -> z + L : £ (2) (77) -> £ (2) (77) 5C {2) (77) = 

Otherwise the action would get an additional contribution ~ dzSC^ by the spatial integra- 
tion. Thus the topology of the line bundle on which 77 lives must be chosen in a way so that a 
surrounding of the compact dimension (z — > z + L) induces a linear symmetry transformation 
of 5C^ 2 \ In both topological sectors the influence of the classical solution ( |158| ) is symmetric, 
thus on the circle one has Os(z + L) = Os(z). Therefore in both sectors one can use the 
following line bundles: 

z -> z + L rj(z + L) = (-) A ri(z) A = 0, 1 

The values of A = 0, 1 corresponds to periodic P and antiperiodic AP BC. This is the Z 2 
symmetry which is despite spontaneous symmetry breaking conserved in the quadratic part 
of the Lagrangian. The BC can be chosen independently in both sectors, all combinations 
(Vac\kink) = (P, AP\P, AP) are allowed since our symmetry principle ensures that no contri- 
butions due to the BC occur. There is no need to use common BC in both sectors (in contrary 



to @,|3J]). This is the physical principle for mode regularization. It is rather simple and not 



restricted to two dimensional theories. In the case of fermions it will become more exciting.^ 

27 Also homogeneous BC are allowed is by the symmetry principle although they are not topological. For 
the sake of simplicity we do not consider them here although in principle possible. 
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4.2.1 Regularized kink mass 



To carry out the path integral in (|159|) one has to diagonalize the quadratic action. In a 



mode regularization scheme this is done by a finite expansion of the fluctuations according to 
eigen- functions of the operator Os(z). For this one has to solve the eigen- value problem 

(-c\ 2 + O s (z)) £„ = co% (160) 

with proper boundary conditions. A finite expansion of the fluctuations according to the 
eigen-functions of (|160| ) 

M+ 

v{x,t)=Y^c n (t)Z n (x) (161) 

-Af- 

leads to a finite, countable set of harmonic oscillators with eigen-frequencies u n (see above). 
The mode number cutoffs M_, M + will be determined for each special case below. Therefore 
the path integral in ( |159| ) reduces to a finite product of harmonic oscillators of frequency u n : 

II dc a l B n (T)Vc n (t)erJT^l-^ 2 J (1 62 ) 

M ~ J JCn„,t' 



The measure B n {T) is the same for each oscillator n, independent of the sector: the measure of 
a harmonic oscillator. Thus for an equal number of modes in the vacuum- and kink sector one 
has the same measure in both sectors. There is a subtlety connected with possible zero modes. 
For the zero mode integration one has to use collective coordinates. In the purely bosonic case 
this is a fairly trivial thing and connected with breaking of the translation invariance by a 
given kink position. As showed above, in the considered order one can omit the integration of 
the zero mode.0 

The residual integration is easily performed. Putting all things of (|159| ) together one obtains 
for the difference between the kink ground state energy and the vacuum ground state energy 
the following one loop kink mass 

h hi M+ N+ \ 

M K = E K -E v = M% + ^^ + %[Y, UJ n- J>? + 5M ^ 5m ) + °^ 2 )- ( 163 ) 

bound \ —M- -JV_ / 



The mode energies arf ' y are given by the eigen values of (|160|) for the kink (K) and vacuum 
(V) sector. 



28 In the supersymmetric case this issue is much more involved, since fermionic zero modes have the same 
origin and should be treated by a common collective coordinates (see below). 
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4.2.2 Vacuum contribution 



For the vacuum the calculations and results are rather simple and in both cases given by 
(A v = 0/1 for P/AP BC): 

BC — quantization Lmk n = (2n + A v )ir (164) 

hm I f (2n + A v )7r\ 2 , , 

groundstate energy : E v = \ ( r J + 1 (165) 

mode number : #y = 2N + 1 + Ay (166) 

,4 (2N + Ay)ll , AyT[ 

energy cutoff : A A = ki = ± T 7 := A + V—. 167 

Lm Lm 



In the vacuum the conversion of the sum (|167|) into an integral is straightforward. Nevertheless 



we use the Euler-MacLaurin formula which will be the appropriate technique in the kink sector. 
So one obtains 



Lm 



Ey= h fhl-Ay)+2 J |W^^+(l + VA^+T) +0(ii). (168) 



Lm 



In the last term within braces, the "surface term", we have already carried out the limit L — > oo. 
4.2.3 Kink sector 

In the kink sector all calculations become much more involved. First we outline the general 
principles. In the following subsections we carry out the calculation for the considered models. 
A very sensitive point is the evaluation of the (potentially) infinite sum. The appropriate and 
mathematically exact tool to do this is the Euler-MacLaurin formula, which is given in the 
appendix (|318j ). 

In both cases the continuum states ( |314| ),( |3T5| ) are asymptotically of the form 



£(q,z — ±oo) = Z ± (q) e iqz = ^le^+^l (169) 

where Z ± (q) are complex valued functions of the momentum q. For SG and 4 , respectively, 
they have the explicit form 

Z±{q) = ±l-iq (170) 

Z%{q) = {2-q 2 )Ti'Zq- (171) 

The absolute values \Z ± \ are not interesting (can be absorbed in the normalization) but the 
argument functions arg[Z ± (g)] =: ^(q) will become very important. Their explicit forms 
depend in a crucial way on the position of the branch cut chosen for the argument function. 

From the asymptotic form ( |169| ) one can see that going once around the (large) space-circle 
one picks up a total phase 

ql + [<p + {q) - <p-(q)] =: q~L + 6(q). (172) 
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Figure 10: The Sine-Gordon scattering phase 5(q) for different branch cuts, (a) cut = Im + , avg(z) G 
(-TT,f]- ( b ) cut ^ Im±, &vg(z) <E (a,a + 2ir], (c) cut = Jm_, arg(z) e (-f, 

where L is the perimeter in z-coordinates from above, i.e. L = yL . The analytical structure 
of the scattering phase 5(q) := V 9+ (?) — V 9 (?) an d its asymptotic values depend on the position 
of the branch cut. By setting up boundary conditions the momenta get quantized as follows 

q n L + 5{q n ) = (2n + A k )tt (173) 

This is a transcendental equation for the allowed g n 's and we will solve them by iteration . 
The constant Ak again determines the kind of boundary conditions {Ak = 0/1 for P\AP 
BC). To apply the Euler-MacLaurin formula ( |318| ) to evaluate the ground state energy 

M+ 



we have to know the addend as an explicit function of n. Therefore we have to resolve (|173|) to 
get q n as an explicit function of n, at least up to sufficient orders in L and the mode number 
cutoff M±. This is done by an iteration, where the second step gives 

L L L L A 

That this is a reasonable approximation is guaranteed by the Banach fixed point theorem. The 
scattering phase 5 is bounded for all values of n (its maximum range is 47r as the difference 
of two angles (|172|) ) and therefore the iteration ( |174[) is a contraction if L > 5 max . This is of 
course true since we are interested in the limit L — > oo. Next we do the explicit calculations 
for the SG and </> 4 model. 



4.2.4 sine-Gordon-kink 



From the asymptotic states ( |314| ) one obtains for the scattering phase 



<K<?) = V 9+ (?) — V 9 (?) = — 2arctang + Cut 

where Cut stands for branch cut position dependent constants. The scattering phase varies in 
range of 2n. Its branch cut dependent shape is given in fig{10]. With the position of the branch 
cut also the discontinuity moves. At the discontinuity the phase jumps by 2n and this is also 
the vertical distance between two neighboring straight lines (2n + A k )ti — Lq n which must 



intersect the graph S for q n being a solution of (|173|) . Therefore in the case of a discontinuous 



phase always one mode has no solution and thus does not exist .0 

29 Note that at the discontinuity the phase takes only one value which corresponds to the semi-open intervals 
for the range of angles 
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From the continuous version of the scattering phase one can obtain by Levinson's theorem the 
number of discrete states as 

S(-oc) - 6(00) 

„ Jrdiscrete 
Z7T 

This is exactly the translational zero mode of the SG kink, which counts as a full mode.0 
Therefore the continuous spectrum in the kink sector is shifted down by one mode relative 
to the vacuum sector. For a discontinuous phase it is a low lying mode (in the sense that 
M± — * 00) which has no solution, whereas the phase goes asymptotically to zero (5(±oo) = 
0). For a continuous phase one has to omit one of the high modes, but now the phase 
takes asymptotically the finite value 2ti. This interplay between the asymptotic values of the 
scattering phase and the omitted modes explains that even though a mode at the threshold of 
the continuum states becomes a "bound" state it is possible to subtract a high or a low mode 
of the continuum. These are the general rules which apply in all considered cases, also for 
fermions. 

We calculate the kink mass for the discontinuous phase with the branch cut at WL + , so that 
the phase jumps at q = and is symmetric.^] For a more general branch cut position the 
calculation is quite analogous but one has to write a little bit more. The case of a continuous 
phase will be considered in the 4 model. We leave the combination of boundary conditions 
of the two sectors arbitrary. For the sum over "continuous" kink mode energies we make a 
symmetric ansatz 

N 

Therefore the "continuous" mode numbers are given as follows 

vacuum : # y = 2iV + 1 + Ay (175) 
kink : # K = 2N + A K . (176) 

Now independent of the BC one must have #y — #x = 1 due to the discrete zero mode in 
the kink sector. For equal BC (Ay = Ak) one can see that this is already fulfilled. In the 
case of different BC one has to add {Ak = and Ay = 1) or subtract (Ak = 1 and Ay = 0) 
in addition one of the high modes. So different BC compensate the effect of the discontinuity 
of the scattering phase. Thus for the kink modes we have 

N 

Vti + l + {A v - A K )my/AF+l. 

~(N+A K ),^0 

For the additional high mode we have already taken the limit L — > 00, but a detailed calcula- 
tion shows that the result is independent of the sequence of the limits L — > 00 and A — > 00. 
Inserting the iterative solution (|174j) and using the Euler-MacLaurin formula we get with the 

30 This is not as obvious as it seems. This will become clear in the case of fermions. 

31 We choose the semi-open interval for angles so that the phase takes the negative value at q — 0. Therefore 
in both cases P/AP BC the mode n = has no solution. 
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variable transformation z 



(2n+A K )n 
Lm 



N 

E 

-(N+A K ),^0 



N 



.K 



mA K + 2j2^n + (Av - A K )my/A 2 + 1 



mAx + 2m 



2tt' 



-Lm 



(2 + Ay)7T 

Lm 




Lm 



8(z] 



+ 1 



+m (l + VA 2 + l) + (A y - ^)mVA 2 + 1 
In the surface term we have already carried out the limit L — > oo . As a next step we expand 



the root y z 2 — j^5(z) + O(-p) + 1, where jj^S{z) is, for all z, a small quantity if L is large 
enough, since 5 goes to zero sufficiently fast. But for branch cut positions where the phase 
takes finite asymptotic values one has to take care at this point. Let us consider such an 
example for the 4 model. Putting all together we get with (|168|) and (|156[) for the kink mass 



/ ,<m I .i-n 



fvtn 

M= M d + — (A K - (1 - A v )) + hm 



(2+A k )tt 



(2-A v )tt 



hm 



{A 



v 



A K )VA 2 + 1- hm (J^ 



Lm 



dz zS(z) 



2tt^TT 



— LmV z 2 + 1 
2tt 



dz 1 
vr y/z 2 + 1 



0( 



0(h 2 



In the second line we have already carried out L — > oo since the integrand is independent of L. 
The surface terms of the Euler-MacLaurin formula have canceled each other. This is always 
the case, so we do not write them down in the further calculations. Doing the integrals and 
taking the limit L — >• oo one obtains 



+ — \A K + A V -1- (A K + A v ) + (A K - A V )VA 2 + 1 + {A v - A K )\J A 2 + 1 



—hm — (2 — 7r). 
2n 

So the divergences cancel each other nicely and we finally obtain for the kink mass 

M K = M cl 



^ + 0(h>) 
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We have been so explicit to show that even in the case of different BC the correct mode 
counting gives the correct and finite result. The calculations for the continuous scattering 
phase are quite analogous and give exactly the same result. 



4.2.5 4 -kink 

In the 4 model everything is straightforward, but more involved since 4 - kink is one "degree" 
higher than the SG - kinkF]. From the asymptotic states (|315|) one obtains for the scattering 

32 SG- and cj) 4 - kink are two special cases of a class of kinks whose zero mode is of the form ~ ~ iR^j - I n 
this sense SG/cf) 4 is of degree I — 1/2. 
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Figure 11: The <p 4 scattering phase for different branch cuts: (a) cut = R_, arg(z) S [— w, it), (b) 
cut ^ R_ V R+, arg(z) S [a, a + 2ir). (c) cut = R+, arg(z) £ [0, 27r), With the position of the branch 
cut the discontinuities are moving. 



phase 



3q_ 

2-<r 



5(q) = -2 arctan ( ) + Cut (177) 



Cut again stands for branch cut dependent contributions. The scattering phase (|177|) takes 



its values always in the semi open interval [— 2ir, 2ir) (if the upper or lower interval bound is 
the open one depends on the convention). For different branch cut positions the phase has 



the form as shown in fig.|ll|. For discontinuous phases two modes do not have a solution and 
therefore do not occur in the sum over mode energies. The information on the discrete modes 
is again encoded in its asymptotic values. Levinson's theorem gives for the number of discrete 
states 

6(-oo) -a(oo) 

rrdiscrete ^ 

111 

These are the zero mode and the excited bound state of the 4 kink ( |315| ). Therefore the 
continuous spectrum of the kink shifts down by two modes relative to the vacuum. 

We now calculate the kink mass using the continuous phase to show how to deal with the 
non-zero asymptotic values of the phase. We again do not fix the BC combination of the 
two sectors to show that there is no need to use the same BC in the two sectors. The above 
consideration ( |176| ) showed that the subtlety in mode counting only depends on the use of 
different or equal BC. Thus, to reduce the notational costs, we choose periodic BC for the 
vacuum and let the BC in the kink sector unspecified. For the sum over "continuous" kink 
mode energies we again make a symmetric ansatz 



JV-l 



-(N-1+A K ) 

Therefore the "continuum" mode numbers are given as follows 

vacuum : #y = 2N + 1 

kink : # x = 2N - 1 + A K 

Now independently of the BC in the kink sector one must have two continuum modes less than 
in the vacuum sector, i.e. #y ~ #x = 2. In the case of equal BC (here periodic, i.e. A K = 0) 
this is already fulfilled. For different BC, i.e antiperiodic BC in the kink sector (A K = 1) one 
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has to subtract in addition one of the high modes. Therefore with the use of (|315|) we get for 
the kink sum 



iV-l 



= -(1 - A K )m + m 2W(y ) 2 + 1 - A K 



Qn ■ 



(178) 



The (dimensionless) momenta q n are given by the iterative solution ( |174| ) and therefore the 
root in ( |178| ) can be written as (L = yL) 



(2n + A k )tt 
Lm 



1 s ^ 2n + AK)7r - 



Lm 



Lm 



+ 1 



(2n + 2 + A K )n 1 (2n + A k )tt 



Lm 



Lm 



[5(2- 



Lm 



+ 2tt] + 1 



where we have shifted j^- from the phase term to the first term. The reason for this is that 
now one can consistently expand the root in the following integral since is even for 

z^ooa small quantity. Without this shift by 2ti the approximation would break down and 
result in a divergence. 



Putting all together ( |156| ),( |167| ),( |315| ) and inserting in the Euler-MacLaurin formula one gets 
with the variable transformation z = ( 2n +^+ A ^) 7r . 

Lm 



M. 



K 



km . t . km 
M d + km— + — (A K -l)- — + hm 



( 



hm 



A K VA 2 +T-hm (J^ 



V 

dz z[5(2z) + 2tt] 

2^~jWTT~ 



■>TT 

Lm 



dz /-^ 

— Lmy z l + 1 
2tt 



+ 3 



dz 



o 27rv^TT 



o( z ) + o(n 2 ) 



In the second line we have already carried out L — > oo since the integrand is independent of L. 
The surface terms of Euler-MacLaurin formula have canceled each other. Doing the integrals 
and taking the limit L — ► oo one obtains 

M K = M d + hm^- + ^ (a k - 1 - 1 - A K + AkVA^+T - A K VA* + l) 

_ 1 . _ 1 

—nm — (3 — 2ix) — nm — — 
2vr l ; 2v/3 

So again the divergences cancel each other by correct mode counting, and one obtains, inde- 
pendently of the BC-combination, for the kink mass 



M K = M d + hm 



4VS 2vr 



+ 0{h 2 



Similarly, other choices of the branch cut positions lead to exactly the same result. This is as 
it should be, since the choice of a certain branch cut position is completely unphysical and a 
purely mathematical convention. So no calculations and results should depend on it. 
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4.3 Hilbert space of the soliton sector 



Since we now know the spectrum of the soliton sector (up to order h) we can consider the 
particle content of this sector. The energy spectrum in the kink sector is of the formF] 



E[{u l }] = M K + u d u d + J2 



where, in contrast to (|108|) , Mk is the one-loop quantum mass of the kink. Thus these are the 



excitations of the kink at rest. For the 4> A there also exists a discrete mode (the zero mode is 
not included). This is the spectrum if the kink is at rest. Like in the vacuum sector, to each 
mode {d, n} and their multiple excitation a state in the Hilbert (Fock) space is associated. 
Single excitations (v n = 1) correspond to the fundamental quanta (in the presence of the 
kink) or even to new particles {y d = 1 or the kink itself). This states are different from the 
vacuum states, as we will show. Thus in addition to the vacuum Hilbert space (vacuum and 
multi meson states) there exists a kink sector Hilbert space. In the vacuum sector, besides the 
vacuum state |0) only continuum states (plane wave eigen-functions) \k\, . . . k n ) exists, which 
correspond to the fundamental quanta of the theory. To perform localized (normalizable) 
particle states one has to built up wave packets of these states. In the kink sector besides 
the continuum modes there exists the kink, which has a localized energy density, and also a 
discrete mode (for 4 ), which is normalizable. Thus the kink sector Hilbert space consists of 
the following elements (particle states) 



1. The lowest state is the kink particle \P) with the momentum P and the energy E = 
V p2 + M K (see dTTD). 

2. The excited state of the kink \P*) (only for 4 ) of momentum P and energy E = 

y/P^ + (M K + UJ d y. 

3. The scattering states \P, k±, . . . , k n ) consisting of the kink particle and n mesons scat- 
tering of the kink with asymptotic momenta P,ki, . . .k n . 

4. The scattering states \P*,k\, . . . ,k n ) consisting of the excited kink and n mesons of 
asymptotic momenta P, ki, . . . k n . 



For the state \P*) one discrete mode u)\ = is excited (see appendix). Higher excitations 
of this mode (v\ > 1) will be unstable against decay into kink and meson since its energy 
v\uj\ = i/imy lie for v\ > 1 above the meson mass m. Note that with the zero mode no 
new state is associated, since it reflects only the collective motion of the kink and is therefore 
contained in the energy-momentum relation of the kink in point 1. 

We have already seen that the kinks are stable under small perturbations (all eigen-values of 
the stability equation are positive) and that the stability or the existence of non-dissipative 
solutions, respectively, is connected with the existence of a topological conservation law (there 
also exist localized finite energy solutions whose stability arises from ordinary conservation 
laws. These solutions are necessarily time-dependent, like the sine-Gordon breather [0]). We 
have also seen, that the different topological sectors are not connected. We show that this is 

33 In this section we use units so that H = 1. 
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also true in the quantized theory. In terms of Hilbert space this is expressed by the following 
two postulates pOU , fl3|: 

I. The kink sector Hilbert space is orthogonal to the vacuum Hilbert space, i.e. for all ampli- 
tudes 

(kink sector\vacuum sector) = 0. 
This can be shown as follows: The topological charge operator^ 



Q = / dx Jo = c 



oo, t) - 4>{-oo,t) , (179) 

is conserved in time since j7o is the zero component of the conserved current 

J^ = e^d v 4> with 3^ = 0. 

Since <\> is hermitian also Q is a hermitian operator and thus its eigen-states are orthogonal. 
In addition Q is conserved in time and translationally invariant (see ( |179| )) , i.e. it commutes 
with the energy-momentum operator V^. Thus, independent of the considered sector there 
exists a basis in the Hilbert space so that and Q have common eigenstates. Therefore 
the eigenvalues of Q are "good" quantum numbers and each state in the Hilbert space can be 
characterized by them. Thus the existence of conserved topological charge Q has, due to the 
existence of a conserved current, analogous consequences as usual Noether charges, following 
from continuous symmetries. Now the action of the field operator on vacuum states are of the 
form 

<j>(x, t) \ki, . . . k n ) = (fiy + r) v (x,t) \ki, . . . k n ) — > <p v + rapidly oscillating terms. 

x— >±oo 



The rapidly oscillating terms do not contribute to the charge Q (see (|29|)), so that for all states 
of the vacuum sector the topological quantum number Q = 0. Whereas the action on kink 
sector states is of the form IHOfl 



(f>(x, t) \P, ki, . . . k n ) = 4>k(x, t) + 7)k{%, t) \P, kx, ■ ■ ■ k n ) — > 0x(±oo, t) + rapidly oscillating terms. 

X— >±oo 

The rapidly oscillating terms in the kink sector are the same as in the vacuum sector, up to 
linear combinations (see appendix), and therefore also do not contribute to the charge. But 
the kink function 4>k gives a non-trivial contribution, so that all states in the kink sector have 
the topological charge Q — 1. Since states with different topological charge are orthonormal 
one obtains for the amplitudes 

(kink sector\vacuum sector) = (Q = 1|Q = 0) = 0. 

Since Q is a conserved operator the kink and the vacuum sector are not only orthogonal 
but cannot evolve into one another. Thus, despite all kink sector states, built around the 
local minimum <pK, have higher energy (in the weak coupling limit, since M c \ = O(j)) than 
any vacuum sector state, they do not decay into the vacuum vector states, built around the 
lower minimum, as expected a priori, purely from energetics. In the quantized theory this is 
apparently due to the existence of the conserved topological charge. 



3 1 



c is a normalization constant, see (26) 
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II. The different topological Hilbert space sectors are not connected by any localized operator. 
Consider any localized physical observable A(t). That is, let 

where a(x,t) is a local function of the field and its derivatives with finite spatial support at 
any given time t. Then the equal time commutator 

[A(t), Q(t)} = lim c{[A(t),4>(L,t)} - [A(t),4>(-L,t)}\ = 

because of the causality condition, i.e. that all commutators of space-like separated operators 
vanish. Thus, any such operator A cannot connect sectors with different topological charges 
Q. This suggests that Q is something like a super-selection quantum number, separating the 
kink sector from the vacuum sector. 



4.4 Continuum calculation (EMC) for bosonic kinks 



In [g^] it was shown that the widely used (see e.g. || |6fl) common strict energy-momentum- 



cutoff regularization (EMC) leads to results for the kink masses which differ from that obtained 



by a mode regularization. In |]39[ the EMC were identified as incorrect by comparing the 
calculated masses with exact results known for the SG breather solution. In recent works 
(|34||)||35|) a remedy has been suggested using an analogy to the Casimir effect, which in 
any physically realistic situation has a natural UV cutoff. There are however several reasons 
why we think that this solution is not satisfying: (i) In a discretized calculation this would 
impose the (as we have seen) unnecessary requirement of identical boundary conditions in 
the topologically distinct sectors; (ii) It also depends in a crucial way on the position of the 
branch cut of the scattering phase which is completely unphysical. The procedure works only 
for discontinuous phases which go asymptotically to zero. However, as we shall see, in the 
supersymmetric case this is impossible because of the presence of half-bound states. 

In our opinion, the deeper question behind this whole issue which has to be solved is how 
to regularize/renormalize two different topological sectors (sectors with trivial and non-trivial 
background) in a consistent way. Therefore one needs a principle that tells one how to regu- 
larize two sectors in the "same way" so that one can compare them in a consistent manner (the 
prescription of a common energy cutoff evidently does not achieve this). The principle must 
determine in particular how to relate regularization parameters (cutoffs) in the two sectors. 
In a mode regularization scheme this is done by mode counting, i.e by the requirement that 
both sectors below a certain energy/momentum have the same dimension in field-configuration 
space. In a continuum calculation, the dimension can be measured by the spectral density. 
But these spectral densities are the quantities to be determined. So one needs an independent 
principle which relates the cutoff in the kink sector to the cutoff in the vacuum sector. We 
show that the requirement that the regularized units in the two sectors can be matched will 
achieve this. 



Required accuracy 

In the discrete case one has two approximation parameters L and A which can be used in 
calculations. In the continuum calculation the only regularization parameter is the cutoff A. 
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This cutoff is fixed in the vacuum sector and by renormalization, i.e. the counter-terms defined 
by renormalization conditions. The important thing is, that the theory should be renormalized 
only once to be consistent. Therefore all other regularization quantities especially in the 
nontrivial sector must be given by A in a unique way. These relations and other approximations 
must be accurate to sufficiently high orders, so that no finite contributions are lost and finite 
errors survive. For calculating the energy they are given as follows: 

Cutoffs: The integration boundaries in the soliton sector, must be correct up to order O(-p-) 
since the mode energies, i.e the integrand are of order 0(A) for high momenta and therefore 
orders O(j) in the cutoff give finite contributions even in the limit A — > oo. 

Spectral densities: The spectral density p = pkink — Pvac measures the (difference of the) 
number of states in the continuum. Since in the kink sector, there are also discrete states, its 
integral should give the negative number of these discrete states. The integral of the spectral 
density must give the correct number of discrete states up to order O(^) so that the error in 
the number of continuum states does not contribute to the energy. An error of order 0(j) 
would result in a finite error in the energy since the wrongly counted modes are multiplied 
with the mode energies which are at the high end of order 0(A). This is the analogue of mode 
counting in the discrete case where it is of course much simpler. 



4.4.1 Vacuum 



To calculate the ground state energy of the vacuum fluctuations one has to path-integrate the 
the quadratic part of the action 

1 f , , , J 2 S 



dtdxi](x, t)— \<j, v r}(x, t) (180) 

TxR 

Therefore one diagonalizes the operator in ( |180| ) as follows (z — ^p, I — 1, 2 for SG, 4 ) 

(-d 2 + l 2 )Uk, z) = cu 2 v (k)Uk, z) (181) 

£v{k,z) = -^e ikz u 2 (k) = k 2 + l 2 (182) 
V 2n 

The energy of a mode is then given as 

hm , , . hm 



E v (k) = —u(k) = -^Vk^TP (183) 

To diagonalize the quadratic part of the action one has to expand the quantum fluctuations 
r](z,t) according to the eigen-functions (|182|) . Here the cutoff regularization takes place as 
follows 

*"=/ ite(A - l «* (M) (184) 

After spatial integration in (|180|) the quadratic part of the action is a continuous set of har- 
monic oscillators a(k,t) with the energies ( |183|) . The continuous set is strictly cut off at the 
momentum \k\ = A. This cutoff characterizes the set of fluctuations which are considered in 
the path integration ( |159| ). In this sense regularization means to restrict the path integration 
on a subset of the continuous functions 

PID = {0[MxT] | v (x,lf) = rf(x,f) = Va(x)} — > PID Reg 
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The regularized subset of the considered functions is characterized by the eigen- functions which 
are taken into account in ( p. 84 ). The subset PlD Reg is obtained by the action of a projection 
operator in the whole space PID which is the identity in PID reg and we therefore call the 
regularized unit: 

5 v (z -z')= [ dk<d(A - \k\)C v (k, z')£ v (k, z) = — [ A dke ik{z ~ z,) = l s W A (z-z')] (185) 

J 27T J_ A 7T [Z- Z') 

The proof that this is a projector and it is the unit in PID Reg is straightforward 

8y = J dy5 v (z - y)S(y - z) = 5 v (z - z) = S v 
5vV — J dz'5y(z — z')r)(z', t) = i], 

where in the second line we have inserted (|184|) for rj. There are two interesting limits of ( [L85|) 



lim 5 v (z - z') = 5(z - z') (186) 

A^oo 

\iuv5 v (z - z) = - (187) 



7T 



The first line is obvious since in this limit ( |185| ) is a representation of the Dirac delta- 
distribution; the second limit provides the diagonal elements which will be needed later. 

4,4.2 Kink sector 



In the nontrivial sector the quadratic part of the action is given by 

1 

2h 



dtdx 7]{x, t)YTo\^KV{ x J t) 



The operator jtx\4> k which must be diagonalized is the stability operator. Its spectrum is given 
in the appendix ( pi 4 |315|) . The mode energies are given as in the vacuum ( |182| ) by 



E K {k) = -yUKik) = ~^f^+P (188) 

The difference to the vacuum is given by the different eigen and continuum states which will 
lead to a spectral density which differs from that in the vacuum sector. The spectral density for 
the continuous spectrum of a differential operator relative to its "free" part, i. e. the associated 
vacuum operator is given by the diagonal elements of their density matrices ( [^4] ,po|): 

/oo 
dz [C K (k, z)£ K (k, z) - C v (k, z)Zv(k, z)\ (189) 
-oo 

This spectral density without cutoff regularization is usually used for zeta-function regular- 
ization |^2|. Setting up the same strict cutoff in both sectors, i.e. multiplying (|189|) with a 



common step-function, one obtains the same spectral density as in |39j 



p com (A;) = e(A-|A;|)^ (190) 
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Usual this spectral density is calculated by starting from a mode regularization (see ||39| and 



references therein). Also in [p4j| this detour has been taken. This is the reason why one 
gets problems with the BC and the branch cut position, since the modifications in j34|] only 
lead to a certain mode number cutoff scheme (certain BC and branch cut position). Other 
combinations cannot be produced by the Casimir trick. 

In the next sections we calculate the correction to ( |190| ) appropriate for an EMC scheme. 



First we construct the analogue to the regularized unit ( |185|) in the kink sector assuming a 
different cutoff Ak- The requirement of a consistent regularization will define Ak as a function 
of the given cutoff A, defined by the vacuum. The consistency of the regularization is given 
by the requirement that the path integration over quantum fluctuations must be restricted in 
both sectors to the same subset PID reg ("regularized path integration domain"). To make this 
notion more concrete is the subject of the following sections. 

4.4.3 Sine-Gordon 

For the sine-Gordon model the normalized eigen states are given by 

& = Tidin < 191 > 

e T tm ^:t k (192) 

The regularized unit is given as 



2tt VWTi 



f . ,. 1 1 f , r ,xe ifc(z 2,) (tanh z' + ik) (tanh z - ik) 



Because the eigen-values (|314|) are symmetric (degenerated) in k we made a symmetric ansatz 
for the cutoff Ak, i-e. we choose the same for positive and negative momenta. For the 
calculation of the spectral density we only need the diagonal elements of the density matrix 
z')^(k, z). Therefore it is sufficient to know A K at the point z' — > z: 

i , a 1 f AK dk tanh 2 z + k 2 ^ . 

The integration is easily performed and one obtains 



x ( \ Ak l f 7i - 2arctanA^ ^ A K 1 ■ sf . , , . 

8 K (z, z) = — + — -2- — = + 7^5{A K ) (195) 

7T COSh Z \ 27T / 7T Z7T COSfi Z 



It is interesting that the scattering phase arises here. In the discrete calculation the scattering 
phase 5 comes in by setting up boundary conditions on the asymptotic states. And indeed 
as shown in |39| the shift in the cutoff is connected with the scattering phase. Note that the 
function 5(Ak) = it — 2 arctan A^ is uniquely given by the integration ( [194j ). Thus there is no 
ambiguity in choosing a certain branch cut position. 

To determine the cutoff A^ we require that the regularized subset PID Reg of paths in the 
kink sector is in a certain sense the same as in the vacuum. To establish this we require 
that the two projectors 5y and 5k must coincide (in the sense of distributions). Since for the 
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spectral density we need only the diagonal elements it is sufficient to use (|195|) and (|18o|) . 



The requirement of the equivalence of the off diagonal elements can be understood as an 
equivalence of the two operators in higher orders of O(j), which however is not needed for 
our purpose. Thus we get 

5 K = 5 V ^A K + I / ( ^ } =A (196) 
cosh z l 2 

This is an implicitly given function for Ak since the cutoff A is given and defined by the vacuum 
and renormalization. Therefore we have to solve (|196| ) for A^- This we do by application of 
the Banach fixed point theorem in an iteration up to sufficient order in A. After the first 
iteration we are at sufficient order 

A '- A -2^k 5(A)+0(l di7i 12 ^ ) (197) 

In (|197p both factors - and ^ are smaller than one (if A > 1) and therefore one has 
independently of z a contraction so that the Banach fixed point theorem is applicable and 
(|197|) is a reasonable approximation. In the limit A — ► oo the unit defined in the kink sector 



(|193|) and the kink cutoff A K ( |197| ) converge to the vacuum quantities as 

lim A K = A + 0(|) 
lim 8k(z, z) — > 5{z — z) 

A^oo 

Thus the kink cutoff A^ approaches A as 4- and therefore this difference can never be ne- 
glected (see section [44]) ■ The second relation is nothing else than the completeness relation 



for the spectrum of the self adjoint stability operator. Relation (|196|) is the analogue of the 
requirement of equal number of modes in the discrete case. This ensures that the considered 
subset of paths has in both cases the same dimension. 

We are now in the position to calculate the cutoff-regularized spectral densities. First we 
define some notational abbreviation to keep the calculations readable. 

e A (|*|):=e(A-|*|) 

5 := 6(A) 

0^(1*1) := e (a 1*1 

\ cosh z 2 
© ( A-i 5) A)(W) = 1 for k e (A - Is, A) ; else 

Respecting the different cutoffs one gets the spectral density in an analogous way to ( |189|) as 



p(k) = J dz [QK K {\k\)G&K -e A (|fc|)ft£v] (M) (198) 

= e A (|Aj|) J dz[&tK-$tv} + 1 dz[e AK (\k\)-e A (\k\)]CA (199) 

=: p com (k) + Ap(k), (200) 

where we have split the spectral density into the conventional part obtained by a strict common 
cut off p com and the correction Ap. The integration for p com is elementary and give the well 
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known result 



b a( |/,|) / ji[ ta ^^ + fc2 -i] 



©a(I^I) 



2tt L jfc 2 + 1 
2 



2vr A; 2 + 1 



e A (|fc|)— *'(*). 



(201) 
(202) 



To show that the spectral density p com is not correct at the order of interest we integrate fl202|) 
over k, which gives 



dkp c 



-— arctan A 

7T 



■1 + -T + °(t5) 
7rA v A 2 ' 



(203) 



In the last step we have expanded arctan A around A = oo. From ( |203| ) one can see that the 
conventional spectral density gives the correct (negative) number of discrete states (1 in the 
sine-Gordon model) only up to errors of order O(j). Next we calculate the correction Ap to 
the spectral density p com . It is symmetric in z and therefore we can restrict the considerations 
to positive z with a factor two. With (|192|) , (|200|) one obtains 



A P (k) = - r°dz[e AK (\k\)-e A (\k\)] 

* Jo 



tanh 2 z + k 2 
k 2 + 1 



(204) 



With a variable transformation x = cosh z (|204Q can be written as 

1 . .... „ 1 f\'~ 



Ap(k) 



1 1 

nk 2 + 1 J x 



dx 



|fc|) - ©(A - |fc|) 



x- 



+ 



(205) 



Because of the step functions in (|205|) the integral is only unequal zero if \k\ < A. In this case 
the first bracket in ( |205|) takes the values 



0...A-|A;|- — >0 



-1...A- \k\ 



2x 2 



< 



Therefore the integration in (|205|) is restricted to values 

gA) 
2(A-\k\) 



x 2 < 



(206) 
(207) 

(208) 



The smallest possible value of a; is the lower integration boundary in (|205|) , i.e. x — 1. This 
leads to a constraint for the possible k values (else Ap = 0) 



1 < 



5(A) 



\k\ > A 



6(A) 



2(A-\k\) 

Therefore the correction is only nonzero if \k\ E (A — -|, A). With ( |208|) and (|207| ) we get 

W=W) . ( \/x 2 - 1 k 2 



Ap(k) = -0 



(A-S.A) 



(1*1) 



1 1 

7T k 2 + 



1 



dx 



x z 



(209) 
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Figure 12: The spectral density correction Ap(k) for the values A = 10, A = 100 and A = 1000 of 
the dimensionless vacuum cutoff. 



The integration in (|209|) can be carried out exactly. The final result for the correction is using 
the abbreviation a A (k) :— '' lAl ■• 



2(A-|fc|)" 



Ap(k) 



e (A _| )A) (|A;|)i (\n[^a A (k) + y/a A (k) - 1] 



fc 2 + l' 



q;a(/c) 



(210) 



Although this expression is rather complicated it is not exact. When solving the implicit given 
function for Ak ( |196| ) after A we had to approximate the solution up to order O(-p-) (|197|) . 
The correction is a smooth function of k and its graph is shown in fig. |l2| for different cutoffs 
A. Let us now verify that the approximate solution (|197Q is sufficiently accurate. As a test 
we show that the integral of the spectral density gives the correct number of bound states at 
least to the sufficient order O(t^): 



dkAp(k) 



2 / dkAp(k) 
'o 



~~J dk i^-[\/ak{k) + a/ a A (k) - 1] - fe2 ^ ^ , 



1- 



a A (k) 



(211) 
(212) 



The logarithm can be integrated exactly by a variable transformation z = 2< j^ A ^ and gives 



5(A) 



71 



[ dz[\n(l + Vl - z) - In y/z\ 
Jo 



5(A) 



71 



(213) 



The second term in ( |212|) we can only estimate as follows 



A- 4 



1 



1 



a A (k) it 



< — I dk 



k 2 



A-4 



— [arctanA — arctan(A — -)] = — — + 0( — 

7T 2 71 A 6 A 5 ' 



(214) 



(215) 



Therefore we can neglect this contribution to the spectral density test since it is of order 0(j§) 
and so does not contribute as discussed in section PO. Thus we have for the complete spectral 



density integral the following result (202), (|21 



J dk(p com + Ap) 



2 1 1 , 

•— arctanA (tt — 2 arctanA) + 0(— 

7T 7T A d ' 



-i + o( 



A 3 ' 



(216) 
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Thus the integral of the corrected spectral density gives the correct number of discrete states 
up to the sufficient order O(-h). The remaining error is only due to our effort to express 
Ak through the the regularization/renormalization defining vacuum cutoff A. Therefore we 
had to solve (|196|) in a reasonable approximation. Setting up A K as the fundamental cutoff 
and expressing A as a function of it one obtains a slightly different correction to the spectral 
density p com (p = p com - Ap(kj): 

e v (k,z)£ v (k,z) 
An analogous (and simpler) calculation as above gives 

Ap(k) = ^Q {AKAK+ i s{AK)) (\k\)\n 
The integral of the spectral density gives now 

r 1 1 

/ dkp(k) = — (— 2arctanA^) — 1 H — arctanA^ = —1 

J 27T 7T 

This is indeed the exact (negative) number of discrete states. Thus, a complete matching of 
the diagonal elements of the regularized units gives even an exact result for the sum rule to 
be satisfied by the spectral density. 



Ap(k) 



dz 



0(A 



K 



2 cosh z 



S(A K )-\k\)-Q(A K -\k\) 




4.4.4 Quantum mass of the kink 



The correction Ap results in a additional contribution to quantum mass of the kink compared 
to the mass calculated with the spectral density with a strict common cutoff p com - With ( |183|) 
and (|188|) one obtains 



AM = ^ J dkVk^+lApik) 

= M^r) J A dkVWTi ^[yfcjkj + ^a K (k) - 1] - -j^L^ 



a A (k) 



(217) 
(218) 



The second term in ( |218|) can be estimated as follows 



7T 



dk- 



1 



1 



I hm 

< 



a^ik) it 



dk- 



VW+i 



= (arsinhA — arsinh(A — -) ^-12? 0. 

7T 2 

To integrate the logarithm we transform to the variable z = 2 ( A ~ k ^ . This gives 



hm 5(A) 

'~7T 2~ 



J dz^j (A - 5 -zf + 1 (ln(l + VT^~z) - In v^) 



Expanding the root and with c?2;(ln(l + y/1 — z) — In yfz) = 1 one obtains for the correction 

hm 

(219) 



AM = -frm-U(A)A + O(-J-) - 
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This is exactly the missing contribution in the strict common cutoff calculation |39 
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4.4.5 The large A limit and comparison with the results of [|34 



We want to demonstrate that the correction to the spectral density obtained by which is 
sharply located at \k\ = A can be obtained by a large A limit keeping the necessary orders. 
Therefore one must respect that the correction Ap (|210|) for varying A is a sequence of dis- 



tributions acting on "test-functions" like the mode energy ( |183| )is one. So we investigate the 
action of Ap on test-functions in the large A limit. As test-functions we use smooth symmetric 
functions ip(k) = ip(—k) like the mode energy ( |183|) is one, which grow at most linearly with 
k and for \k\ > A we can assume that they vanish fast enough to be a test-function, since we 
are considering large but still finite A. This symmetry is not really necessary but so we don't 
have to treat the k < 0-domain separately. Thus we consider the following expression for large 
A: 

lAp(<p) = J dkAp(k)ip(k) = 2 dkcp(k)Ap(k) (220) 
First we calculate the logarithmic part of ( |220|) . With (|210| ) and a variable transformation 

2(A-fc) , 

z = - L -^ — - we get 

5(A) / 1 ^(A-^)[ln(l + v / r^)-lnv / ^] (221) 
Jo 2 

J dz[ V (A)-^\A)^z + 0(5\z 2 )}[\n(l + ^/l^~z)~\n^~z] (222) 
S M J(A) + i,<(A)^ + (i)] = -i^Wo^). (223) 



TT 

5(A) 



The action of the second part of ( |210| ) we estimate as follows 

2 f A 1 / 2(A - k) 2. .... f A 1 



A-f A-f 



2 |^(A)l^ + 0(-^) (225) 



7T 



A 3 V A 6 ' 



Therefore get for @2g) with (223),(g2g) 



7T A z 

Thus in the large A limit the distribution Ap approaches 

Ap(k) -> -^5 D (A - \k\) = --5(k)5 D (A - |Jfe|) 

7T 7T 



This is exactly the additional term in the spectral density obtained by |34j using an analogy to 
the Casimir effect. Therefore we understand the results of as a large A limit of the smooth 
correction Ap (|210|) . But the calculation in |34| suffers from the problem that it works only 
for certain branch cut positions of the scattering phase. Other branch cut positions lead to a 
divergent result for the kink mass. This ambiguity cannot be fixed without the use of mode 
number cutoff considerations. In our case the occurrence of the function tt — 2 arctan A, which 
equals the scattering phase 5(A) with a certain branch cut position, in (|195|) is completely 
unique. It simply comes from a uniquely defined integral. 



80 



4.4.6 Robustness of the procedure 



We now investigate the z-dependence and the stability of the relation between the cutoffs A^-, 
A fll 9 71 ) and the spectral density correction Ap (|204T) . For this purpose we approximate the 
factor cosh -2 z by a rectangle of width 2b and high a which is also symmetric around z — 0. 
The area under cosh -2 z is given as J cfccosh -2 z = 2; nevertheless we shall leave the area of 
the rectangle unspecified for now. Therefore ( |197|) changes to 



cosh -2 ^ -> aQ(b- \z\) (226) 
A K = A-aQ{b- \z\)^ (227) 



Here we must require that | < 1 to ensure that (|197Q is still a contraction for finite A and the 



Banach fixed point theorem is still applicable. With this the correction (|2Q4f ) changes to 



Am n [ b dz -ae(b-\z\)+k* + l 

Ap(k) =-e (A -|5,A)(|A;|)y — -^-^ (228) 

= -©(A- W (l*l)i[6-^] (229) 
The integral of the spectral density correction is now given by 

1„ f A „ r , a& 



JdA;Ap ="-2/ ^-^y] (230) 



A-f<5 



1 5 

= ab[S(A) - 2[arctan A - arctan(A - a-)} (231) 

7T 2 

= -I[aW(A)-2o6^ + 0(^)] (232) 

In the last step we have expanded the second term around A = <x>. The first term in the bracket 
corresponds to ( [2131) and the second term to (|215|) . This coincides with the undeformed result 
(|213| ) if and only if ab = 1 and thus when the area of the rectangle lab = 2 equals the area 
under the function cosh -2 z. Analogously, the correction to the quantum mass of the kink 
does not change if ab = 1. The analogue of (|218| ) using ( |229| ) reads: 



AM =-/ dkVWTlAp(k) = -— [ A dk ( bVWTl - ; ab ) (233) 

2 i 7T J \ y/k 2 + 1 / 



hm < b [AVvTT - A J (A - a 6 -) 2 + 1] + ab 6 -^J(A - a^f + 1 (234) 



7T V 2 L V 2' 1 4 

+ -[arcsinhA — arcsinh(A — a—)] — a&[arcsinhA — arcsinh(A — a,—)]) (235) 

a^oo hm ,.11 , hm , 

ab(- + - + + 0) = ab (236) 

7T 2 2 IT 



This gives again the correct result ( |219|) if and only if the area of the rectangle equals the area 
of cosh -2 z, i.e. if ab = 1. It seems to be arbitrary where to locate spatially the modification 
of the kink spectrum as long as the average is the same. Therefore the equality between the 
two projectors Q1961) does not have to be a strong relation, i.e. a identity between operators, 
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but rather a weak relation, i.e. an identity between their action on states, that is, a relation 
between distributions. 

All the above considerations work in full analogy for the 4 model and give also the correct 
results. 

Discussion 

In conclusion we can say that we have a working principle that relates the regularization 
parameter of the nontrivial Ak to that of the trivial sector A in a way that both sectors are 
regularized in a consistent way, i.e. so that one can really compare them with each other. 
This is the main point in our opinion: to find such a relation so that one can regularize the 
nontrivial sector consistent with the vacuum sector, in which also the renormalization is fixed 
and defines the physical parameters (mass, coupling,..) of the theory. Again the big advantage 
of our principle is, besides that it gives the correct results in a consistent way, that it is not 
restricted to two dimensions or supersymmetric theories. It is to expect that it also works for 
fermions. This work is in progress. 

Nevertheless further investigations are in order. Especially the identification of the regularized 
units has to be investigated from a mathematical point of view. For the "off diagonal" term 
these are very complicated integral equations, similar to Fredholm equations, for the unknown 
function A^(A). It would be worth to investigate this further to find out under what conditions 
a solution for A K (A) exists. It is conceivable that in general a strict cutoff- function (step 
function) as used in our ansatz does not solve this problem. 

5 Fermions 

Until now we have only considered bosonic fields and treated them in the path integral as 
classical functions. The nature (physics) of fermionic fields is completely different. While 
classical Bose fields are found in nature (e.g. electromagnetic waves, gravitational fields, etc.) 
classical Fermi fields are not, at least not in the same sense. From the quantum-point of view 
a collection of a very large number of bosons in more or less the same quantum state, i.e. 
a coherent state, can be described by a classical (on-shell) field and also observed as such 
fields, if one does not look too closely. The same cannot happen for fermionic fields since the 
Pauli exclusion principle forbids more than one fermion per state. This "strange" behavior is 
respected in description of "classical" (not operators) fermionic fields by a "strange" algebra 
of the fermionic degrees of freedom, namely Grassmann algebras. Heuristically this can be 
obtained by the formal limit h — > in the anticommutation relations of Dirac fields: 

{ij a (x, t),^(x', t)} = h5 a p5(x - x ) 
{ip a (x,t),if;p(x' ',*)} = {ipl(x,t),ipl(x',t)} = 0. 

By this point of view classical fermionic fields are functions over the space-time, parametrized 
by (x, t) which have their values in a Grassmann algebra. This limit does not describe the 
physical world in an approximative sense as mentioned above. To describe the (quantum) 
dynamics of a system we not only need the functional dependence of the degrees of freedom 
(DOF) on the the space-time parameters but we also need the dependence of dynamical 
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quantities like the action or the Hamiltonian on the fields (coordinates for finite systems). Since 
fermionic DOF are elements of a Grassmann algebra we need a generalization of operations like 
variation and path-integration to "functions" on a Grassmann algebra. For this we consider 
some basic properties. 

5.1 Grassmann calculus 
5.1.1 Grassmann algebras. 

A finite dimensional Grassmann algebra £/at(K) over the field K can be constructed from a 
set of N elements {ai, . . . , ajv}, called generators which fulfill the following algebra: 

{a h aj} = ViJ = 1, . . .,N. 

This relation is invariant under general linear transformations — > G^aj, where the matrix 
entries G^ G K. The whole algebra is a 2 N dimensional vector space in which the ordered 
products 

1 

{ ai \i = l...N} 

{(Haj\i < = 1 . . . N} 

{aiCija k \i < j < k; i, j, k = 1 . . . iV} 

0,\(l2 ■ ■ ■ GtTV, 

form a basis. Concrete realizations of this algebra are for example the exterior algebra of forms 
over a iV-dimensional vector space or the algebra of iV fermionic excitation operators acting 
on a Fock space. The algebra is the direct sum of an even and an odd part, 

Q = Q+®Q-, 

and thus a Grassmann algebra is a Z 2 graduated algebra. The even part consists of all linear 
combinations of basis-elements which consists of an even number of generators and analogously 
the odd part is the linear hull of basis-elements consisting of an odd number of generators. 
Both, Q±, have special features. 

The even part Q + is a commutative sub algebra, i.e. 

for f+,g+eg+ f+g+ G G+ 

[/+,<?+] = 0. 

Therefore one can define functions on Q + and multiply and add them in the usual way. For 
example an element of "degree" two 

N 

W = A,//,r/ r 

i,j=l,i<3 

35 K will be mostly equal to C 
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where the coefficients Ay G K are ordinary numbers and antisymmetric in the indices one 
can define an exponential function, e.g. 

^W m G Q+, (237) 
m 

where W° : = 1 per definition. For a finite dimensional Grassmann algebra this series truncates 
at 2m > N. Also usual function-equations like 

e s e T = e s+T for S,T G (238) 

are meaningful. Only the existence of an inverse element is not guaranteed in general. But 
for the exponential (|237| ) even the inverse exists, i.e. 

e m e -iw = L 

The dynamical quantities like the action (Lagrangian) will be such even functions, so that we 
can work with them in the usual way. 

The odd part Q- whereas is not closed under multiplication (the product of two odd numbers 
can be even) but the elements of Q- are nilpotent, i.e 

ip e G- (ip) 2 = ipip = 




The above considerations are also valid for the infinite dimensional case N — ► oo and this is 
the case of interest for quantized fermionic degrees of freedom. For pure "classical" (on shell) 
considerations one can describe a system of N fermionic DOF within the framework of Grass- 
mann mechanics (pseudo-classical mechanics) as elements of an iV-dimensional Grassmann 
algebra (||Tg|1,||17||). This is not possible for quantum considerations within the path integral, 
where the fermionic DOF do not become operators, even for finite degrees of freedom (not 
field theory) as we will see. To describe the evolution of a system we are of course interested 
in Grassmann valued functions, i.e. objects of the form 

x f(x), 

where B is the parameter domain. With regard to path integral quantization we have already 
chosen an infinite dimensional Grassmann algebra. For a ^-dimensional (pseudo) mechanical 
system thus one has to consider functions of the form 

<i{t) = 

k 

where {f l k } form a complete set in a infinite dimensional function space. The use of the sum is 
a priori a symbolic notation, but will coincide with concrete expressions due to regularization. 
Analogously one obtains for fermionic fields 

ip(x,t) = J~]if) k (x,t)a k , 

k 

where again the the functions {ip^x, t)} form a complete set in an infinite dimensional function 
space. The index set {k} in the case of fields is of course larger than for finite DOF. Of 
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particular interest are spinor fields i/; a (x,t) on a Z)-dimensional Minkowski space, so that the 
underlying function space is the set of square integrable functions ijj : M.d C 2[D/2] in which 
the components {ip a ^(x,t)} form a complete set. The hermitian adjoint field writes as 

k 

where the {a^} are independent of the {a^}, so that the algebra Qoq is generated by the 
infinite set of generators {ak.a^k = 1,...}. By this the fields ip and ift are treated as 
independent degrees of freedom. Usually one decomposes the spinors ipk into positive and 
negative frequency parts as well as according to different spin. This is absorbed in the "master" 
index k. When one studies the Dirac equation in the sense of "first quantization", i.e. a one 
particle wave equation, one solves the Dirac equation for the components ipk(x,t). These 
solution, as for example for hydrogen-like atoms, are in the sense of path integral quantization 
"classical" solutions of the system. Thus "second" quantization is the path integration of 
quantum fluctuations around these classical solutions. 



5.1.2 Variation and integration 

The Lagrangian of a theory is a composite object of the degrees of freedom. In the bosonic 
case thus it can be treated as an ordinary function depending on the fields or coordinates. To 
adopt Lagrangian methods to Grassmann valued fields we have to consider functions defined 
on the Grassmann algebra, i.e. 

L : D{Q) — > I{Q) 

9 -> L(g), 

where D(Q) } 1(G) are the domain and the image, respectively. In general each analytic function 
can be generalized to a super- analytic one ||14|| . For a general Grassmann number of degree N, 
which includes L and g, which are of the form (to make this expansion unique the coefficients 
have to be antisymmetric) 

/ = fo + fi^i + —fijaidj H h —Jix^Oii ■ ■ ■ ai P , (239) 

z. p. 



one can define a norm II II as follows 14 



N 

= i/oi + EE i^i 2 - 

p=lii<..i p 

Thus one has a topology on Q and therefore the concept of "being close to". In the following 
the Lagrangian L will exclusively be bilinear in fermionic DOF. Thus we restrict ourselves to 
the consideration of Lagrangians of the form 

M/cn fa, 9a, 9a) = faD a /3g/3, 

where D a p is a matrix valued differential operator and f a ,g a are tuples (e.g. spinors) of 
Grassmann valued functions of the form 

fa = 5^(/a)iOt , 9a = J^S^A- (240) 
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Thus the Grassmann algebra is generated by the set {a;, bi} where a priori all a«, bi are different. 
Thus under a variation 5f a = the Lagrangian changes by 

i 

5L = 5f a D a pgp = -D a pgp5f a , 

and analogously for a variation of g a , where derivatives must be partially integrated in the 
action as usual. Thus the variational calculus is very similar to bosonic DOF, the only thing 
one has to care for is the ordering of the DOF. One can also define differentiations w.r.t. 
Grassmann variables, which act as derivations on the algebra, but these are rather formal 
operations, i.e. they are defined purely algebraically (see for example [18|1,|^0]]). The reason 
for this is that it is not possible to define differentiation as a limit of differential quotients, 
since the inverse of a Grassmann number, especially for Grassmann numbers like ( |240|) , is in 
general not defined (the inverse exists only if the body fa is unequal zero fi~4| ). 

One can also define a formal integration on the Grassmann algebra, which like differentiation 
is purely algebraic. The so called Berezin integral is defined by the following axioms ( [|i"8H , 
p0[ ): On a N- dimensional Grassmann algebra Q^, generated by {a^}, the linear functional 
/ dai '■ Q+ — ► IK are defined as follows 

1. / dai(l) = 



2. J dai(oi) = 1 

3. {dai, dj} = for i ^ j 

4. {da^ daj} = for i ^ j 

With this definitions one has for example 

day ' ' / dajyf :— / dai ■ ■ ■ da^f = f\ 



,N, 

where /i,...,jv is the "highest" component in the expansion analogous to ( |239| ) of /. Of particular 
interest are integrations of exponential functions of the form 

„ r n 

da x dbi . . . da N db N e~ L XkClkbk = / da x db x . . . da N db N II e ~ Xkahbk (241) 

/ k=l 



n r 

n / da k db k (l - \ k a k b k ) (242) 

NT N 

II A fc / da k a k db k b k = U X k . (243) 

:=1 J k=l 



In the first two lines we have used (|238|) and (|237|) . In the last line we have applied the rules 
1, 3 and 2. 



5.2 The Grassmann oscillator, fermionic boundary conditions 

We consider a hermitian fermionic oscillator. Its Lagrangian is given by 

L — - (a ] a - a ] a) - uja)a (244) 
= -a ] (id t -ui)a + -a(id t +u)al. (245) 
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The fermionic DOF er , a are elements of the infinite dimensional Grassmann algebra which 
is generated by the infinite set {ak,dk}- Thus afl(t),a(t) are of the form 

a(t) = J2fk(t)a k , a\t) = J2f*k(t)^. (246) 

k k 

5.2.1 Variation principle 

There exists a fundamental difference between fermionic and and bosonic DOF. Besides being 
anticommuting fermionic DOF are first order systems, i.e. the e.o.m are first order differential 
equations. In the Lagrangian this is reflected in the fact that the velocities occur linearly. In 
the canonical formalism this leads to so called constraints. For the variation principle this 
results in the need of introducing surface terms for the action to be able to define a consistent 
variation principle ^3|. Variation principles fit perfectly to the principles of quantum theory, 
since they fix initial and final positions (vanishing variation) rather than position and velocity 
at the same time. Thus they lead to boundary value problems rather than to initial value 
problems. But for first order systems this is problematic, since fixing the initial and final 
values overconstrains a first order differential equation. For field systems this generalizes to 
spatial boundaries. Therefore one needs a modified variation principle which is consistent with 
first order systems and leads to the classical e.o.m. This is done by fixing a linear combination 
of the boundary values rather than each of them separately, for example 

a{t') + aa{t") = £ = const => 5[a{t') + aa{t")} = 
a\t') + aa\t") = £ = const => 5{a\t') + aa\t")} = 0, 

where £, $ are constant Grassmann numbers. But for this variation principle to give the 
correct e.o.m. one has to introduce surface terms in the action. For the action 

S= / dtL + aUa{t')a\t") + a\t')a{t")] 
Jt> 2 

this variation principle leads to following e.o.m.: 

5 a) : id — tua = with BC a(t') + aa(t") = £ 
5a: itf + ua^O with BC a)(t') + aa\t") = 

As one can see for the classical "paths" a, a) only the boundary term contributes to the action: 

Sa= l -W)S} + a\t')i}. (247) 
For Dirac fields this generalizes as follows: The variation principle is defined as 

1>\b» - r (M) ^|^ = => 5 l^\ B ,, - r ( ^| B J = o 
^\ B » - r^U; = tl M => 5 l^\ B ,, - r^lsj = o, 

where B^,B^ are the boundaries, spatial or of the time interval, and £b m ,£b are constant 
Grassmann spinors, for the particular boundaries of the "direction" B^. are constant ma- 
trices0. It can be shown that for classical fields, i.e. those fulfilling the classical e.o.m., again 

36 Depending on the explicit form of the Lagrangian the matrices IV^) have to fulfill certaint relations, so 
that the above variation principle give the e.o.m. 
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only boundary terms contribute to the action and this contribution are again proportional to 
the constant spinors, i.e. 



Thus the fields fulfilling "homogeneous" boundary conditions, i.e. £ = = $ give no boundary 
contributions to the action., even if they do not fulfill the e.o.m. For unbounded space-time 
regions one does not have to care about these things, since because of natural boundary 
conditions all contributions vanish at infinity. 

5.2.2 Spectral function and energy spectrum 

As can be shown in the holomorphic representation of the Grassmann algebra of a, a\ the trace 
of the time evolution operator and thus the spectral function is given by the path integral of 
antiperiodic paths Ipifl . Thus the boundary contributions to the action are zero and one 
obtains for the oscillator 



Tr 



e~* HT 



K(T) J Va ] VaeT^T dtL . (248) 



antiperiodic 



The Lagrangian L is given by ( |245| ) and K(T) is an appropriate measure constant. To perform 



the path integration we have to diagonalize the action. For this we determine the coefficient 
functions in ( |246| ) so that they solve the eigenvalue problem^ 

(id t - Lo)f k = t k f k , f k (t + T) = -f k (t). 

The solution is easily obtained and given by 

hit) with p k = (2fc + 1)?r , (249) 



e k =-(pk + cu) , fc = 0,±l,±2,... (250) 



The complex conjugate coefficient function in (|246|) automatically fulfills 



(id t + u)fl = -e k fk, 

where is given by (|250|) . Inserting this into the Lagrangian fl245|) the action can be written 



as 



S= [ dtL = J2 

JT 7, 



k 

For regularization k takes only a finite number of values, i.e. k = —N, . . . , N. The measure is 
defined as 

Va ] Va := K(T) U I dd k da k , 



37 Since the classical solutions do not contribute to the action, this is equivalent to an expansion arround 
around these classical solutions. 
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where K(T) will be defined by a suitable normalization condition. With (|243|) one obtains for 
the trace ( jjgp 



N I % 

K(T) II / da k da k -a k a k e k 



-N 



N i 

k{t) n 



h 

T 



— > lim 



Nh(2k+1)7T 
N —i 

k(t) n 



n (1 uT 

-n \ + (2k + 1)tt 



at ft (2A;+ l)?r 



cos • 



(251) 
(252) 
(253) 



The closed formula for the second product in fl252|) is given in §27§. As in the bosonic case, 
the measure constant K(T) is purely kinetic and does not exist by itself. We normalize it as 
follows 



lim 



A' 



k(t) n - 



T 



Nh(2k + l)n 



Thus we get for the trace 



Tr 



e i 



HT 



E 



e * 



e 2 + e 2 . 



(254) 



The spectrum of the fermionic oscillator consists only of two levels. From ( |254|) one obtains 
the ground state (lowest level) 



En 



and it is same as for an bosonic oscillator but with opposite sign. 
5.3 Mode regularization including fermions 

Now we consider the supersymmetric extension of the SG and 4 model, respectively. As 
mentioned in the introduction we will not stress the supersymmetry of the system (up to some 
fundamental properties) and mainly concentrate on the influence of the nontrivial background 
on fermions in the regularization/renormalization procedure. 

5.3.1 Classical properties 

The supersymmetric extension of Lagrangians of the form 



is given by 



(255) 
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where ip is a Majorana spinor field and V((f>) is related to the original potential U((p) via 
V = 2\^U. For SG and 4 it is given by 

W) = vf " t' (256) 

K scW = ^3i n V^ (257) 

The associated action is invariant under the (rigid) SUSY transformations 

— > + 50 : 5<fi = eip 
i) -> V + <ty : # = - U(0)]e 

where e is a constant Grassmann spinor. The classical equations of motion are 

□0 + U(0)U'(0) + ~WV"{<f>) = 

[i£ - V'(<f>)] $ = 

and there are the following classical (kink) solutions 

fermion vacuum : = 0^ CT -0 = (258) 
fermionic zero — mode : = 0^ CT ip u = —0$ K a Pa e (259) 

where 4>K a are the (anti) kinks ( |148| ). The second solution can be obtained by a SUSY transfor- 
mation of the first one and was first given by f|6]] , |47]] . The projector P a = |[H — ai^ 1 } acts on 
the constant Grassmann spinor e. From this one can see that the ground state {0^,^ = 0} 
is invariant under the half SUSY transformation with parameters P a e = 0. 

For the following calculations we choose a Majorana representation 7 = a 2 7 1 = i(J\ of the 
Clifford algebra. With this choice we have 7* := 7°7 1 = —a%. The intertwiners for spinors are 

$ = = ^ 

The Majorana condition therefore simply becomes ip* — i ) - 



5.3.2 Vacuum sector and renormalization 



The trivial (vacuum) solutions are given as 

SG: i/j v = 
4 : i) V = ® 



rn 



2X 



(260) 
(261) 



Expanding the action around these solutions leads to standard perturbation theory (Feynman 
graphs) and one obtains in a minimal renormalization scheme (m 2 = m\ — (Sm 2 ) susy ) the 
following counter-terms 



SG : (5m 2 ) 



2 " -h 



7m 



susy 



in 



dk 



(8m 2 



susv '"2ttJ VP + m? 



= h- 

90 



dk 



The renormalization conditions are that the bosonic seagull loop for the fermionic two-point 
function vanishes {SG) and the bosonic and fermionic tadpole do not contribute (0 4 ), respec- 
tively. Because of the change of 5m 2 also the counter-term contribution in the kink sector ( |156|) 
is different in the SUSY case. In both cases one has 5M susy = hm J Q A ^-^= := 5M B + 5M F 

SG: SM^-hmJ^^f (262) 
4 : 5M F = -2hm (263) 

The quadratic part of the expanded fermionic Lagrangian is given as 

^ = - V'(4> v )]<il> + 0(Wv) (264) 

In both models one has V'((f>y) = m so that for SG and (f) 4 the quadratic fermionic Lagrangian 
is given by 

Cf = ^[i?-m]i(> (265) 

Boundary conditions 

Applying our symmetry principle on ( |261| ) we get with the ansatz (A = 0, 1) 



M-L/2) = HML/2) (266) 

for the change of the quadratic Lagrangian when transported around the compactified dimen- 
sion 

x — > x — L 5£ = 

Therefore we have to use P/AP BC in the vacuum sector so that the action gets no boundary 
contribution. By contrast, twisted (anti)periodic BC in the vacuum sector as used in p3|| , 
ip(—L/2) = (—) A/ y ir tfj(L/2), induce a boundary contribution to the quadratic action S^ = 
J TxL dtdxC^ of the form ~ J L dxdC = J L dx2rmpip and therefore also to the energy. 



All the following results do not depend on the choice of our symmetry BC (|266|) as it should 



be for BC that do not induce boundary contributions. This is analogous to the bosonic case. 
5.3.3 Spectral function 



In the action expanded around the vacuum ( [26 1| ), S^[i],ipiijj\, no interaction between the 



bosonic and fermionic fluctuations occurs. So we can calculate the fermionic contribution to 
the spectral function separately: 



Tr v e-* H1 \^ = J VW,e*™ + 0(h 2 ) (267) 

antiperiodic 
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Here one has to integrate over fields which are antiperiodic in time |^4|. We diagonalize 
the spatial part of the operator in the action ^[id t + i"i±d x — 7°m]^, which gives (with our 
representation of the 7's) the following spectrum 



ip±(k,x) 



±-(k — im) 



e ikx u± = ±uu = ±Vk 2 + m? 



The BC quantize the momenta as 

ij ± {k,x + L) = {-) A tfj ± {k,x) -> Lk A = {2n + A)>K 
so that a finite (symmetric) expansion gives for the fermionic field 



N 



tp(x,t) = ^ (a n (t)tp+(k n ,x) + b n (t)ip_(k n ,x)) 

—(N—A) 



(268) 



(269) 



(270) 



which now automatically fulfills the BC ( |266| ) , i.e. if)(x + L,i) = (—) A ip(x,t). The time 
dependent coefficients in ( |270| ) are Grassmann-valued functions. 



Majorana condition 

Now we have to set up the Majorana condition if)*(x,t) = i/)(x,t). From fl268|) and (|269|) we 
see that 



ip*±{k n ,x) = ip T {-k n ,x) 
_h.A _ uA 



So the Majorana condition for (|270|) gives 



<(0 = b-( n+A ){t) 
K^) = a-(n+A)(t) 

These two conditions are compatible and therefore the Majorana condition for the field is 
fulfilled for all times 

N 

4>(x,t) = { a n(t)lp+(kn,x) +<(*)^+(fcn,aO) 

-{N-A) 

Inserting this field in the action in (|267| ) one obtains 



N 



o(2) 



-(N-A) 



dt 



~(a*d n + a n a* n ) + u{k A )a* n a n 



This is the sum of 2N + 1 + A Grassmann-oscillators with the frequencies Up = uj{k A ). Note 
that a complex conjugated pair forms one degree of freedom. 
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Ground-state energy 



The measure in (|267|) is 

VipVip = (K(T)) n Tlda* n da n 

and as in the bosonic case independent of the considered topological sector. And also as in 
the bosonic case there exists a subtlety due to zero modes in the nontrivial sector. But now 
it weighs much more as in the bosonic case (see below). Performing the path integration one 
can read off the ground state energy of the vacuum: 




(2n + A)tc 



2 



ground — state energy : Ey = — — y( — j +1 (271) 

mode number : # y = 2iV + 1 + A (272) 

energy cutoff : A A = k& = ^ ( 273 ) 

So one has up to the sign the same ground state energy as in the bosonic vacuum sector, as 
it would be expected by supersymmetry. 



5.3.4 Kink sector 



The treatment of the kink sector is analogous to the vacuum sector, but more involved and 
with some additional subtleties. For the semi-classical calculation (one loop) we expand the 
Lagrangian around the stable kink ground state^ 

{0 = 0^ = 0} 



The influence of the kink in the quadratic Lagrangian (g6J ) is thus given by V'((J)k) which 
reads in our dimensionless variables z = for both models (/ = 1, 2 for SG, </> 4 ) 



Vi((j) K ) = mtanhz 
Therefore the fermionic quadratic action is given by 



o(2) 



- / dt — / dzij)\z, t)[id t + ry*<9 2 — ^mtanh. z]tp(z, t) 

Z It Tfl 



(274) 



(275) 



To perform the path integration fl267|) , but now for the kink sector, we diagonalize the spatial 
part of ( |275| ) , where we normalize the eigen-states properly, so that the factor j- in ( |275| ) is 
canceled. The eigen value-problem to be solved is thus 





-if* 









T 




T 




= OJ 




. p _ 




. p _ 



(276) 



with the operators 



A] 



d z + I tanh z 
—d z + I tanh z 



38 If one expands the action around the other classical configuration where tp = ip zero ^ the fermionic and 
bosonic fluctuations interact already in the quadratic action. 
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The operator in Q276|) is self adjoint w.r.t. the scalar product 



< x\i> >= i dzxty 



if the surface term iXi^i + X^^li^ - (XiV'i + X2' ? /'2)_l/2 vanishes. This is true for the following 
spin structures (anti/periodic and twisted anti/periodic, A = 0, 1) 



(P/AP) : ^(-L/2) = (-)ML/2) ~ 
(TP/TAP) : V("£/2) = H A lML/2) 



p 
p 



T 

b 

2' 



p 



b 

2' 



P 



b 

'2' 



The extra minus for (TP/TAP) is due to our metric signature (+, — ). The coupled system 
of differential equations ( |276|) can be decoupled by expressing the lower component p through 
the upper component £. So one obtains 

AjM = =■ E(i)Z -> eigen values u± =: ±uo = ±j^E^ (277) 



P± 



±p = ±( 



(0 



--M) iff E^O 



(278) 



The differential equation for £ is the same as for the bosonic fluctuations and the lower compo- 
nent p is algebraically related to and thus uniquely determined by £. The case E — (the zero 
mode) must be investigated separately and is given by the classical solution (|259|) . Therefore 
in both models the eigen-functions come in pairs (except for the zero mode) 



± 



m r— 

with u>± = ±—^/E 



(0 



For periodic boundary condition it is necessary to change the eigen-basis and work with "parity'" 
eigen-functions 



^±o, ^ u = (^=Ae«)iff^/o 



which form the pairs of solutions 

u± = 



a 9 

±p9 



With U ± = ±fy/Efi 

with u± = ±^^EJ) 



The explicit expressions for ifj± and u±,(j)± are given in the appendix and will be needed to 
determine the scattering phases. The discrete states fall off fast enough so that they fit in all 
considered spin structures. 



Boundary conditions and symmetry principle 

The influence of the kink background on fermionic fluctuation V'{4>k) is now m both cases 
an antisymmetric function fl274j) and lives therefore on a line bundle with the topology of a 
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Mobius-strip (like the kink itself). Thus one obtains with an ansatz for ip for a surrounding 
of the compactified dimension 

x^x + L: V'{(j> K ) -> -V'{4>k) 

ip -> (-) A Tip 

For the different spin structures the Lagrangian, when transported around the compactified 
dimension, gets the additional contributions 

P/AP : T = 1 5C = V'{(j) K )i)il) (279) 
TP/TAP : T = 7* 5C = (280) 

So for the TP/TAP - spin structures one does not pick up a BC-contribution to the action 
integral. This reflects the residual chiral symmetry of the quadratic part of the expanded 
action. 



5.3.5 Sine Gordon 



First we consider the SG- model with the TAP- spin structure 



fa 

V>2 



-L/2) 



fa 
fa 



(L/2) 



(281) 



(because of the extra minus due to our metric convention). But as we will see the TP- spin 
structure is automatically also treated by TPA. By setting up the BC for the modes ip± the 
whole field, expanded according to this modes, automatically fulfills the BC. 

For the ip + modes (see appendix) one gets for the two components in ( |281|) 



le 



;[kL-\ 



i[kL+6+] 



k + L + 0" 
k + L + ^ 



71 

2nn + - 
2 

2n'vr 

2 



(282) 
(283) 



where ± = arg(±l — ik) are the arguments of the asymptotic (, + (z — > ±oo). One can also 
absorb the factors ±i and the addend ±| in ( |283| ) in the angles 6 I± , but this is pure convention. 
The two quantization conditions (|283|) are consistent if 6 + + | = — {6~ + |) + litm in the 
considered momentum regime. 

The analogous expressions for the ^_ modes are given by 



i[kL+e~ 



-le 



:[kL+e+] 



k'L + Q- -- 
k~L + 6 + -- 

The two quantization conditions are consistent if 9 + — | 



2nir - 
2nn 



2 

TV 



(284) 
(285) 



+ 2 r nm for each k. 



Quantization phase 



We choose 6 + for our quantization conditions. Its graph is given in fig.[L3|. We consider the 
continuous phase (a) for which 6~ is also continuous and the consistence equations for the 
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Figure 13: The quantization phase 9 + for the branch cut positions: (a) cut = Im +7 arg(z) S 
[— 37r/2, 7r/2]. Also plotted are the momenta-evens k£ (solid), k~ (dashed) with positive solution, (b) 
cut = R + , arg(z) G [0,2tt], 



quantization are fulfilled for all momenta simultaneously which makes mode counting much 
simpler. The quantization condition for the modes ip± are 



71 



because of the symmetry of the continuous phase the momenta are related as 

~K = *C„ (286) 

This is the reason why we did not include the | into the phase, otherwise in this relation also 
an index shift occurs. The modes if)± are related to each other by complex conjugation (see 
appendix) so that with fl286|) one has 

So we expand the full fermion field as follows 



M+ 

1p(z, t) = do(t)M*) + Yl M^+On , Z) + <(t)lj)l(k+, Z)) 

fc+>0 

M- 

+ J2 {b n (tW-(K, *) + K(tW-(K,z)) ■ 

fc">0 



(287) 
(288) 



The first term is the zero mode. Due to the Majorana-condition its Grassmann coefficient is 
real. Even if we would complexify ip$ by a complex normalization factor, g?q would depend 
linearly on d . By our choice of the basis and Grassmann valued coefficient functions the 
Majorana condition ip* = ip is automatically fulfilled for all times. 

Inserting ( |288| ) in the quadratic action ( |275| ) one obtains 



o(2) 



dt 



M_ 



M+ 



dodo + 7;( a * n dn + a n d* n ) + tu(k+) (X n (l ri 



This is the sum of Grassmann oscillators except for the zero mode which has to be treated 
by collective coordinates. The path integral measure is up to the zero mode the same as 
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in the vacuum sector if one takes equal numbers of modes. Neglecting for the moment the 
subtleties due to the zero mode one can read off of the spectral function the ground state 
energy contribution due to the fermions: 

h m+ h M- 

E§ = ~ E^-tE^ + * ( 289 ) 



Let us now verify that TP BC give exactly the same result. For TP BC the l.h.s of ( |281|) is 
multiplied with —1. Therefore only the relations fl283p and ( |285| ) are interchanged so that ip + 
has now the momenta from ?/>_ in (|288| ) and vice versa so that in the energy ( |289| ) only the 



names +, — are interchanged. To see this directly was the reason why we used the somewhat 
"lengthy" basis in ( |288|) , with another choice, one sum over positive and negative momenta 
does the same job. 



Mode counting 

If we apply Levinson's theorem to the (continuous) phase one obtains for the number of bound 
states 

»-(-°°)-»-(o,) = i (290) 

Z7T Z 

Thus the fermionic zero mode is a half bound state. Note that this has nothing to do with 
our convention to not include the | in the definition of the scattering phase in (|283|) since this 



constant cancels in the difference in ( |290| ). Therefore the continuous spectrum shifts down 
only by a half mode relative to the vacuum. This cannot be compensated by a discontinuous 
phase, where an integer number of modes does not have a solution. This information can only 
read off the asymptotic values of the scattering phase. For the continuous modes one has 

vacuum : # v = 2N + 1 + A 
kink : # K = M + + M_ + 1 

With the ansatz M + + M„ + 1 = 2N + 1 + A one has to subtract in addition the energy of 
one half high mode (note that the mode energies are negative). The fermionic contribution to 
the kink mass is therefore (using (|271| ) , ( |289| ),( p6^ )) 



M F = E% - El 




J2V(k+y + 1 - ^^fc) 2 + 1 + 5Mp 



2 ^ v v „, 2 

1 

For the calculation of the sums we use exactly the same techniques as in the bosonic case 
(iterative solution for the kink momenta and Euler-MacLaurin). Independently of the splitting 
of the modes M_ and M + and of the vacuum BC (A = 0, 1) one obtains for the fermionic kink 
mass 

M F = ^ + 0(h 2 ) (291) 

Z7T 
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This is the expected and correct result. Note that if we would have counted the zero mode 
as a full mode the result would be divergent. Of course it is unsatisfactory that we have to, 



fall back on the Levinson theorem, like Graham and Jaffe |49|], and cannot produce this as 
an implicit result of mode regularization which validate the Levinson theorem. We think that 
a proper treatment of both zero modes , bosonic and fermionic, with for example collective 
coordinates and the associated path integral measure will give the desired result and validate 
the Levinson theorem. This is work in progress. 



5.3.6 4 -model 



Next we consider the 4 model with the TP- spin structure. For the ip + modes (see appendix) 
one gets for the two components in (|281|) : 

e i[kL+(a+-<p-)] = _i ^ k +l + ( a + _ = (2n + 1)tt (292) 
e <[fcL+( v +-«-)] = _! ^ k +l + fo+ _ a -) = (2n' + 1)tt (293) 

where a 1 * 1 = arg(— k =f i) and <^ ± = arg(2 — k 2 =f 3ik) are the arguments of the asymptotic 
state p + (z — > ±oo) and £ + (z — > ±00). The two quantization conditions (|293|) are consistent 
if (ip + — of) = (a + — (p~) + 27im in the considered momentum regime. Here again we must 
be more careful with the choice of the branch cut position as in the bosonic case. Doing this 
in a consistent way (the same for all angles) one can see that it is not possible to choose the 
same phase at — L/2 for upper and lower components. 

The analogous expressions for the ip- modes are given by 

e i[kL+(a+-<p-)] = 1 ^ k -l + ( a + _ = 2mr ( 2Q4 ) 

The two quantization conditions are consistent if (y? + — a~) = (a + — tp~) + 2ixm. Again the 
the TAP BC only interchange the two relations (|293|) and (|29o|) . 



Phase shift 



There are two branch cut positions for which the consistency of the quantization condition is 
fulfilled for all momenta simultaneously with m = 0, i.e. (f + — a - ) = (a + — f~) =: 5s(k) 
(S stands for SUSY to distinguish between the bosonic phase 5). This makes mode counting 
comfortable. Their graphs are given in fig.0. We choose the discontinuous phase. As one can 
see for both momenta k^ the mode n = has no solution. For k — > 00 63 goes to zero but 
for k — > — 00 5s approaches the value it. It is not possible to choose the branch cut so that 
the phase is zero for k = ±00 as in |Q,||3^|. This can only be achieved if one, inconsistently, 
chooses different branch cuts for the angles a and cp. This finite value at k — > —00 carries 
the information of the half bound state as we will see. Because of the symmetry of the 
discontinuous phase and the quantization condition one has for the momenta 



Lkt + 5 s (k+) = (2n + l)K 
Lk~ + 5s{k~) = 2mr 

/bit" — — k rr, 



(296) 
(297) 
(298) 
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Figure 14: The phase shift 5s and the momenta evens of the quantization conditions (k£ solid and 
k~ dashed lines) for different branch cut positions: (a) cut = R_, arg(z) € [— 7r,7r). (b) cut = R + , 
arg(z) e [0,2tt). 

Proceeding as in the SG case (i()±(k) = ^t( — ^) anc ^ Majorana condition) we get for the full 
quantum field 

M+ 

(299) 

where we have chosen a more convenient representation than in the SG-case since we already 
know that TP - BC give the same result. The quadratic action is again the sum of harmonic 
oscillators and in full analogy to SG (also with respect to the subtleties connected with the 
zero mode) one obtains for the fermionic kink ground state energy 



-Af_,^0 V 



Mode counting 



Applying Levinson's theorem to the continuous phase in fig. 14 one obtains 



<S g (-oo) -fls(-oo) _ -1,1 

p, Trdiscrete ' n 

Z7T I 

Thus again the fermionic zero mode counts as a half mode (bound state), the excited bound 
state on the other hand counts as a full mode. From (|299| ) the difference between these two 
modes becomes clear since the excited modes form a complex conjugated pair of degree of 
freedom (also in the action), in contrast to the zero mode. 

The residual calculations are quite analogous to the S'G-case. For the fermionic 4 - kink mass 
one obtains 

Again this result is obtained independently of the combination of the (allowed) BC for the 
two sectors. 
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5.3.7 Anti/periodic spin structures P/AP 



Since P/AP- BC induce an additional contribution to the action (|279|) they are unacceptable 
for a regularization. But they can be of interest for a nontrivial topology of the universe. As it 
is shown by |50|] a nontrivial topology can violate the CPT symmetry but the derived effects 
vanish in the large limit of the compactified dimension and it is therefore questionable if such 
effects will be measurable. The calculation of the kink mass with P/AP on the other hand 
gives a result that differs from that for TP/TAP BC by a half low-lying mode |39| even in 
the limit L — > oo. The reason for this is that in the P/AP- spin structure the zero mode is 
counted as a full mode. For simplicity we consider the SG- model but all considerations and 
results are analogous for the 4 model. 

Quantization phases. The processing is analogous to the TP/TAP calculation but for 
P/AP BC one needs the parity eigen states u±,<fi± (see appendix). The parity eigen-states 
are of the form 



u± 



t 9 
±p 9 



cu 



where the components are given in the appendix. In u± the upper component £ 9 is an even 
and the lower component p 9 a odd function. For <j>± the situation is reversed. Of interest are 
their asymptotic forms which are given by 



u± 



i 9 {z 
p°(z 
C(z 
p u {z 



±00) = iN q (± sin qz — q cos qz) 
±00) = N q \J q 2 + 1 sin qz 
±00) = N p (± cos pz + p sin pz) 
±00) = —iNp^/p 2 + 1 cos pz. 



Periodic BC gives no constraints for the even components £ 9 and p u . The odd components 
must vanish at ±L/2. Thus we get from the asymptotic states: 

qL 

periodic P : u±(—L/2) = u±(L/2) =>• — = nir, 



M-L/2) = ML/2) 
The second quantization condition can be resolved as 



tan 



pL 



V 



pL + 8p(p) = 2nn with 5p(p) = 2 arctan -. 

P 

Thus the eigen-states u± have freely quantized momenta. For antiperiodic BC the odd com- 
ponents are not constrained but the even components must vanish at ±L/2. This gives 

antiperiodic AP : u±(—L/2) = —u±(L/2) =^ tan = q, 



M-L/2) = ML/2) = 
Again we can resolve the nontrivial quantization condition: 



pL _ (2n + l)ir 
~2~ ~ 2 ' 



qL + 5ap(q) = 2ni\ with 5ap(q) = — 2arctang. 
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Figure 15: The quantization-phases for (anti)periodic BC. (a) For periodic BC the momenta of 4>± 
are quantized w.r.t. to a nontrivial scattering phase 5p. (b) For antiperiodic BC only the momenta 
of u± are quantized in a nontrivial way by Sap- 



For both cases (P/AP) the scattering (quantization) phase is given in fig{L5|. The full Majorana- 
fermion field can be written as 



Ah, 



if)(z, t) = d ip + ^ { a n{t)u + {q n , z) - a* n (t)u* + (q n , z)) (300) 
l 

+Y, {a n (t)4> + (p n , z) + <(t)0>„, z)) . (301) 
i 

Because of the linear dependence and the symmetry in the quantization conditions only the 
positive modes occur in the expansion. The modes n = do not occur since the associated 
momentum has no solution or the eigen-mode is the trivial solution. In both cases Levinson's 
theorem gives one for the number of bound states. But here ( |301|) only half of the continuum 
modes contribute to the bound state. Therefore it is not completely clear how to really count 
this mode. Nevertheless only the naive counting of the zero mode as a full mode gives a finite 
result. This result differs exactly by one half low lying mode from the values calculated above 

(HD m- 

js . r hm hm 
Mp =M d A 1 

This result is obtained independently of the combination of vacuum- and kink- BC, i.e. 
(Vac\Kink) = (P, AP\P, AP). This is evident since the additional contribution to the ac- 
tion (|5.3.2| ) is in both cases (kink P, AP) the same. 

As we have made clear, the origin of this subtlety is tight up with the requirements of a 
proper treatment of both, bosonic and fermionic zero modes. This suggest a connection to 
the Atiyah-Singer index theorem for Dirac-operators in a topologically nontrivial background 
and different topologies of the spin structure. 



5.4 Aspects of derivative regularization 

In this section we repeat the calculation of [38|] for the SUSY-kink mass but in a more general 



way. The derivative regularization scheme developed in |38| is a proper method for the cal- 
culation of the kink mass. We will show now that it is indeed insensitive to the subtleties of 
mode counting encountered above and that it can therefore be used as a benchmark for other 
regularization methods which might be inevitable in the calculation of other quantities than 
the soliton mass. 
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Figure 16: The fermionic scattering phase 5f of ref. [38] with the skipped modes n = —1,0 (TP/TAP 
solid/dashed lines). It corresponds to the branch cut at M + and has a discontinuity at k = 0. The 
phase jumps by 4n so that two modes do not have a solution and fall out. 



Like |38| we shall concentrate on the SUSY-0 4 model. Since the derivative regularization 
involves the derivative w.r.t. the mass parameter m we re-introduce the physical momenta 
which are related to the dimensionless ones as k p h ys = rnk dim ,-i ess . We also make use of 
supersymmetry so that the vacuum energies cancel each other. Thus the SUSY- kink mass is 
given by 

M susy = (E B — E B ) + (E F — E F ) = E B + E F 
For the evaluation of the sums we again use the Euler-MacLaurin formula (|318|) . For ease of 

comparison with m 

we transform the integration variable n to the physical momentum k 



according to the quantization condition (TP/TAP-BC) 

Lk{n) + 8(k(n)) = (2n + A)n 

where the scattering phase 5 is different for bosonic and fermionic fluctuations as will be 
specified below. Sums are thus evaluated through 



n 



/" C dk 

dnf(n)= J ^.[ L + 5'(k)]f(k) (302) 



We have omitted the surface terms since they always cancel each other. The integrand f(k) is 
now, for both fermionic and bosonic fluctuations, given by the derivative of the mode energies 
w.r.t. the mass m. Including the measure one has: 



[L + S'(k)) f(k) = —u K {k) = __ + -VWT^S'(k) + O(-) (303) 
dm v/c 2 + m 2 m L 

In principle the scattering phases for bosonic and fermionic fluctuations must be chosen in a 



consistent way, i.e. one has to select in both cases the same branch cut position. In |38| for the 



bosonic fluctuations phase (c) of fig.|Tl] where chosen, i.e. the branch cut position cut 



For the fermionic scattering phase [g^] used the phase shown in fig.[L6]. For negative k the 



phase differs by 2tt from the phase (b) in fig|L4], which is also obtained for a branch cut at 



The difference comes only from the fact that in |38|] the lower component p of the spinor was 
calculated with the asymptotic values of the upper component £ in (|278| ). Since the difference 



is 2ir this has no physical meaning and the phase used here also gives consistent quantization 
conditions like the phase in fig.|14]. 

The choice of the scattering phase in |3£| thus suggested that the fermionic zero mode has 
to be counted as full mode. We show now that the derivative regularization is completely 
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independent of mode counting arguments and therefore also consistent with the half-counting 
fermionic zero mode. 

For this we leave the magnitude of the discontinuity in the bosonic and fermionic scattering 
phase arbitrary and assume that the modes i/ BF + l, ...vb,f — 1 do not have a solution. Then 
the derivative of the supersymmetric kink mass reads 



dM. 



susy 



d(SM t 



susy i 



dm 



dm 



h 

+ 2 



N B+ 

E 



dujK 
dm 



■ VB- 



Np_ 



E 

,^l/' F + l,..M F - 



duop 
dm 



(304) 



From ( |303| ) one can see that the evaluation of the sums can be split in two parts. The second 
term in (|303|) is independent of L and only involves the derivative 5' of the scattering phase 
which is - for all branch cut positions - the same continuous function. So for this part one can 
perform the limit L — > oo in the integral in (|302|) . Together with the counter-term contribution 
this yields 



d(SM t 



susy ) 



dm 



h t 
h h 

2tt 



+ [ A f v/P+^ (V - 5 F >) 
2m J_ A 2tt 

+ ^ r p- ^ (s B ' ~ 5 F ') 

2m J_ A 27r 



2tt y/k 2 + m 2 



(305) 

h 

•s (306) 



This result is independent of the considered numbers of fermionic and bosonic modes as long 
as the difference between the mode numbers is "small" relative to the highest modes Nb,f ± so 
that for L — > oo this difference vanishes. 



Now we have to show that also the L- proportional part of (|303|) is independent of mode 
counting arguments. First we investigate the required accuracy for our integration boundaries. 
In the following we use the abbreviation A = ^ ZI . Since the integrand is of order O(k) one 
has to respect the integration boundaries in ( |302| ) up to and including order 0(4). Therefore 
we can use the following expressions for the scattering phase- contribution to the boundaries 
(we use the iterative solution for k n ) 



For the different boundaries this gives 



1 2n + A 1 



v, v : 
N b ,f ± « N : 



7T 



x ( 2N+cb,f±+A 



7T 



i5(0 ± ) + O( 

= + o 



L 2 - 



L 2 < 



(307) 
(308) 



In the second line all numbers Nb,f±°£ ( P04j ) can differ by an arbitrary amount cb,f± as long 
as it is not of the order of N, so that —j^- is not of the order A. So the residual integral of 
the difference for bosonic and fermionic contribution is of the form 



A+^&B 



+ / 



-a B 



A^oo 




A+|fe F 



dk 



Lm 



2?r v / fc 2 +m 2 



UPb -a B + a F - (3 F ) 



(309) 
(310) 



The interesting thing is that in the limit A — ► oo the integration boundary-differences jbs,F, 
-tolb,f do not contribute because of the derivative regularization and so the numbers cb,f ± 
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can really be chosen arbitrarily, not only for their phase-contribution (|308|) but also for the 
other part of the integration boundary in (|302|) . Without the additional derivative of the mode 
energies these terms would be of order 0(A) instead 0(j), and there would be potential linear 
divergences which cancel each other exactly only by a correct mode counting, especially the 
half counting of the fermionic zero mode. It is this improvement of convergence by derivative 
regularization that makes this scheme so robust. But for a complete proof of the independence 



of mode counting (not only at the high end) we also have to investigate the contribution ([ST 
further. From fl307| ) and (|302|) we have the following expressions 

a B = {2v B + A B )ir - 5 B (0+) 
a F = (2u F + A f )tt - 5 F (0 + ) 
p B = (2v' B + A B )n - 6 B (0-) 
/3 f = (2v' F + A F )ir - 5 F (Q.) 

Independently of the combination of BC for fermionic (A F ) and bosonic (A B ) fluctuations we 
obtain for (|3"To|) 



, 6 F (0+) - S F (0-) 

K V F — Up) 



2tt 



, 5 B (0+) - 5 B (Q-) 
\y B — v B ) 



2tt 



(311) 



This expression vanishes independently of branch cut positions and discontinuities of the 
scattering phase. This can be seen as follows: The jumps of the phases are always integer 
multiples of 2n, since different (consistent) conventions can only differ by an amount 2tt for 
angles. So the phase terms in ( |311| ) is an integer number which is equal to the number j^ om it 
of the modes which have no solution (for the two phases in fig{L4] it is for the continuous and 
1 for the discontinuous one, and for the phase in fig.|l^ as used by [[38|] it is 2). On the other 
hand the difference of the mode numbers u, v' , which are the first/last modes with a solution 
at the discontinuity of the phase, gives one more than the omitted modes #om%t- Thus for 
(|311| ) one obtains 



vf^omit 1 ^omit) vf^omit ^ Permit) ^ 

So also near zero-momentum the derivative scheme is completely independent of mode counting 
and therefore insensitive to the subtlety connected with the half-bounded fermionic zero mode. 
The only contribution to (|304| ) comes from fl306|) , so that 

dM susy h_ 

dm 2ti 

Integration w.r.t. m gives therefore the correct result for the SUSY-kink mass 

M susy = -t^ + 0{h 2 ) 



where the integration constant is fixed by the normalization M{m — » 0) = ||38 . 



5.5 Discussion 

We have shown that the derivative regularization scheme is very insensitive to subtleties con- 
nected with mode counting. The reason for this is that the potential linear divergence in 
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the difference of the sums in the mode energies is by differentiation converted into vanishing 
contributions. Without the derivative this linear divergence is controlled by correct mode 
counting and the asymptotic values of the scattering phase which is very sensitive on the 
branch cut position. Also the information about half bounded states is encoded in the asymp- 
totic values of the scattering phase and can as such never be seen in derivative regularization, 
because only terms with the derivative of the scattering phases contribute in this scheme. The 
derivative 5' is the same for all branch cut position. As can be seen from ( |306| ) it contains the 
logarithmic divergence which combines and cancels with the counter-term and the supersym- 
metric part of the fermionic and bosonic kink mass contribution. The non-supersymmetric 
contribution which despite SUSY gives a correction to the kink mass is completely given by 
the the derivative of the counter-term contribution. 

Because this scheme is so stable against subtleties connected with mode counting it is a good 
benchmark for regularization schemes that might be required to calculate quantum corrections 
to other quantities then the kink mass. 



6 Conclusion 

We have seen that stable non-trivial classical solution play an important role also in the 
quantized theory. Of special interest are static topological solutions, since they become new 
particle states in the quantized theory. Their stability is guaranteed by the existence of a 
topological conservation law. Although we have considered only two-dimensional theories, 
the conclusions are dimension independent. The existence of a topological conserved current 
is sufficient for the existence of a Hilbert space sector, which is independent of the usual 
vacuum sector. Of course it is much more involved to find non-trivial solutions in higher 
dimensions and as mentioned (Derrick's theorem) it is not possible within a simple scalar field 
theory but one needs gauge fields for a possible existence of topological non-trivial solution. 
Another important feature of solitonic solution is that they are non-perturbative results. That 
is that they are proportional to the inverse coupling constant and thus have an essential 
singularity in the weak coupling limit. From this is clear that they are not traceable within a 
standard perturbation theory. Thus they become important in the strongly coupled regime, 
where standard perturbation theory is not applicable. Especially the possible duality between 
ordinary quanta of the quantum field theory and bound states of solitons make them very 
interesting for perturbative quantum theoretical considerations in the strongly coupled regime. 

The discussion of the last decade has shown that the quantization procedure in the presence 
of a non-trivial background like solitons is a highly non-trivial issue. For the renormaliza- 
tion/regularization procedure one has to compare the trivial and non-trivial sector, respec- 
tively. The crux is to find a regularization scheme so that the two sectors can be compared in 
a consistent manner within the regularization procedure, i.e. one needs a consistency relation 
between the "cutoffs" of the two sectors. We considered mode- energy cutoff- and derivative 
regularization schemes. We excluded dimensional regularization since the existence of our con- 
sidered solitonic solutions depends on the dimension. Also we did not consider zeta function 
regularization, although it is proved to be working perfectly to regularize functional determi- 
nants, but it is incompatible with symmetry transformations including fermionic degrees of 
freedom like supersymmetry |5l] . 



One of the crucial points in mode regularization are the boundary conditions which one sets 
up in the two sectors. As we have shown they can be derived from a very simple symmetry 
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principle. The power of this symmetry principle, besides its simplicity, is its generality. It 
is not restricted to two dimensions or supersymmetric theories. We were able to show that 
all combinations of boundary conditions, allowed by this principle, give the unambiguous and 
correct result for the quantum corrections of the kink masses. To compare the two sectors 
consistently one has to consider an equal number of modes in the both sectors. This is 
equivalent to the requirement that the two regularized Hilbert spaces, or path integration 
domains in functional-language, have the same dimension and are thus isomorphic. There is a 
profound subtlety connected with the occurrence of zero modes and their counting. Especially 
in the case of fermions we could resolve this with an additional knowledge gained by the 
Levinson theorem. This is still an unsatisfactory state within the regularization procedure, 
but the realization of the correct boundary conditions through our symmetry principle and the 
need of half counting of the fermionic zero mode is an important step towards the complete 
resolution of this issue. There seems to be an attractive connection to index theorems of Dirac 
operators in a non-trivial background in connection with non-trivial spin structures. This of 
course has to be investigated further and an important point will be the proper treatment 
of the fermionic zero mode within the framework of collective coordinates and respecting the 
constraints, which are always present for first order systems like fermions. 

For the energy-momentum cutoff scheme we have given a heuristic principle to find a relation 
between the cutoffs in the different sectors. We think that also for the energy-momentum cutoff 
regularization scheme the consistency between the two sectors is given by the equal dimension 
(isomorphy) of the different (regularized) Hilbert spaces or path integration domains, respec- 
tively. It is impossible to set up in a consistent way a common strict cutoff since the presence 
of the kink changes in a nontrivial way the density of states which determine the dimension 
of the space in the continuous case. Our heuristic principle of the "equality" of the regularized 
units is of course only an idea of what really happens. Especially the space dependence of the 
cutoff is somewhat strange. That the procedure nevertheless works well seems to be a result of 
the decoupling of the infrared and UV modes. Also the stability against deformations of the 
space-dependence indicates that the space dependence information is redundant and that a 
more fundamental principle should not include it. Nevertheless it is questionable weather it is 
possible to base the energy-momentum cut off regularization on fundamental principles which 
are completely independent of mode counting arguments, at least for the discrete modes. A 
more constructive approach, relying on the requirement of isomorphic Hilbert spaces, relates 
the kink cutoff to the vacuum cut off by setting the integrals over the difference of the densities 
equal to the number of discrete states, i.e. 



dkp reg (k) = dk dz [9{\ K - \k\)C K (k, z)£ K (k, z) - 6(A V - \k\)C v (k, z)£ v (k, z)\ = 



But in spite of the still open problems in energy-momentum cutoff regularization we have 
showed in an continuum calculation, i.e. independent of any boundary conditions, that there 
is no possibility to set up a common strict cutoff in both sectors. 

In addition we have shown that the derivative regularization scheme, developed in |38|], is 
completely independent of mode counting arguments. Most subtleties in the regularization 
procedure are connected with the potentially linearly divergent Casimir-like energy contribu- 
tions. This is completely circumvented by the derivative regularization. Thus this scheme 
provides a very robust cross check for other schemes and principles. The disadvantage of the 
derivative regularization scheme is that it is not applicable for anomaly considerations and 
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thus, although very comfortable and consistent, one has to use and consistently formulate 
other regularization schemes in addition. 

Summarizing the results of this work we can say that the mode regularization scheme is the 
one with the most advantages. The invention of the symmetry principle to find the correct 
boundary conditions, the independence of the combination of the allowed boundary conditions 
and the realization of the half counting fermionic zero mode are very encouraging results for 
further investigations of the important issue of non-perturbative quantum field theory 



7 Appendix 



7.1 Stability equation 

The Sine-Gordon and 4 kink are members of a family of kinks whose zero mode of the stability 
equation is T}ip(z) = <P'k( z ) ^ eos | 1 i , where I = 1 corresponds to the SG- kink and I = 2 to the 
4 - kink El. The stability equation for this family is given by 



-di 



cosh z 



and can be solved using supersymmetric methods |]42|]. The operator 0\ can be factorized into 
Oi = A\Ai with 

Ai = d z + I tanh z , A\ — —d z + I tanh z 

The discrete spectrum (bound states) consists of the zero mode 

1 



' cosh' z ' 

and a set {m = 1,...,/ — l}of excited states 

1 



E = 



Z,i,m{z)=N Um A\{z)...A\_ m+1 {z) 
The continuous spectrum has the form 



cosh'"" 1 z 



E, 



Lm 



I 2 -(I- m) 



aq,z) = N l (q)Aj(z)Al 1 (z)...A\(z) 



2tt 



E l (q) = q 2 + f 



(312) 



where the factors N are proper normalization constants (f dz^*(q, z)£(q', z) = 5(q — q')). The 
eigen- functions form a complete set. The mode energies of the quantum fluctuations around 
the kink are therefore given as 

u>f = m\[~Ei 
For SG the explicit form of the spectrum is: 



OJ, 



K 




&z) = N 



cosh z 



CO 



K 



(q) = m\j q 2 + 1 £(q, z) = N(q) (tanhz — iq) e 



iqz 



(313) 
(314) 
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For the <jr model the spectrum is: 

1 



lu* = = N — p- (315) 

cosh z 



^ \/3 . , sinhz 
cf = m^- ^=^—5- (316) 
z cosh z 



^( 9 )^^7T «,,,) = W t^,-l (317) 



7.2 Fermionic eigen modes 
Complex waves 



- model u = ^ (^) = ^ ' 



coshz 
Po = 

w±(Jfe) = imVfc 2 + 1 £{(k, z) = N k (tanhz-ik)e ikz 

p{k, z) = -iN k VWTTe ikz 



4 - model uj = £ (z) = 1 



°cosh 2 z 



p Q = 



■ -v/3 t , x , r sinhz 
wi,± = ±m— = Ni — 2 



2 N cosh z 

p 1 (z) = -iAr 1 - 1 ! 



cosh 2 



^±(7,-) = ±/;q/(^) 2 + l i :) .V / ,(3tanlr > : 1 /,- 2 /3/,-tanli :) f ,7 '' : 



/o(ib, 2) = -iN k Vk 2 + 4(tanhz - iA;)e ifcz 
Parity eigen-functions u±, (f>± (continuum states): 

SG — model : u± : i g (k, z) = iN k (tanh z sin kz — k cos kz) . . . even 

p 9 (k,z) = N k Vk 2 + lsin kz . . . odd 
0± : z ) = ^fc(tanh z cos kz + k sin fcz) . . . odd 

p u (k, z) = —iN k Vk 2 + 1 cos kz . . . even 

4 — model : u± : £, g = N k [(3 tanh 2 z — 1 — k 2 ) cos kz + 3A; tanh z sin kz] 

p 9 = —iN k \Jk 2 + 4(tanh z cos kz + k sin fcz) 
0-t : = *-Wfc [(3 tanh 2 z — 1 — k 2 ) sin fcz — 3k tanh z cos fcz] 

p u = N k Vk 2 + 4(tanh z sin kz — k cos fcz) 
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7.3 Euler-MacLaurin formula 

The Euler-MacLaurin formula is given by 

N rN 



with 



£/(n) = J dnf(n) + - (/(«/) + f(N)) + S n 



a ■_ B * ft | B ± f (3) , ■ B 1p f (2p-l) ,N , p 

n '~2! ; 4! ; (2p)! p ' P 

where 5^ are the Bernoulli numbers and the rest is 

i 

The functions Ck(x) are the modified Bernoulli polynomials. 
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